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X-RAY INVESTIGATION OF DISTURBANCE FIELDS DUE TO INDI- 
VIDDUAL DISLOCATIONS IN Si AND Ge 


By J. AULEYTNER 
Institute of Experimental Physies, Warsaw University, Institute of Physies, Polish Academy of Sciences 
(Received August 28, 1960; translated paper received December 6, 1960) 


A rather simple method of X-ray detection of disturbance fields from individual dislocations 
on silicon and germanium crystal surfaces, using a divergent beam from an X-ray point focus 
and a spectrometer with oscillating film parallel to the crystal surface, is described. X-ray mapping 
of the surface and optical mappings are compared. Linear dimensions of the disturbance fields, 
in the case of a silicon single crystal with a density of less than 10* dislocations per cm?, were 
found to range from 20 to 120 u. The increase of intensity of the diffracted X-ray beams within 
the disturbance regions was measured, and a discussion of the contribution therefrom to the 
intensity of the G.-M. counter recorded reflections is given for the case of a double crystal spectro- 


meter. 


The purpose of the present investigation is to show that the disturbance fields due to 
individual dislocations can be detected by a rather simple method, provided dislocation 
density is not too high. Other purposes were: to present an analysis of the diffraction patterns 
obtained, i.a. by comparing them with microscopically obtained pictures of etched single 
crystal surfaces, to investigate the dimensions and density of the disturbance fields, and to 
determine the type of dislocation from the geometry of the latter. 

The linear dimensions of the disturbance fields produced by individual dislocations in 
the atomic lattice of a crystal are known to attain some tens of microns (Lang 1958, Newkirk 
1958, Bonse and Kappler 1958, Borrmann 1958, Barth and Hosemann 1958). Within the 
regions of such disturbance, the X-ray reflection power of the crystal is greater than in non- 
disturbed regions. This makes it possible to use X-ray as one of the ways of detecting indi- 
vidual dislocations. 

Batterman (1958) investigated the intensities and shapes of diflracted beams by Ge 
single crystal surfaces for decreasing dislocation densities and found no variations of mean 
X-ray reflection power below 10° per cm?. Bonse and Kappler (1958), too, used a double 
erystal spectrometer; the diffraction patterns, however, were recorded on a photographic 
film. At the points corresponding to disturbances, they either obtained an increase in 
blackening of the pattern or a weakening thereof (i.e. white-black contrasts), according to 
fine setting of the erystal with respect to the Buonochrometie and almost puo pau 
X-ray beam. The authors explain this effect by assuming a lattice drlorioation to arise in 

‘the dislocations region as a result of local changes in the lattice constants and orientation of 


the lattice planes. 
(371) 
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However. they do not consider the effect of such deformations on the X-ray reflection 
power. The problem seems to be of interest, and hence the next purpose of the present 
investigation will consist in the experimental determination of the ratio pr intense ee the 
beams diffracted within the regions disturbed and non-disturbed by individual dislocations. 

In rendering visible the disturbance fields due to individual dislocations. reference was 
made to observations relating to the explanation of the effect of beam divergence on the 
width of the diffracted line for simultaneous oscillation of the crystal investigated and the 


film (see, Auleytner 1958). 


Fig. 1. 


Consider the two extreme beams OA and OC proceeding from a point focus at O (Fig. 1). 
Assume the crystal at position EC fulfill Bragg’s condition for the beam OA. The reflection 
deriving from this site on the crystal is recorded at point S, on the film. In order that Bragg’s 
condition shall be fulfilled for OC, the crystal has to be rotated by the angle æ. The angles 
OAO’ and OBC’ are now equal, as both are subtended by one and the same are (Bragg’s 
collimation condition). At the same time, the film-holder coupled to the erystal rotates by 
a, and the reflection recorded at S} shifts to the position S’, whereas the one deriving from 
point B is recorded at Sy. As the angles 57,445,, AOB’ and g are equal, the directions AS, 
and BS, are parallel. When placing a plane film coated on one side with emulsion in a po- 
sition parallel to the crystal surface, in the case of an “ideal” erystal, a mapping of the 
surface in the ratio 1:1 would be obtained in the direction perpendicular to the vertical 
rotation axis of the spectrometer support. In the direction of the crystal axis, an extended 
image would result due to divergence of the beam. This deformation should be small if 
the distance between crystal and plate is small and if the plate is set parallel to the crystal 
surface. If the film is relatively close to the crystal, the distance between the Ka, and Ka, 


lines is small, affecting the change in shape of the reflections but insignificantly, since here 
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the angular dispersion is exactly one half that for a resting film. Owing to the two foregoing 
properties, the method is well adapted for disclosing the disturbance fields due to individual 
dislocations and other surface details!. 

The use of characteristic radiation shortens greatly the time of exposure as com- 
pared with other methods, thus e.g. if the characteristic radiation of copper is used. 
In order to reduce the continuous spectrum and f radiation, it is advisable to include a suf- 
ficiently homogeneous absorption filter in the path of the primary beam (preferably made 
of rolled foil). The time of exposure depends on the angular range of crystal oscillation. 
The mapping of large single crystal surfaces requires a considerable angular range of oscil- 
lations and the time of exposure increases correspondingly. During experiments described 
in the present paper, the angular range was 2° when mapping an area of the order of 1 cm? 
and 8° for one of 20 cm?. In order to obtain satisfactory resolving power in disclosing micro- 
defects it is necessary to utilize one-sided fine grained films. Agfa-Dokumentenfilm, Fein- 
korn-Emulsion was found to be the most satisfactory of the ones available. The use of one- 
sided films of small grain-size leads to exposures several times (about 8 x) longer than for 
the usual X-ray films. Hence the importance of using a quasi-point form X-ray source of 
a power of about 5 KW/mm?2?. In the present investigations the X-ray tube with a focus of 
the smallest effective linear dimensions of about 10u was used. The time of exposure ranged 
from 15 minutes to 4 hours according to the reflection power of the crystal. 

In order to obtain a mapping of the surface with beams diffracted from a given family of 
lattice planes, the crystal should be set up on the support so that the direction of 
the normal (n) to the family of planes shall lie in the plane perpendicular to the axis of 
rotation Z of the support (Fig. 2). Such a position of the crystal can be obtained by rotating 
it about the X-axis. If the crystal is now oscillated, one observes a displacement of the dif- ’ 
fraction line parallel to the Z-axis in a direction perpendicular to that of the Z-axis on the 
fluorescent screen. 

The surface of the crystal should coincide with the Z-axis. The position of the crystal 
should be such that Bragg's condition for the wavelength of the characteristic radiation 
utilized in the experiment shall be fulfilled when the spectrometer support is in the position 
corresponding to the centre of oscillations. 

'The film is set parallel to the crystal surface at such a distance that it shall not intersect 
the path of the primary beam as the crystal is oscillated. However, in order to obtain a sharp 
image, the distance between the film and the crystal should be small. 

Fig. 3 shows an X-ray mapping of a germanium crystal surface of block structure, 
as obtained by the method under consideration (magnification 12 x ). 

The broad white line along the image points to the existence of a boundary between two 
regions mutually rotated by an angle of the order of several degrees. The narrower white and 
black lines correspond to low-angle boundaries between blocks of different orientation. 

Fig. 4 shows an X-ray mapping of a Zn single crystal surface with srongly developed 


1 The author is indebted to Dr Barth (Fritz Haber Institut, Berlin) for informing him that Guinier proposed 
a similar method using a divergent beam for investigating block structure in crystals. The present paper, however, 
brings an analysis of the advantages inherent in such a method and discusses the possibilities it provides of 


obtaining images of disturbance fields due to individual dislocations, 
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macro-mosaic and clearly discernable domain structure with preserved privileged directions, 

due to slip. a 
As already mentioned, in order to investigate the resolving power of the method, 
J 'u i ` sili i f a density of 102—10? per cm? 

a mapping of the surface of silicon crystals with etch pits of a density pe 

was obtained. In Figs. 5a and 5b, this mapping is compared with the microscopic image 

(Fig. 5a) of the surface, at the same magnification. Fidelity of the mapping is seen to be 

g. 96 ; 


satisfactory. The light spots on the X-ray image (dark on the positive) correspond to etch 


Z 


Fig. 2. 


pits, as the beams of the rays diflraeted within these regions are partly absorbed (shade 
image). Fig. 5a is a photograph obtained with a microscop of small magnification. The 
pits were slightly over-etched in order that the image should be more distinct at the magnifica- 
tion utilized. 

The foregoing photographs prove the method to possess sufficient resolving power and 
strength for the detection of perturbation fields of linear dimensions of about 10. 

In order to detect the disturbance fields due to individual dislocations, the silicon 
erystal surface from which the microscope image of etch pits had been obtained was polished 
to perfection and etched with a smoothening etching solution (one not yielding etch pits) 
in order to eliminate any scratches. Subsequently, the edges of the crystal were etched 
in a solution yielding etch pits. The central part of the crystal was left unetched by the 
solution, With surfaces thus prepared, an X-ray mapping was obtained (Fig. 6, negative). 
The sites of the etching pits are rendered by light spots on the background of the reflection. 


In addition, the unetched region presents local discontinuous strengthenings of the film 


Fig. 3. Fig. 4, 
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having the same direction of ordering and density. The boundary marked by light spots 
on the X-ray image due to individual etch pits possesses the direction [110]. Since the 
normal to the crystal surface investigated subtends an angle of 12° with the direction [1 11]; 
the dislocation boundary observed most probably constitutes a low-angle lineage of (113) 
planes. These planes comprise the direcions [121] and |211] of edge dislocation lines of 
Burger vectors <101> and <011>, respectively. The unetched parts present ordering of the 
dark spots in two directions, namely, the one coincident with ordering of the light spots 
(the boundary already discussed), and the direction [211]. The latter direction can be inclosed 
in the plane [011]. This plane contains i.a. the following directions of dislocation lines 
that can occur in silicon and may here be taken into account: [011] and [100]. Thus, once 
more, we have ordering of edge dislocations. Screw dislocations inclined at so considerable 
angle with regard to the surface investigated could not possibly yield such a strong effect 
of variation of the diffracted X-ray intensity (Newkirk 1958). 

The spot density in the rows in different directions is the same. Thus, crystallographic 
analysis of the X-ray patterns obtained made it possible (in the present case) to determine 
the type of dislocation locally deforming the crystal surfaces. 

Fig. 7 presents an X-ray mapping of the surface of a silicon single crystal of low disloca- 
tion density. The surface investigated had been polished to perfection and etched with 
a smoothening solution so as to avoid eflects from scratches on the surface. The photograph 
exhibits local strengthenings of blackening and two very poorly visible boundaries 4B and 
CD (whose directions are marked by the arrows on the margin). On etching with a pit pro- 
ducing solution, the surface was X-ray mapped once more. Some of the light spots deriving 
from pits correspond exactly to sites of previous strengthenings of the blackening. However, 
image of etch pits cannot be ascribed to all the sites yielding strengthening of the diffracted 
X-rays. Hence it is to be concluded that the etching solution utilized failed to disclose a certain 
type of disturbance (Fig. 7b — points a, b, c correspond to the same points in photo- 
graph 7a). 

The linear dimensions of the disturbance fields observed in the crystals investigated 
range from 20 to 120g. The sharpness of the dark spots depends on the diffracting plane. 
The sharpest image is provided by reflections from the planes (110) and (112). Obviously, 
these are planes wherein the disturbance fields due to individual dislocations, are large 
enough to yield a distinct effect in the X-ray diffraction image. 

Germanium single crystal surfaces of low dislocation density investigated likewise 
revealed, in addition to point-form perturbations, others having the form of straight lines 
and curves, or even sometimes of loops (Fig. 8a and b, arrows). Investigation by etching 
proved these disturbances not to be due to scratches on the surface. Images of this kind 
are probably due to dislocation lines running almost parallel to the crystal surface. Their 
curvature may be the result of deformations arising in thermal processing o the crystal, 
e.g. in the process of cooling (doubling of the lines arises from resolution of Ka, and Ka, 
doublets). Theoretical evaluation of the increase in reflection power of the crystal at the 
regions disturbed by the individual dislocations would be no easy matter; indeed, it would 


be hardly possible at all in the case of real crystals, just as primary and secondary extinction 
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cannot be evaluated theoretically for real crystals. The foregoing experiments lead to an 
evaluation of the intensity ratio of the diffracted beams at the disturbance and non-disturb- 
ance sites. For the Si erystals investigated, the maximum increase in intensity due to 
a disturbance amounted to 10%. 

It would seem that this increase in blackening cannot be ascribed to an effect of superpo- 
sition of images deriving from the surface of regions of different orientation or diflerent 
lattice constant, The intensity of spots would then be much greater. 

Moreover, the results of measurements of the intensity within the region of spots explain 
why no changes in intensity of the diffraction lines recorded with a double-erystal spectro- 
meter are observed in practice if the sample presents a. dislocation density of the order of 
10* per em? or less. 

The ratio K of the intensity of the beam diffracted on the area S of a crystal surface 
of dislocation density o and the intensity of one diffracted on the surface of a "perfect? 


crystal of analogous field of the face investigated is 


- JS+opl 
T et u OT 
IS 

wherein / denotes the increase in intensity of the beam diffracted at the points where the 
dislocation lines interscet the surface, / —A/,. and p is the area of an individual perturbed 
region. 
IS + opl » S+ opk 

Ig 5 


the quantity A can be determined by photometrization. In the case of the silicon single 
crystals investigated, the maximum value of / amounted to 0.1. If the dislocation density 
is 10* per cm? i.e. 100 dislocations per mm?, and p —0.01 mm?, then k=1.1. Hence, the 
total increase in intensity due to the presence of dislocations amounts to roughly 10% of 
the mean intensity of the reflection at the most. At lower dislocation densities and the usual 
sensitivity of the G.-M. counter devices, it is hardly possible to record the effect of the : 
variations of dislocation density below 10? per cm? on the mean reflective power for X-rays. . 

The author intends to continue the present investigation in order to compare the X-ray - 
and optical mappings for different types of defects and diffusion of admixtures. 


Conclusions 


l. The divergent beam and oseilating film method with quasi-point X-ray source was. 
shown to present the following advantages rendering it applicable in the investigation of | 
the disturbance fields due to individual dislocations: 

a) In the plane perpendicular to the axis of oscillation a 1:1 exact mapping of the single 
(ideal) crystal surface is obtained. 

b) The dispersion of the spectral lines is one half the value found when using the 
stationary film method, providing for relatively sharp mapping of details of the crystal 


surface when utilizing monochromatic Ka. Kæ, Cu radiation and a small distance between 
film and crystal. ) 
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c) The exposure required for mapping the surface of single erystals is considerablv 
shorter than when using other methods. | | 

2. A comparison of the X-ray and optical mappings thus obtained, for the case of 
etched surfaces of Si single crystals, proved the present method to be sufficiently powerful 
for use in detecting individual disturbance fields of linear dimensions of about 10u. 

3. The X-ray mappings of the (111) surface of selected Si single crystals that had been 
polished to perfection and etched with a smoothing solution present local increases in blacken- 
ing of the reflection, having the shape of spots of linear dimensions of from 20 to 100u. 
Some of these spots from rows corresponding to the direction of the low-angle boundaries 
detected by etching. 

4. The density of the spots was found to be approximately equal to that of the etching 
figures, and the positions of some corresponded exactly to those of etching pits, the mapping 
being carried out previous to etching. 

5. From Conclusions 3) and 4), the presence of the spots is considered to be due to 
that of disturbance fields within the crystal lattice as the result of individual dislocations. 

6. In the case of Si, the ratio of the intensity of the beams diffracted within the disturbed 
regions and of the mean reflection intensity was determined, and the local increase in intensity 
due to the presence of disturbance was found to amount to 10% of the mean intensity of 
the reflection at the most. 

7. Proceeding from Conclusion 6), the effect of disturbance fields on the total intensity 
of the reflection recorded with a G.-M. counter or an ionization chamber was discussed 
for the case of dislocation densities of the order of 10* per cm? and less. From these considera- 
tions, changes in dislocation density at values of less than 10? dislocations per cm? are not 
easily recorded with the usual counter devices. This explains the results of Battermann’s 
experiments which revealed no changes in the mean reflective power for X-rays on transi- 
tion to dislocation densities below 10? per cm?. 

The author wishes to express his thanks to Professor Dr L. Sosnowski for his discussion 
of the results and for his kind aid throughout the present investigation. 

The author is indebted to Mr. B. Kolakowski, B. Sc., for preparing the etched crystal 
surfaces, and to Mrs. Grabowska, B. Sc., Mr. Kozielski, B. Sc., and Eng. Kobus for providing 
the silicon and germanium samples. 

The author is also indebted to the technicians Mrs. Z. Furmanik and Mrs. B. Kunicka 


for their valuable help during the investigation. 
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ON THE THEORY OF TRANSPORT PHENOMENA 
IN SEMICONDUCTORS 


By J. KOŁODZIEJCZAK 
Institute of Physics, Polish Academy of Sciences, Warsaw 
(Received October 25, 1960; translated paper received December 17, 1960) 


On the assumption of arbitrary energy band structure, the electric current and heat flux 
are computed. The transport equation is solved by the Mc Clure method. The theory of certain 
thermo- and galvanomagnetie effects is developed on the assumption of a spherical though not 
necessarily parabolie structure of the bands. The computations are for the general case of transport 


by of many kinds of carriers. 


Introduction 


'The more recent methods of measuring directly the effective mass, as e.g. by cyclotron 
resonance or Faraday’s effect, lead to the exact determination of the energy band structure. 
These measurements showed the band structure in real semiconductors to be highly involved 


h 2p2 


In the case of 


and to differ from the simplest case as described by the expression ¢ = ne 
zm 


^ germanium or silicon, the relationship between the energy and the wave vector is quadratic, 
but the isoenergetic surfaces are not spherical. Indium antimonide presents the inverse case, 
namely, spherical energy surfaces but no parabolic dependence on the wave vector. Such complica- 
tion in the energy band structure cannot fail to affect the phenomena relating to carrier transport. 
The theory basing on the assumption of spherical and parabolic bands, if applied to the present 
cases, would lead to considerable error. 
The present paper was conceived when investigating theoretically indium antimonide. 
The general assumptions concerning the structure of the energy bands admit of a high degree 


of generality in the description of transport phenomena. 


$ 1. Solution of the Transport Equation 


In its most general form, the equation of transport is 
De re (1.1) 
dt dt [at \dt f con 
letermining the total variation of the distribution function due to external forces and to 


ollisions. 
The term in eq. (1.1) accounting for the field equals 


s e — — v grad, f — k grad; f (1.2) 


(379) 
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whereas, on assuming elastie collisions; weak external forces, and isotropie scattering, the 
term accounting for scattering becomes 
r 
d — Af 
j D: tsa = ! (1.3) 
T 


dt | con T 
— the distribution function at thermody- 


f denoting the carrier distribution function, f, 
namic equilibrium, and 7 — the relaxation time. Eq. (1.1) is now rewritten as follows 
(1.4) 


) £ 
Ug 


^. df 
Í — foe) + D(is, hy. ke, E) TI 


wherein F is an external force, and v — the velocity of the carrier. 
The distribution function f is expressed by means of the function fọ as follows: 
(1.5) 


— 


In the presence of an external electric and magnetic field, F is the Lorentz force 
(1.6) 


Throughout the present investigation, account will be taken solely of effects linear in the 
electric field strength; thus, the range of electric fields is restricted to that wherein Ohm’s 


law holds. 
Assuming that, in the erystal under consideration, thermodynamic equilibrium obtains 


locally, we have 
grad, f = grad, fo 


Substituting (1.6) and (1.5) in eq. (1.4), retaining only terms linear in the field strength. and 


df, 
Jo grad, T, 


taking account of 
grad, f, = —— 
o ryo dl 


(1.7) 


rad OS oa 
- v grad, T — 0 


ei x T nas 
(ux IH) orad, Or = e. p — — | — 
n T ; dT VT 
Eq. (1.7) will now be solved by Me Clure’s method [1]. Accordingly, the transport equation 


we have 
en 


(1.8) 


(1.7) can be rewritten as follows: 
PELA Pie n. 
dE Le a ara sei) 
with 
E A hea 
Dear 
1 SIVE M 
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The electron velocity v can be expanded in a Fourier series 


b= (nene 


n= —cOo 


(1.9) 
© denoting the cyclotron frequency. The solution of eq. (1.8) is also required to possess 
the form of a Fourier series, thus 


D (mem 
Substituting (1.10) in eq. (1.8, we have 


(1.10) 
cpm T9 Ev (m) — Dv(m) g grad, (1.11) 
1+ imot 
The current density 
E 2e; " 
J= Om): | d,kfv. (1:12) 
by (1.5), (1.9), (1.10), (1.11) and (1.12), assumes the form 


(1.13) 
In deriving eq. (1.13), use was made of the fact that the integrals containing /, vanish 
The vector M is given by the Fourier coefficients v(m) as follows 


x M Yi um) [e;E — D grad, T] 


eiwt(n-+m) 
l + imor 
The vector M is now replaced by its mean value over one period S 


(1.14) 
5 
» ane 
Mg, = s Mdt (1.15) 
0 
The respective value is non-zero for n=—m only. 
Finally, the components of the current density are obtained in the form 
, al A 
Ji = oa E — On (1.16) 
dx} 
wherein gy is the electric conductivity tensor 


O; is the thermoelectric tensor, given by 


dfo zz 
Fi gx aa a Lcd IE 


(1.17) 


(1.18) 
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Tee TEE EE EEE. Ne 3 en 
The tensor Sy is given by the Fourier coefficients in the series expansion of the velocity 


+00 


out y vi(— m) v(m) (1.19) 


1+ imocrt 
m ——0o 


The tensor ©,, is rewritten as the sum of two components 
R 1 


O; = 9; Bo, (1.20) 
with 
2e dfo ER 
94 = : —— Si 1.22 
Vi (22)? fas de T 7 il ( ) 
l "M EA 
1 2: 1:22 
ATAT 2) 
The heat flux is computed similarly: 
QUEE Ü 1.23 
y= Qn dskfev (1.23) 


A computation analogous to the one carried out in the case of the current density yields 


Qi = xaEi— £u A (1.24) 
wherein 
£p = ga Bxa xa = TO, (1.25) 
and 
qiu = ey | dk | — EL T F Si (1.26) 


By imposing appropriate conditions on the current density, heat flux and temperature 
gradient, the solution of the system of equations (1.16) and (1.24) can be made to account 
for all thermo- and galvanomagnetic effects. 


§ 2. Spherical Isoenergetie Surfaces 


If the energy bands present the spherical symmetry, the description of the effects under 
consideration is greatly simplified. Henceforth, the energy will be assumed to depend solely 
on the absolute module value of the wave vector, the form of this dependence being arbitrary. 
Semiconductors exist wherein the energy bands possess spherical structure, or present 
a deviation from sphericity so small as to be negligible. This is the case e.g. in indium anti- 
monide (Kane 1957). The conductivity band therein is described in k-space by spherical 
surfaces; however, the energy-wave vector dependence is not quadratic. 


ees 
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The general assumtpion will be 


Eli ky, ky) = elk) (2.1) 


The velocity of the carrier is determined as follows: 


U; = 


l de 
h ok; a 


by (2.1), this can þe rewritten, thus 


Il aa YE: Jb es 


on ok ahh dh 122) 
With the notation 

JW dé ee 
mE . h?k dk’ (2.3) 

the velocity finally assumes the form 
hk; 24 
v; = 4 
1 m* ( s ) 


Thus, the carrier can be considered to be a free one if the quantity m* is understood to 
represent its mass. In the general case of a band of the type of eq. (2.1), the mass m* need 
not be constant, but may depend on £A. 

We now proceed to compute the components of the tensor Sy on the assumption. of 
an arbitrary spherical band. i 

The coefficients v(m) in the Fourier expansion of the velocity are found by solving 


the equation of motion 


Taking account of (2.4) and of the fact that the absolute value of the wave vector is time- 


independent in the case of motion as described by eq. (2.5), we have 


hk 
ul) = (7-1) = 5 
E _ hk 
pev. neci. (2.6) 
hkg 
v3(0) = m* 


The magnetic field has been assumed parallel to the 03-axis, 4j is the projection of 
the wave vector on the plane perpendicular to the field vector, and kz is its projection on 


the latter. All other coefficients vanish. 
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The components of the tensor S are given as follows: 


2 A 
s < M LUE ent Y 
ae eugiat a 
2 2 
EN h?ksvo 
Te 5 €: P 5 

EL 2 2 \ m* |] 1+ 7202 

1 \2 
; WI 
DLE m* (2.7) 


Hence, in the case of spherical energy surfaces, the longitudinal effects are seen not to 
depend on the magnetic field. By (2.7) and on integrating over the angles, the tensor com- 


ponents of 6, v. and q are obtained in the form 


t 1 €; uj 
dil = O99 FT A EE 5 
= 32? \ 1 + H?u? 


2 
i ed 1 U; 


i 
et er ee 

x we AV H*u? 

i ei 3 
033 = z Qui? (2.8) 
Bi. — pi, ho zu; 

11 s 


22 372 \1 + Hu? 
Hh ziu? 


gt qt 
01» — — 0: 5 
x "27 342 \1+ Hh? 
gt ko 2 9 
O33 = za Qili) (2.9) 
37a 
"que 2 
i i Tks zu; 
( = Go, = L—— | -——S 
ii decr Ste, Ng Hu 
i ; ATR] Zw 
q12 => qz a 2 D 9 9 
712 [21 3a3e; \1 + H2u? 
ba et (d 2.10 
Q33 = y3. (Aili) (2.10) 
OTTO, 
wherein 
oo 
df oi 
(A) = SU) la(s) des, (2.11) 
d de; 
0 
A PN ; eiT; & : é 
u; is the mobility equalling weirs and e; — the carrier charge +e. The index 
m? 
1 0 


"P" denotes a given energy band. Henceforth, the electron charge will be taken negative 
and that of a hole positive: e, — — e, 0, — +e. T is the Kelvin temperature, ky — Boltzmann's 
constant, and H -— the magnetic field strength. Í 
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$ 3. Galvano- and Thermomagnetic Effects 


The number of states per cm? of the crystal is 


1 1 i E - 
hy = UO | | b | k?(e;) de = A [1> (3.1) 
372 ; DIT“ ; 


In considering carrier transport occurring in many bands simultaneously, the respective 
tensors given by (2.8), (2.19) and (2.10) should be added together. 


Thus, the electric conductivity in the absence of a magnetic field can be expressed 


as follows: 
UM gi- | ) ei Cui» ) ein;u (3.2) 
> 71 = iNi Ui 9.4 
1 3m 1 t 


wherein u; is the mobility averaged over the state density. 


u; = > (3.3) 


With respect to the concluding remarks of the preceding Section, e; 


; and u, change their 


sign simultaneously: thus, the conductivities due to different bands are added together in 
all cases. 


A. Isothermal Hall Effect 
The isothermal Hall effect is determined by the following conditions: 


meg n 
dx, E Oka = 


If these are used in solving the system of equations (1.16) with respect to £j, we have 


b= A (3.4) 
OH 
ESRA (3.4a) 
wherein R is Hall’s konstant, defined as 
il Loe 
R : (3.5) 


z i 12 i 12 
H [201 + [3 es 

1 1 
is the electric conductivity in a transversal magnetic field, given by the expression 
i i 12 
[zou] + LE cil 

A i 

i 
È ci 
i 


On 


(3.6) 


OH = 
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By (2.8). Hall’s constant for a vanishing magnetic field assumes the form 


3a* Se; Qui» 


= 3.4 
Ro = "Ii Cun 8.7 

On the other hand, 
Qul» = 3z?n;u?; (uj? = 3a?nju; (3.8) 


A" E 
We have u? — ru; with r; — 
Hence, Hall’s constant is finally obtained in the form 


—2 
>} rj€e;niuj 


i l —2 
a [3 ein;ui]? | g? 25 Suse: G1 


If all species of carriers are non-degenerate, the coefficients r; are equal. Albeit, in the case 


of materials such as e.g. indium antimonide, wherein the effective masses of the electrons 
and holes differ considerably, the electron gas is degenerate whereas the hole gas is clas- 
sical throughout the intrinsic region. Here, the coefficients r; accounting for the electrons 
and holes, respectively, will differ. In the case of a phonon-acoustic mechanism of scattering, 
the difference amounts to 18%, whereas, in that of scattering on ionized admixtures, it 
may amount to as much as 90%. 

For a very strong magnetic field (H— œ), and by (2.8), eq. (3.5) yields the following 
formula for Hall’s constant: 


= [> en]? (3.11) 


B. Magnetoresistivity (Isothermal Effect ) 
The relative variation of the electric resistance in a magnetic field is given by 


Ao O — OH 


me S (3.12) 
with ø denoting the conductivity at zero magnetic field. 
By (3.6) and (2.8), we have 
il ur F 
Ao Xx cia) | «tl (tra) | 
"dat (3.13) 


OH 21 e; Cui» 2: kid —) 


1 


For weak magnetic fields, such that u;H<1, we have the expansion 


Go: 3)^ Qu) — H? <u?) (3.14) 
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The magnetoresistivity now assumes the form 


y M eniu; x einzu -— -I2 eirin; w]? 


g E: 
mE ol 
OH [x e;nju;]? (3.15) 


wherein 


me <u> <li>? „_ <u> <1;> 
en RSEN (3.16) 


For strong magnetic fields, i.e. such that Hu>1, we have 


T [2 eini]? 
-=> =] = - aly 
OH » eini; 2 qiein;/u; 6 s 


with 


1 
(2) Me (3.18) 


C. Isothermal Nernst-Ettingshausen Effect and Thermoelectromotive Force in a Magnetic 


Field 
The effects are described by the following conditions: 
dT 
Ze dx, — 9 


The system of equations (1.16) can be written in the form 
Oy E, + 035 E = Oy, D On = ate Otk 
2 dT 
— djs E, + 035 E, = — Oy, 5— d On = 2: Oir (3.19) 


Solution in E, and E, yields 
3,0530 + ~ O13 2; 01» dT 


T5. ete >= 3.20 
: [x + 35 UE di, ) 
x u x os x 022 oh dT 
Be He (3.21) 
I I» FAI dE ur ole dx 
Let 
Ben (3.22) 


di 
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denote the Nernst-Ettingshausen effect; then the constant Pyg for the case of many species 


of carriers is given by the expression 


m 01» n Du -5 oi » 01, + T DE Has 0 — 5255 E sp 


[es] Dial (3.23) 


2 = = 
Pug = 


For a single species of carriers and vanishing magnetic field. the constant Py, assumes 


the form 


Fee en (u» (zu? d > <su> (3.24) 


deriving the foregoing formula, use has been made of the relationships: 


i €i MAE €i dT T, 
n die: - CER Aon Sw 
The assumption is that — —! is the same for all indices "7". This assumption is by no 
e; € 


means always fulfilled exactly. In consideiing e.g. the different species of carriers to consist 


of electrons and holes, and assuming zero energy at the bottom ot the conduction band, 


the Fermi level for elecirons is 7; = €, = € and for holes is 7; = ^ — — Č —£r, with £4 denot- 

ing the activation energy for thermal transitions from the valence band to that of con- 
: i Vh dG TE e 

duction. If ep depends on the temperature, then the quantities "^ will differ as 


Cj C 
between electrons and holes. 
The thermoelectromotive force is given by the longitudinal electric field of eq. (3.21), 
as follows: 


dV 
Ge (3.26) 
= d E,da, 
Ir a is the same for all indices "7", the longitudinal field (3.21) assumes the form 
€; C » 
1 Qi Či i ] i ? Qi Či i i 
Di VIT D On 2 Era 2 y 1a yy Vio — M Punt dT ldt 
ee ae — i l —' Sas E Me 
[Soi]? + [X ci]? dx, ^ epa. 
Deak 
With Cg 
el da | ] 
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the following expression is obtained for the thermoelectromotive force in a transversal 
magnetic field: 


er ae 3 | 
gm. ; 3.29 
[Y e]? [E oil? Ec 


By eqs. (2.8) and (2.9), and for a vanishing magnetic field, the thermoelectromotive 


l ko ; (zu; &; z 
x (0) = = Duc Zl dme fs (3.30) 


On the other hand, for H— co, we have 


uz koni E R | 


x (oo; SET (3.31) 


force is 


If, in the foregoing considerations, one species of carriers consists of electrons and the 


other of holes, the distribution function fọ in eq. (2.11) is 


k 
Vndc QE ct 
l-4-e^T 
for the electrons, and 
l + 
toe Sag ire for the holes; Mee 
It AF d) 
quite generally, 
E 
Jo = re, f 
ieee (3.33) 


however, assuming zero energy for the bottom of the conduction band, the Fermi level 
for the electrons is ¢,=¢,=—¢ and, for the holes, — £;—0,— —6— er, which, on substitution 
in (3.33), yields eqs. (3.32). 

The author wishes to thank Professor Dr L. Sosnowski for directing the present inves- 


tigation, and for his valuable discussions and critical remarks. 
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The interaction of magnetostriction with internal stresses is one of the principal factors 
influencing the magnetostriction curves of nickel in low magnetic fields. Steinberger (1933) 
found a 10% increase in the magnetic induction flux density at H=0.77 Oe and a pressure of 
10.000 atmospheres in pure annealed nickel. The present authors attempted to establish to what 
extent this effect is due to pressure dependence of magnetostriction. For this purpose, the initial 
magnetization of nickel both under tension and high hydrostatic pressure was measured. 


I. Theory 


In this section we shall attempt to express the initial susceptiblity of strained nickel 


by means of its magnetostriction, tension and crystal anisotropy. 


Stability of the domain magnetization vector in a crystal under tension. Under a tension 
of several kG/mm?, in a well annealed ferromagnetic, the direction of the magnetization 
vector is principally determined by the magnetoelastic energy and the anisotropy energy 
(see E. W. Lee 1955). We shall consider crystals of cubic symmetry only. The lowest 
(nonzero) terms in the expansion of the magnetoelastic energy. consistent with cubic 


symmetry, are 


202 , 202 , „22 o p > > à 
E, = — 3 oroo (Pi + Xab + 2303) 3 0Àq (01 %e/ Da + 06513 + X395) — (T-1) 
wherein: c is the applied (tensional) stress f, Pa, 3 — the direction cosines of the normal 
to the stress surface, with respect to the orthogonal crystal axes, o4, %, & — the direction 
cosines of the domain magnetization vector with respect to the same axes, A499. 444, — Magnet- 
ostriction constant of the material. 


For the anisotropy energy, we have 
7 2 2,2 c 
Ek Kaia, + oS + Aka) (1-2) 


where K, stands for the first magnetic anisotropy constant of the material. 
(391) 
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Using spherical coordinates we have 


Pp, — sin O cos P a, — sind cos G 
f» = sin O sin ® % = sin Ü sin q 
Ps = COS [2] Ae = cos Ü 


The part of the crystal energy of interest here is 
BARE FEK (1-3) 


The conditions for domain stability are: 


IE IE IIE g*p JE TE | 
2 =); ^9 — 0; = 9 = 2 ~ ^ ) >0 (1-4) 
ep Id ep“ DIE Ipod 


The introduction of a magnetic field changes the condition for domain stability. The addi- 
tional magnetic energy is: 


Eg = — HI, cos y (1-5) 


where H is the applied magnetic field, /, — the saturation magnetization at T°K, y 
angle between the direction of the magnetic field and the domain magnetization vector. 
The magnetic field is assumed parallel to the stress. Domain stability conditions are now: 


SQ En) _ 9, 


Ip 


E+ En) 


EIS 


=() “ete, <.. See (1-4) ... (1-6) 


For small values of the magnetic field, 


Ip p“ 
OE 92E E 
on (S2) — Po) + E E AK — bo) (a 


where the values of g, and J are determined from the equlibrium conditions at zero magne- 


tic field. From eqs. (1-6), (1-5) and (1-7) we obtain the expression for the initial susceptibility: 


VEI lI Id cosy 


yen ita os 
I ie p p 9 fe 9 cos y 2 cos y s o?E d cosy : 
90? [o go J MID] ID I 90? |} 9 o9 
8*EY| [IE ag:g \* 
Ip? Jo V90? o. 19989 (I-8) 


On account of the difficulties of exact calculation, approximations are necessary. In the 


case of random distriburion of erystallites in polyerystalline material, the theory of Becker 
and Dóring (1939) involves three assumptions: 
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l. magnetostriction is isotropic, 

2. the applied stress is so considerable that internal stresses and crystal anisotropy 
can be neglected, 

3. only the lowest (nonzero) order terms in the expansion of the free energy in a power 
series of the directional cosines of the magnetization are relevant. 
Becker's theory leads to the expression 


12 
le A et 1-9 
Mo 3(— As) o ( ) 
for the initial susceptibility, where A, is the saturation magnetostriction at T^K. The exper- 


: . LT 
imental data (4) show a linear dependence of upon ug —1-—4 — for 0278 kG/mm? in an- 
o 


nealed samples. 


Some applications of the expression (1-8) 


The temperature dependence of the initial susceptibility of annealed nickel wire under 
tension is, as yet not clearly understood. Taking into account that wires have a pronounced 
texture, and that stresses below yield point are not so large as to justify neglecting the aniso- 
tropy energy, we carried out calculations for nickel, which had been experimentally investi- 
gated by Brockhouse (1953). We used the third assumption of the Becker-Döring theory 
(see above). We take into account two orientations of a crystallite 


1 
y3 


A crystallite is strained along the [111] axis. Let us consider one of these equilibrium direc- 


l. By = By = By = 


7 
tions with 9 = T the values of 9 depends upon ø. Below, we give a table of calculated values 


&—a (which are equal to y upon the foregoing assumptions) for stresses varying from 
10 kG/mm? to 20 kG/mm? in Ni at 0°C. œ is the angle between [111] and [100] axes, equal 
to 54^40' and I, — the domain magnetization vector. In our calculations we put A=5 x 10* 


érg/cm2, 7,—499 (cgs), Ao~ —92 x10 $, Am —19 x10 


o(kG/mm2) 10 ee | 14 | 16 | 1,99 50255) 


9—« 79°A" 80°2° | 80°35’ | 81°18" | 81745" 82°18” 


I 
For equilibrium directions we have: 


IE = sin (à — a) | 3444 sin(9 — a) — 3K cos 20 sin (9 + a) — 5 K sin 29 cos (0 + 2| 


90? 
(I-10) 


4 
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Obviously. 


sin(9—a)~1 cos (O+a)~—1 cos y = T3 [sin (cos @ + sin g) + cos 9| 


d 
sin 202 —1, sin (9 —- a) «£1, an cos y —0, —-cosyzl cos28<1 (I-11) 
dq dü i 
292p | 
It can be shown that — = 0 for all values of 9. Thus, by substituting (I-11) in (1-8), we have 
Ipod : 


TET ET 
Zah, Lok, 


(1-12) 
and substituting the corresponding values for 2,4, and A, for nickel at 0°C we gel 


u 0.23 x+ 0.3 (1-12a) 
X 


where 
c 
TANT kG/mm? 


x 


Brockhouse measured the initial susceptibility of nickel wire which had been subjected 

to a seasoning stress of 19.5 kG/mm?. The texture of this wire is nearly [111] In Fig. 8 of 

the paper by Brockhouse (1953) the initial susceptibility is plotted against the reciprocal 
N 

stress — . We computed 1/z from this graph. The values of 1/x as plotted against ø lie on 
y i 

a staight line. 


The dependence of l/x on g can be expressed as follows: 


Wea UE 

— =0.33%+0.32 at 77°C (1-13) 

x 
This expression is of a form similar to (I-12a). The considerable divergence between the 
values of the coefficients of x in expressions (1-13) and (I-12a) can be explained by the 
departure of the real texture of the wire from an ideal [111] texture, and by a contribution 


from the Z,9 term in the expression for magnetoelastic energy. 
2 m m0 A 
2. Pı=Pp=-) pg 


A crystallite is strained parallel to the edge of an elementary cell. The magnetoelastic 
energy is given by the expression 


Be = 2 Ayo cos? Ü (1-14) 
The equilibrium direction coordinates are: Ü =7/2; d 
A 208 2h 
owe k=; vd; RLE Sn: LE =) 
4 2 op ODOM 
Then, by (1-8), we have 
x = UT 1$ a l f Ni Val 
‘Salto acto creo or Ni at 0°C . (1-15) 


995 


Brockhouse determined the texture of annealed nickel wir es by neutron difraction, finding 
that the texture may be described as a [111] texture (65 +10%) and one in which the [100] 
axis subtends a small angle with the axis of the wire (25 4-109). 


We therefore shall assume that the texture of the wire is a composite one: 65% nearly 
[111] and 35% nearly [100]. Let us consider the expressions (I-12a), (1-13), (1-15). It can 
be seen that expression (I-13) for nearly [111] texture represents an interpolation between 
expressions (1-12a) and (1-15). Similarly, we may express the initial susceptibility for nearly 
[100] texture by 

Hio 7 
0.52x — 0.22 


Thus, for o=10 kG/mm?, we find: 


1 1 


FREE 
E 6.331032 ^ 7 ^ $52 — 032 


— 2.16 (2.62) 


N 


and, for ¢=20 kG/mm?, 


1 = 
== (Oy, /o 0.664082 | 35% 194 — == 1.08 (1.3 1) 


The value of x for [111] texture is taken from expression (1-13). In brackets are given the 
values of x calculated on neglecting the anisotropy energy. It can thus be seen that the 


linear dependence 


i 
ee ee utes 
Ho z 


cannot be considered to represent a decisive argument validating expr. (1-9). 

According to Fig. 5 of the paper by Brockhouse, at 28°C the values of x are 2,25 for 
10 kG/mm? and 1.05 for 20 kG/mm?. It is known that, at 200°C K=O and thus for o= 
—]0 kG/mm? we obtain x= 2.8. Here, we take into account the difference between the 
values of 4,94 at 200°C and 0°C amounting to +15% as found by Corner and Hunt (1955). 

From Brockhouse's Fig. 5 we find x=2.75 at ca 200°C. Our interpretation of Brock- 
house's data can be regarded as satisfactory taking into account the gross approximation 
made. Thus, for the nickel wire used by Brockhouse, we can estimate that the anisotropy 
energy accounts for a ca 15% decrease in the values of the product xo with respect to the 
value for the hypotetical, ideal case, in which there is no anisotropy energy. 

We also computed the initial susceptibility of a crystallite oriented in the direction 


determined by 


f,—0.6821  f,—0.4421 f, — 0.5783 


The angle between the axis of the wire and the [111] axis is now 10°. The calculation shows 
to what extent the influence of the anisotropy energy upon initial susceptibility can depend 
on experimental conditions. Eq. (1-6) was solved by numerical computation. We used 
herein second partial derivatives in the form d by eq. (1-7). This approximation is the 
better, the lesser the values of p—q and 9—d,. We give below a table of the coordinates, 
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total energy and initial susceptibility for possible stabile orientations. There is also a second | 
E x : : Oo | 12059 1o 7a 1 | 
set, of possibile orientations given by equations q'—180*-—g; 0'—180*—0. Values of 


the total energy, cos y and initial susceptibility are the same for both sets. 


Stress fe | " 9 | T | aes iid! | » 
kG/mm? | | | (10? erg/cm?) | 
10 67.25: | 131°34’ 0.1189 | — 0.08 | 5.39 
10 131°29° | 60°52’ 0.1739 + 0.80 | 2.08 
20 | 90°20 | 197° | 0.0018 +13.15 2.66 
20 26°30 | 70°40’ | 0.1477 --16.89 1.96 


Neglecting the anisotropy energy, we get for 


o=10kG/mm2, $g9-100*30; 0 —180'45'; cosy=0.1074 x= 3.94 


II. Experimental Part 


For the purpose of measuring the initial magnetization of niekel under tension and high 
hydrostatic pressure we built a device permitting to strain a nickel element with a force 
practically independent on pressure. A full discussion of the dependence of force and ten- 


sion upon pressure is given in the last section of the present paper. 


Straining appliance and nickel element 


The principle of the straining arrangement is easily understood. A distended 
spring strains the nickel element 4 by means of a block for suspension of the nickel 
element, distance washers and dises (1,2,3). This appliance yields uniform stress in the 
parts of the nickel element having the greatest magnetic resistance. However, it does 
not yield a uniform magnetic field in all parts of the nickel element. We thus are unable 


to make on absolute determination of the value of the initial susceptibility. 


High Pressure Apparatus 


In our experiment, a normal high pressure apparatus with conical external self-support 
was used. The interior diameter of high pressure chamber was 25 mm. In the high pressure 
vessel an additional hole 1.5 mm in diameter was drilled from the manometric to the prin- 
cipal chamber. In the manometrie chamber a manganin pressure gauge for continual use 
was mounted. ‘The bottom of the high pressure chamber was closed by an electrode plug 
with eight conical electrodes. A chromel-alumel thermocouple was introduced directly into 
the chamber by means of two of these eight electrodes. Temperature could be measured 
to 0.05 C. Pressure was measured with to an accuracy of 5 kG/mm?. i 


IAE 
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Measuring Circuit k 


The measuring circuit consisted of three parts: 

l. The field coil circuit including a field coil (200 turns of enamel insulated 0.08 mm 
copper wire) a decade rheostat, milliammeter and current commutator. The field coil pro- 
duced a magnetic field of intensity H ranging 3 to 9 Oe. 

2. The search coil circuit of the nickel element consisting of: the search coil (2000 
turns of 0.05 mm enamel insulated copper wire) and a decade rheostat for adjusting the 
critical régime of a ballistic galvanometer. The sensitivity of the ballistic galvanometer 
used at critical régime was such as to give 5Gs per scale division for our nickel element. 

3. Control circuit including a 10 mH mutual inductance standard. The search coil of 
the standard was in series with the search coil of the nickel element. The control circuit 


was used for eliminating error arising from changes of resistance in the search coil circuit. 


III. Results 


Determining the value of [Mo 


We assume the B(H) dependence in the form B—aH —bH? ... (MI-1) 
according to Raleigh's law for small fields. 
We can write the expression (II-l) in the form: 


x — a'i bi (III-2) 


where a’ represents the sum of four readings of the ballistic galvanometer in scale divisions, 
and i the intensity of the field current. However, (III-2) is not identical with (IL-1) 
because F is not strictly proportional to i. The value of a’ is (with the foregoing restrictions) 
proportional to ug. The value of b’ proved to lie below the threshold of observation. There- 


TABLE 1 


lI. Spring's stretch 3.19 mm 8— ca 9.15 kG/mm? 


l 
fea) "e x(i) — x(50 mA) 50 mA 
Pressure : | T 
75 mA 100 mA 125 mA 150 mA Xy Xy — Xy 
0 414 1.0 0.8 0.7 0.2 411.3 4-0.7 
2500 47.4 0.8 06 ta +0.8 469.7 —0.9 
5000 47.4 0.6 —0.4 +0.1 -1.0 410.8 +0.2 
7500 47.5 —0.6 —0.3 0.0 +0.5 471.2 +0.6 
10000 48.2 2.1 —2.3 —1.5 —2.1 470.4 —0.2 
7500 46.9 0.0 F0.7 +1.4 +23 | 470.5 —0.1 
5000 41.2 0.3 —0.4 +0.6 +1.4 470.7 +0.1 
2500 47.1 0.4 +0.2 +0.8 +18 471.1 +0.5 
0 47.3 0.9 —0.8 "0.3 4-0.5 469.9 —0.7 
—6.7 —4,7 +0.7 +5.7 9470.6 
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TABLE 2 


2. Spring’s stretch 3,65 mm o=ca 10.4 kG/mm? 


an) x(50 mA) ee: 
Pressure — — - 
75 mA 100 mA 125 mA 150 mA xy Xy — Xp 
0 43.1 -0.2 —0.1 +0.5 +1.8 432.9 —3.7 
2500 43.1 | "+04 +0.3 +11 H2.2 434.9 —1.7 
5000 43.7 —0.6 0.0 —0.2 10.6 436.7 +01 
7500 43.9 —0.6 EY, E07 10:3 437.0 40.4 
10000 44.1 —0.6 +0.3 —0.2 —0.3 439.4 +2.8 
7500 43.5 E03 706 | 208 +12 437.2 +0.6 
2800 43.2 4- 0.5 +1.1 +1.8 +2.8 438.2 +1.6 
0 43.6 —0.7 0.0 +0.1 +1.1 436.5 —0.1 
—2.1 +2.0 +3.2 +9.7 9436.6 
TABLE 3 
3. Spring’s stretch 4.40 mm o=ca 12.5 kG/mm? 
50 mA : 
(ata) x(50 mA) en 50 mA 
Pressure 
; 15 mA 100 mA 125 mA 150 mA Xy Kp — iy 
0 36.6 +0.3 — 0.6 + 0.6 +0.2 366.5 —2.8 
2500 36.75 +0.2 — 0.1 0.0 -- 0.8 368.3 —1.0 
5000 37.4 —1.3 —].l E07 —1.5 369.8 +0.5 
7500 31.2 —0.1 —0.4 —0.2 4- 0.2 rules +2.2 
10000 37.4 —1.8 —2.4 —].7 —2.3 365.9 —3.4 
7500 SY? — 0.1 —0.9 —0.5 0.0 371.0 +1.7 
5000 37.0 +0.1 —0.5 +0.2 —0.4 370.2 +0.9 
2500 lew —0.3 — 0.9 + 0.8 —0.7 BY A ies) +2.0 
0 372 —1.0 —0.9 —0.6 —]1.1 369.5 +0.2 
—4.0 —7.8 —2.1 —4.8 9369.3 
fore, instead of x, we used the value 
Xy = Xs9mA T X75mA T X100mA T N125mA + X150mA T € (111-3) 


where & symbolizes a correction term taking into account readings of the control circuit. 
We carried out measurements at pressures of 0, 2500, 5000, 7500, 10000 kG/cm?, in increas- 
ing runs. The measurements were made for the three values of the stress ø of ca 9.15 
10.4 and 12.5 kG/mm?. 

Results of the measurements are shown in Tables 1, II and III. Instead of values of x 
for different values of the magnetizing current intensity, only values of for i=50 mA and 
values of the expression: 

D 


AUN EOS 50mA 


are given. 
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. ~ , 
The data in Tables I. Il and III do not show any systematic dependence of Xy upon 


the pressure. We therefore introduce the average value 


ET 1 a4 RÀ { Tí fi d E. E: 
Aye 8 [Xy, + + AS 9500+ ^ YS 5000 + TAN 7500+ T ^ X 10000 
RE d 1 "Ws E 
+ Ay 75004 + Y55000; + YS 25001 Foul (111-4) 
4 3 , —/ 
The greatest difference between xy and Xy does not exceed 1%. 
IV. Discussion 
Corrections for the straining device take into account: 
a) changes in its dimensions, 
b) changes in spring rigidity. 
The straining appliance is made of three metals: 
a) the spring — of steel, 
b) the nickel element — of nickel. 
c) other parts — of brass, 
It is evident from Fig. 1 that 
Le+h,=ty; and Ald AG — Aly (IV -1) 
where: 
I, is the length of the active part of the spring 
lj, is the length of the active brass parts 
In; is the length of the strained portion of the nickel element. 
Because A/=—Ikp with & for the linear compressibility and p for the pressure, we have: 
, T 
E bye Kr E bor hy) P i UE E lyi byi P (IV -2) 
, "Ve ^ B 
where Al pr represents the additional stretch of the spring, due to the pressure p. 
Thus 
pw i 
p > dk lorkor — Inikni (IV -3) 
Substituting (IV-1) in (IV-3) we have 
Aber 
NES die + lorkor — Inikni (IV-4) 
According to Hooke’s law 
ns Op — 6 
pr p — Oo 
- (IV-5) 


(Alspr)o e 99 


where Al pr stands for the principal stretch of the spring, and 0, — for the stress at pressure 


9° . . . . d B 
p -if the rigidity of the spring does not depend on pressure. Substituting for the quantity 
67,2? r x í 

k”? 1/3 of the value of the well known bulk modulus at O ata values of Loo lor. bi om aa 
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into eqs. (IV-3) and (IV-5), we have, at 10000 ata 
o +o 
nen 05%. 05% 
To 


for stresses of ca 9.15; 10.4 and 12.5 kG/mm?. 


Other corrections 


Corrections for changes in dimension and for pressure reduction of the area of the 
element strained are one order of magnitude lesser than the ones introduced above. However, 
we cannot neglect the changes in spring rigidity. On the ground of Bridgman's data (1949) 
it can be assumed to amount to 2.2% at 10,000 kG/em?. The increase in spring rigidity 
at constant stretch leads to an identical increase in the stress. Taking into account both 
the correction for compressibility and. spring rigidity, we find for the stress of the nickel 
element an increase of 7/10/09 — 1.695 — 1.795. | 


V. Conclusions 
We can put 
xy (o: p) = k(p)f(o) (V-1) 


In the foregoing equation we assumed that /(c) does not depend on p in the first approxima- 


tion. In our experiments we found that 
X: (010000 ata)’ 10000 ata) = xy, (Oo ata); Ô ata) (V -2) 
From 1 and 2 it follows that 
k(10000 ata) x /(o<10000 ata)) = [t(9 ata) X f (o(o nbl 0 ata) 


Ak Af (0g) Ak = k(10000) — (0) 
—— b where j 
k J (Go) Af (o) = f (010000) — Foo) 


If 0,0009 were equal to Og Xy; (10000 ata) would exceed xy (0 ata). 

We can put xy (10000 ata) = xy (Coama; O ata) and Aas = 25(6 ata; 10000 ata) — 
— Xy (oaa; ata) 

Axs Ak = Af(ag) 

xx k f (aq) 

We can obtain the values of the function —df(o)/f(c) (or rather: —A[f (o9)]/f(cg)) from 
Tables I, II and III and eq. (V-1). 


We have 


aS Eus 


For ie — 1.6% — 1.795 we found — 
(B) 


If we put domm -= A “* (we neglect the dependence of H upon x and oc) 


x 
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we have = = 1.2% + 1.3% where Ax = x (ay; 10000 ata) — x (09; 0 ata) 
; A 
Ax = Apo Au Ko since u, e 30 
“ Hol Mo 
Ax Ao 


If we put 


— + — we have 
oO 


ae 16% +1.7% 
x 
Thus, we can conservatively evaluate the increase in the value of xo of nickel at a pressure 
of 10,000 kG/em? under tension as amounting to 1—2%. 

We have here also into account the uncertainty from assuming pressure independence 
of Xy. The change in xo cannot be given a strict interpretation, as the texture of the element 
is unknown (see ““Theory’’...). It is possible that the effect of the increase in xo is due to 
a decrease in the anisotropy energy as suggested by D. Gugan and G. Rowlands (1958) 
who, from data published by Steinberger and others, computed a 10% decrease in the 
value A, for nickel at a pressure of 10,000 ata. 
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Effective formulas determining the electromagnetic field for the case of Fraunhofer 
diffraction within the framework of Kottler, Severin and Vasseur theories are derived. The prop- 
erties of an arbitrary plane diffracting screen are shown to have no effect on the distribution of 
the electromagnetic field energy density near the centre of the diffraction image, in the theories 
considered, as in this case the formulas yield identical results which is not true in other cases. 
The author suggests that measurements should be carried out with the aim of deciding as to 
which of the three theories yields results in the best agreement with reality. 


Hitherto only a restricted number of exact solutions of diffraction problems are known. 
Because of the importance of these problems, which are related to the present rapid develop- 
ments in radar, millimetre and sub-millimetre wave techniques, and with respect to the 
considerable difficulties involved, various approximate electromagnetic theories become 
increasingly useful. Such theories already played a role in classical physical opties, where 
they yielded not only means of verifying known results, but, moreover, were used in analyz- 
ing these features of the scattered radiation which quite naturally lay outside the scope 
of Kirchhoff’s scalar theory. 

The present paper brings a discussion of the approximate electromagnetic theories 
of Kottler (1923), Severin (1946), and Vasseur (1952) together with a comparative study of 
the results yielded in each separate case. 

An infinite plane screen S having an aperture 2’ of contour /' is assumed to be situated 
he plane z=0. The waves incident theoreon come from the half-space z —0. The diffracted 
ield is investigated within the half-space z>0. 

In Kottler’s theory, the diffraction problem is solved not as a boundary one but as 
| saltus problem. In particular, Kottler assumed the field components to have a discontinuity 


saltus) at the surface of the screen S, so that e.g. Bye Easy Eo with E, denoting the 
lectric field of the wave incident from the region z<0. The Kottlerian electromagnetic field 
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in z>0 is given by the formulas (the factor e" is omitted): 
5 Ae! LH EET PT 
ET Ani o BE TI ei. 
= | (5x Ej) Gds + 
Xx. 
| MN " 
+ Trik’ (sH,) grad G ds (la) 
Eat 
- l + OG , 9H, » 
EE m | Ho en E en = 
Mer RER 
ue vue | (SX Ho) Gds — 
v7 
Id bee 
Hur fen, grad G ds (1b) 
J 


wherein G = e*/r, r is the distance of the point of observation from the current 
point of the surface X, and the operator grad acts on the coordinates of the point on 2; 
k=2a//. The other notations are: 5 — the unit vector tangent to the curve J’, and n — that 
of the normal to & ~ directed parallel to the z-axis. The directions of s and n form a right- 
handed system of axes. Àj and li, are the electric and magnetic fields of the incident wave, 
respectively. 

From the point of view of physics, formulas (la, b) yield a deseription of the electro- 
magnetic field arising from the electric and magnetic charges and electric and magnetic 
currents distributed over the screen S (the respective formulas containing integrals over S and 
I can be transformed exactly to the form (la, b)). Kottler's formulas fulfill Maxwell’s equa- 
tions and represent the exact solution of the electromagnetic saltus problem. It should be noted 
that the formulas can be applied to an aperture of arbitrary form, eut in an arbitrary (not 
necessarily plane) screen. The theories proposed by Vasseur and Severin, which are to be 
discussed later, can be applied in the case of the plane problem only, as shown in Fig. 1. 

In Vasseur's theory, the diffraction problem is solved as a ‘one-sided’? boundary 
problem, i.e. in this theory, on the surface z= +0 of the screen S the condition nxE=0 
is fulfilled. Thus, the screen S is a one-sided perfect conductor, as it were. From this point 
of view, the Vasseur type sereen would seem to be physically more realistic than the one 
of Kottler's theory. Vasseur proposed the following formulas for the electromagnetic field 
in the region z>0: 


-= l S M m 

Ey — P" (n x Eg) x grad Gd X, (2a) 
z l ir f y =g >, N E l "m 
Hy = x J (— ik(n x Ey G+ (nH,) grad G} dX — 5x (sE,) grad G ds. (2b) 
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The notations are those already used. The foregoing formulas fulfill Maxwell’s equations. 
From eqs. (2a, b) it will be seen that Vasseur assumed the electromagnetic field on the 
aperture X to be the non-perturbed field of the incident wave. This approximation the more 
justified the greater values ka (ka>1) assumes (a denotes linear dimensions of 2). From 


M 


nI 


the point of view of physics, the field (2a, b) is produced by the magnetic currents distri- 
buted over the surface J. 


Severin proposed the following formulas for the electromagnetic field in the region z>0: 


7 1 SUPRA. Ty +7 il E Ih m d 
Es = Dx | us x Hy] G+ (nE,) grad G} d£ + za] (sH,) grad G ds (3a) 
5 P 
er (Pee RE 
fg | (n x H,) x grad G dd (3b) 


Here, the condition 7 x E=0 is not fulfilled on the screen S; but in conformity to real 
condition the tangential component of the magnetic field on the aperture 2 is equal to that 
of the electromagnetic field of the non-perturbed incident wave. Severin's formulas fulfill 
Maxwell's equations. The field they describe can be interpreted as due to the electric currents 
distributed over the surface JY. 

Formulas (la, b), (2a, b) and (3a, b) are rather involved. Even in the most simple cases, 
they require approximate computation. In order to provide a starting point for an eflective 
comparison of these formulas, we shall now proceed to generalize them to the case of Fraun- 
hofer diffraction. Assume the incident wave to be a plane one and the point of observation 
(in the region z>0) to be remote from the plane z —0. Thus, the problem consists in trans- 
orming formulas (la, b), (2a, b) and (3a, b) in such a way as to obtain new formulas describ- 
ing the field at great distances from the surface S+2. 

Assume a plane wave to be incident on the screen from he half-space z«0: 


By = Ae, 


Hy = (px A) e*Fe (4) 
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with p denoting the unit vector determining the direction of propagation, and A being the 


amplitude vector (in the general case, complex), (4p) —0. Consider formula (2b): 


E meet = ee Y ~ l -p M 
De = I —ik (n x E) G+ (nH,) grad G} dZ— X (sE,)) grad G ds 


9 5 
AT, SITU a 
S T 
eir eikr E 
Transition to Fraunhofer diffraction is carried out as follows. grad G = grad — ~ ik — — 
r rg 


> 


: : br LE 
and, in the distant zone, the vector — is assumed to be equal to the vector & directed towards 
r 


the origin of coordinates (Fig. 2). The quantity r in the exponential is replaced by rọ + Ra, 


with rg denoting the distance of the point of observation from the origin. and R — the 


ly 
id 


radius-vector of the current point on the surface (Y (Ay). The term 1/r is replaced by 1/ 


To. 
Now (by (4) 
| METTUS ME et Be : ; 
ex em ik(n x Ey) G+ (ny) grad Gi dX 
. . a "7 
ik eikro | > " aR > > à E ET => Br. „2 > > l 
ERNE 1$ x | nenn (ote) dd'+ a| (px A) | ne*R Ota gy (5a) 
on ry IO, : | 
— x 
>» A ikra 44 
T: | (SE) grad G ds > — lem. | TE "n (5b) 
2zik . 2% To | K H - 
2 X 


Putting p exp [i R(p- a)| in the relationship / 


gds = [n «grad y d V, we have immediately 
p x 


pa 


ARGS >? PNA A - amc in 
| Mig Bak Uae tk (p+ a) x f nec ra) ID 
P 4 


(5c) 


v 
u 


Eqs. (5a, b, 9) together with (2b) now vield 


Hy — B(Ax Fax A)F] - a[A(( +a) x FY} a) 
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wherein 


= | ien OTT DT ee 


Similarly, we obtain 


Ey = Bia x [Ax F]) (6b) 


Generalizing the formulas of Kottler and Severin for the case of Fraunhofer diffraction, 


we have 


-- & [A [(p 4- &) x FI} (7b) 
BE = B(F x (px A) - a(FA) —a@ (px A) [(p-- a) x F]I). (8) 
HE = B[Fx (px A) x & (8b) 


> > 


The notation is the one used hitherto. The indices A, S, V attached to the vectors E and H 
denote the field of Kottler’s, Severin’s and Vasseur’s theory, respectively. , 
Eqs. (6), (7) and (8) describe transversal electromagnetic fields, since the following 


relationships hold: 


P =) 
foe (9) 
=a Ae 
in the case of each of the fields considered, as is readily verified. 
Moreover, it is obvious that 
iir. E 
Kk $0y t Es) (10) 
x — Ey Hs) 


The existence of a quantitative relationship between Kottler's theory and those of 
Vasseur and Severin was to be expected a priori, considering the sources to which the fields 
of eqs. (6), (7) and (8) are due. 

In the general case, for arbitrary &, all terms in eqs. (6), (7) and (8) are of the same order 
of magnitude, and thus should all be taken into consideration. The situation is entirely 
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different. however, when considering the energy distribution in the proximity of the centre 
of diffraction, i.e. for & e — p. 

Let Gx, Gs, and Gy denote the field energy densities in the cases of the theories 
of Kottler, Vasseur and Severin, respectively. By (6). (7) and (8), taking into account the 
equation. deseribing the energy density 

| 


eo. — A HII* 
4 = qs (EE*+ HH”), 


and with the relationships of eq. (9), which yield EE*=HH*, we derive immediately that, 
for x— —p, 


EIS ee 
K= s= óy = m |4|? |B|? |F B a3 u 


(a, is the z-component of the vector x). From eq. (11) it can be concluded that, in the case 
of Fraunhofer diffraction, the properties of the screen S have no eflect on the electromagnetic 
field energy density distribution at the centre of the diffraction image. Moreover, since the 
fields described by eqs. (6), (7) and (8) are given by continuous functions for all &, the 
foregoing conclusion still holds if x differs but little from — p.. Hence we conclude that 
eqs. (6), (7) and (8) for x & — p conform to reality since, as is readily verified, they yield 
in this case results that coincide with those of Kirchhoft’s scalar theory, which has been 
confirmed experimentally. 

In order to gain insight into the situation as it is when a differs considerably from — p, 
we shall deal with the following simple case. Assume a plane circularly polarized wave to 
be incident perpendicularly on the aperture X. We investigate the energy density distri- 
bution in the plane x—0. Simple computations yield 


o l > 9 Jn 19 Z 

QR = 324a lal“ LB]* 121° (a — 1)*, (12a) 
o o l | 2 PES Am) 2 € 
OS'S Oras ha |A|? |B]? |F|* (a3 + 1) (12b) 


Thus, in the case under consideration, the theories proposed by Vasseur and by Severin 
both lead to the same result which, however, diflers from that of Kottler's theory. It is note- 
worthy, on the other hand, that eq. (12a) is here identical, to a multiplicative constant, 
with the corresponding formula of the scalar theory. Fig. 3 brings a graph of the ratio 


on m" 
Uu Versus ats 


OS 


Ihe function w is seen to diminish monotonously from 1 to 1/2 as æ, varies from —1 (at 
the centre of the diffraction image) to 0. Hence it is to be concluded that the divergence 
between the results of the scalar theory and those of the electromagnetie theories of Vasseur 
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and Severin are the greatest when the vector & is almost perpendicular to p, i.e. at great 
distances from the centre of the image. With the existing possibilities of the measuring 
techniques, the decision as to which of the formulas (6), (7) or (8) yields a description of 
the energy density distribution in the diffracted wave that is closest to reality should come 
from experimental results. 


Fig. 3 


The author thanks Professor A. Rubinowiez for his interest in the subject of the 
present investigation. 
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Introduction 


To F. Kottler (1923) is due the correct formulation of Huyghens’ principle for electro- 
magnetic waves. His theory is the conclusion of the course along which Larmor and Tedone had 
worked towards a solution of the above problem. Kottler's idea was so successful because he 
took into account the electric and magnetic charges on the edge I of the diffracting aperture. 
This was a step forward in comparison with the ideas of Larmor and Tedone, who had assumed 
that the diffracted image was due to charges and currents of the electric and magnetic type distrib- 
uted appropriately on the surface S of the diffracting body. If the surface S is closed, then, of 
course, Kottlers edge effect vanishes and we have the Larmor-Tedone formulae exactly. It is 
now necessary to specify in detail why we think Kottler's formulae to be the correctness. Namely, 
his solution of diffraction problems is at the same time a solution of the complete system of Max- 
well's equations. A theory of diffraction is correct only if its final formulae fulfill this system of 
equations. But this, clearly, is not a sufficient condition for the correctness of the theory. Such 
is also the case of Kottler's formulae. Kottler solved the diffraction problem as a “saltus problem"; 
however, if the problem is to be stated correctly, it must be considered as a problem with boundary 
conditions. Kottler’s formulae can be presented in the form of surface integrals extending over 
the surface of the diffracting aperture and of curvilinear integrals over its edge Z”. In these integrals, 
the nonperturbed electromagnetic field appears and it might seem that the presence of the dia- 
phragm has no meaning as faras the magnitude of the field on the surface of the diffracting aperture 
is concerned. From the physical standpoint this is wrong. Thus, we have here an analogy with 
Kirchhoff’s scalar theory. Nevertheless, it is a well known fact that, on the basis of Kirchhoff's 
theory, we can obtain a correct description of the variability of the light intensity in the neighbour- 
hood of the shadow boundary limit. Thus, we can also conjecture Kottler's theory to be the correct 
one in this region and its results to represent naturally a good approximation consistent with the 
results of Kirchhoff’s theory. If this were not the case, Kottler's formulae would be wrong, since 
they would then fail to describe correctly the diffraction phenomena. 


It is the aim of this paper to clear the problem. 


$ 1. Kottler's Theory 


Let us assume that, at the point L, there is an oscillating electric dipole whose momentum 

> > : 
s parallel to the z-axis. Now, let Fand H denote the magnitude of the electric and magnetic 
elds, respectively. Kottler's theory gives the following formulae for the values of these 


(411) 
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vectors at point P: 


+ er ] e Z2 
E(P) = E*(P) + — grad, — uP(L, P) + x .— uP(L, P) -+ 
la OZ, 4r 
er: EINST pie 
fees (ds x Eo) w (Q, P) + T gradp í (dsH,) w (Q. P). 
r r 
- n tk — 
H(P) — H*(P) — ve [grad; uB(L, P) x x|+ 
ie E VD l 5 FJ > 
des (ds x Hy) w (Q, P) — nen gradp I (ds Eg) w (Q. P) (1.1) 
f P 
wherein 


m = : l EN n : 
FE EN E Py E gradp toll, P)+ or ug(L, P), 


c 
3 
Cx 


E x EE = th > 
H*(L, P) - HL, Py = [grad u (L, P) x 2]. 


within the light cone, and 


E*(L. P) = H*(L. P) = 0 


in the shadow. Here. 


1 * e-ik(e r) (r dt 
uB(L, P) 2 — | —— —— -— - ds (1.2) 
4x. ro or+or 
» 
and 
€ —ikr po 
w(Q, P) = ——— 
TPQ 
X denotes the unit vector jn the direction of the z-axis, t—-one directed along the tangent 


yo I . > tote ` - 1 3 1 1 
to the curve J’, r — the distance between the point of observation P and a point Q on the 
, tg . . . . . 
edge I’ of the diffracting body, o — the distance between the light source L and point Q, and 
> „—ikR| : . 
ug(L, P) =e R, where R denotes the distance between L and P. 
M PR? N 1 ERO . A B 5 c TM 
Fhe terms in eqs. (1.1) containing the function uP(L, P) will be referred toas "Kirchhoff^s 
ene Sud peel ur B ; Nee 
terms", as the function ug(L, P) describes the diffracted wave already existing in the 
Kirchhofl theory. 

Rubinowiez (1917) was the first to transform Kirchhoffs’ surface integral into a curvi- 
linear integral, thus materializing Young's ideas according to which the diffraction image 
arises through interference of the waves "reflected. from the edge of the diaphragm with 
the wave emerging through the aperture without changing its direction of propagation. 


Kottler's formulae can be considered as the formal realization of Young's idea for electro- 


magnetic waves, 
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Method of approximation 


It is easily seen that all the integrals appearing in the formulae of eqs. d l) can be 
transformed into integrals of the following kind: 


fe Ex E ds (2.1) 
Zi 


where =r +ọ, and g(s) is a function of the length of the are of the curve 7, assuming different 
values for different integrals. If we assume the function g(s) to be everywhere regular, then, 
with the aid of the stationary phase method applied by Rubinowiez to the scalar case, we 
can find the approximate value of the integral (2.1). The latter is 


m s $m 2.2, 
Ese) ds eo aes) | / ee 4 (2.2) 
tOr 


there and henceforth the index v, tells us that the values of the corresponding functions 


r 


are taken at the point s=s, in which £’(s)= 
By simple calculations we can verify that 


: #2 l l l x f 

CU m sin? (o, ds) s + — [eos (r, K) + eos (o, K)] (2.3) 

r 0 K Š 

where A denotes the radius of curvature directed along the principal normal to the curve I’. 

Now, in the neighbourhood of the shadow boundary. we cannot apply this method 

g ) PP) 

to the integral appearing in the formula (1.2) because of the singularity of the function 
o(s) in this region (at the shadow boundary ro + ro = 0). 

This difficulty, however, can be avoided (Rubinowicz 1924) and the following formula 
is obtained for the integral (1.2) in the neighbourhood of the shadow boundary within the 
light cone: 

- iun =) V2kn- (0--r—R) 
(de en x uA e 


up Ry2 ez" dv 


Ihe integral (1.2) in the neighbourhood of the shadow boundary (in the shadow) can be 


expressed by the formula: 


=i (uno. 32) Ver cr- RO 
€ 4 


Qu t m RR > 
ne R y2 | 


je 


v? 


dv (2.5) 


"inally, we note that total light wave within the light cone can be written as follows: 


—ikR E (a = —Vkn t (o--r— R) 


m A efr iiu. Wat (2.6) 
u = uB - R Ry? i 


Cc 


$3. Approximated formulae for diffracted dipole wave 
In the present section we shall apply the foregoing approximation method to obtain 
the approximate values of the vectors E and H as given by the formulae (1.1). The entire 


discussion. will be for the x component of the vector E. From eqs. (1.1) we have 


1 3 E 1 * > e-ikr 
7 (P) = ES SUP. + — Š =— Tad EP PAS gne fas E r E^ 


4móOxpóoOzp ` ix. 
T 
1 e f; E e` ikr 
rea ds Ł " 3.1 
Amik xp = (ds Ho) r G.M 
r 


wherein 


É i N) e tke E. k 2 e ike 
1) = — frac T 
0 5 9 920 0 An o 
| | (3.2) 
d. T e tke = | 
= —— | grad; x 
0 - = 0 


In the formulae for Ej, and Hy we carry out all operations indicated and take those terms 
only which are proportional to the highest powers of k. We drop the other terms as being 


very small. Simple caleulations yield 


= k? e-ike Jo a Ks eg uM 
E, - — ——- grad; o + « — 
: Ag 0 oz,” Aa 0 
= k? eike ^ 
H, = —— [grad o x x] (3.2a) 
dsc) oe 7? 


We now once more take up formula (3.1). We already know the approximate value of the 


function u” (L, P) in the region considered (i.e. in the neighbourhood if the shadow bounda- 
. . LI . - r . . I 
ry), since it is given by eqs. (2.4) and (2.5). Thus, we consider the integral 


DOREM ers 5 m e ihr 
(ds x E), — = by Eg. — ds — | t; Es, ds. 
r i F ^ ia 


« 


r pP r 


We can apply formula (2.2) to either of the integrals appearing on the right hand side of 
the above formula; now, however, the role of the function g(s) is played by suitable expressions 


calculated. from formulae (3.2) or (3.2a). On carrying out the calculations we obtain 
I Beh uos (ers 1) Zi S 
dsxE = — IXE, : 
T ( 0)x r r kt" | ols! (3 3) 
y 


where E, is given by (3.2) or (3.2a). 
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We proceed similarly when caleulating the fourth term appearing in the right hand 
4 
[o 
side of formula (3.1). Putting the operator —— under the sign of the integral and applying 


exp 
formula (2.2). we find (on dropping the term l/r with respect to ik) 


On EN, ed ^ 2 iin 1) Ir Qa >> 
ro (ds H,) —— = — ik E” (tH); (3.4) 


Oxp. 7 r xp 


> 


H, in the above formula is given by the (3.2) or (3.2a). 
Taking into account (2.4), (2.5), (3.3) and (3.4) we find the following approximate value 
for the component E, given by (3.1) 


V2kn (y +0 ER) 


re 


"S 
—1—3t n 
t^ 57 3 7f iR D abi 
e © € 1 v 
E,— E; + | De 100 ne 
2 L 4 y 


Ly le de 3 27 n 1 Ej Letra) In Ir LH.) 
| Age 7; Bu alt “Act r kt" xp \ g 


Phe term E; (and similarly H; for the vector H) appears — as we know — only when we 


we within the light cone; we then have to take the sign minus before the expression 
Q og ER 
x, 92L v 
If we are in the shadowed region, this expression should be taken with the sign 


jlus. 
We shall refrain from giving at full length the calculations leading to the other compo- 


ients of the vectors E and H, since these are carried out in a similar way; we will limit 
»rselves only to presenting the results. 
Thus, eqs. (3.1) assume the following approximate form: 


Li (se) ET 


= =r ro) € H —iL y 
MEL Ac c [e = 2 = d SIE 
E= E* +k ES x RR | € J : 
ler —— —i (n 7) Imm 
uU e ( 1) Oe i l € 4 276 
z — — ——————— tH) grad , 
Er P kt” [ ol l Art T ke" ( o) gradp r : 


Ñ = Be a | ennt ue ec UICE i] 


(3.5) 
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é 
The expressions acted on by the differential operators grad, p. 


or grad, are put into 
brackets with the subseript v. The latter denotes that we first are to take the value of the 
mentioned expression at the point P». and then proceed to carry out the differential opera- 
tions indicated. 

As to formulae (3.5), we would like to stress two points. Thus, at first sight one could 
be led to think that in the case of the formulae (1.1) the same reasoning is applicable as 
the one applied by Rubinowicz (1924) in the scalar case. This means that we could inquire 
whether use can be made of the fact that. in the shadow boundary, the following relation 


is fulfilled : 


PODER = RD 
EOL Et = Et 


where E(P and EC? denote the x-component of the electric field in the light cone and in 
the shadow, respectively. On choosing such a standpoint. all further calculations would 
be quite similar to those of Rubinowicz, with the sole difference that we would obtain expres- 
sions more complicated since the formulae for E, contain an operator of double differentiation. 
However, it is impossible to proceed along these lines. as the integrals obtained are divergent 
in the boundary. 

In connection with the formulae (3.5) we have to discuss once more question, namely, 
whether operations such as differentiation and the method of stationary phase commute, 
In order to ascertain that the answer is positive (toa good approximation), consider the 
integral 

ded pan Je Rr FF (r, o)ds (3.6) 


Tr 


# 


In our case, the function /(r, o)=ljro, and we may multiply it by the corresponding 
of the tangential vector (the one tangent to the curve 7"), which, however, we omit here as 
unimportant. The integrals appearing in the formulae (1.1) are exactly of the form (3.6). 

Firstly, we apply to the integral (3.6) the stationary phase method and then differentiate 
the expression. Thus, applying formula (2.2) we find 


^ GR Ber | 
I = e (e RAT, 0) Te 4) (3.7) 
CXL CG v 


It should be remembered that the value of the expression acted on by the differentiation 
5 1 
ə 
operator 5— is taken at point P, in which 
8x, 


di © ER (3.8) 


As it is our intention to differentiate the expression in eqs. (3.7), we have to take into 


account that, by virtue of (3.8), r as well as o (and also Z") depends on x, in another way 
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than subsequent to applying approximation calculation. Since 


e Vx 
"EU 


Yp) + (x — Zp)”, 


yp)? + (s 


xr) + (y Zr)", 


the dependence of r on x, is implicite in the variables x, y, z and the dependence of o on 


Xr — in these variables and in the variable Xr, itself. We must take this into account when 


differentiating the corresponding expression. This is clear also from another point of view. 


In differentiating with respect to xz, we implicitly assume the coordinate x, of the light 
source to be varied. This variation, clearly, affects the position of point P, on the curve J’, 


which is assigned by the dependence of x, y, 


s On xr. 


After these remarks we carry out the differentiation in (3.7), obtaining 


Dre : 
I = fo; r) er e 4 ^ e iko+r)|, + e tke ir) FE 
k£ i» € XL e Xr, 


e 2 


(r 0) 2 MIS | 


Now, let us consider the term 


A) 
9 gik(o--n), 
e xr, 
We have 
e-iKr-Fo)y — — jke iketr)|, = (o d r) " 
CXI, eX 
; NS 2. AO 
= A me de 
a ike-ikG--0), X XL Eu Y jit. C) | 73 E 
o Q Ox 0 xL 
E VOCUM UNE — nk: X — Xr, 
I I 
yh ew p^ Qr. men oa. 0 » 
PE E Fara Eam 
ex BES à 
Grouping the terms in brackets and taking into account that 5 = 5.7, eic. we find 
o Xr SEE XL 
2 [tè ti? ðs ara 
T gcik(etnl, = — ike-ik(4n], al = 
Ox, 0 ? Jus DENT 
ind hence, by (3.8), 
2 X — Xf. 
e ROEN — ike Fety -— —- (3.10) 
Ox, O- jr 


t should be noted that the result of this differentiation is the same as what would have 


een obtained had we assumed o to depend on x, only through the variable x, and not 


hrough x, y, z and r to depend on xz. 
Substituting (3.10) in (3.9), we have 


^ .3T 
Tikes ern Ef 9) 1 / ee 4| + 


0 kt" " 


(3.11) 


-+ eik(e+r)|, 9 
OXL 
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Since k is a very large number, we shall obtain a good approximation for Z if we retain only 


those terms in (3.11) which are proportional to the highest power of k. Finally, we find that 


I = ike-iketn ~—*E gu 9) m 7 (3.12) 


0 ke ly 


We now return to formula (3.6). In order to find an approximate value of the integral J, we 
shall first differentiate this integral and then apply the stationary phase method. We shall 


then have 


2J 


TA. 
` , x— X | NU gr HESS 
I=ik eth (et —— fr, o)ds ip e-ik(etrn) Er yo ds. 
0 4 CO CXL 
r 
Applying formula (2.2) to either integral, we find 
x— X, ACD 
I = ike ETTE fr, 9) 1 Te t + 
0 Ra 3 
7 : rer (ed 
| -ik(o*) I 90 2x -i7| 
E s -i a 
JO CXL kÈ ly 


here, on the basis of the same argumentation, we retain only the first term; this yields the 
following result: 


I = ike-ik(r-+0) ae 0) oe and (3.13) 


The value of the integral / calculated in the two different ways (3.12) and (3.13) is 
thus found to be the same. A similar result could also be obtained for the other differential 


f 2 e re i o ə 2 9 2? 
operators e.g. for > „ and also for the operators Aas ^ ^ š 
l = Oy, > Oz, Ox, Oz, Oy, Oz, oz 


In this last case, only the terms proportional to A? should be retained in the corresponding 
expressions. Finally, we can formulate the following conclusion: differentiation and the 
method of stationary phase do indeed commute, provided we retain in either case only the 
terms proportional to the highest power of k. Hence it follows that, when calculating the 
terms i 

grad, | (...)ds, 
r 
we can first take the integral by means of the stationary phase method and then differentiate, 
or else we can proceed inversely. The results will be identical if the values E, and H, are 
taken from (3.2a). 


In the formulae (3.5), terms of the following type also appear: 


V2kr-r+o—R) 


2 


ô. eikR ag Ed v? 
grad; = —— Pus rs 
(e ZL R A " 
V2kn-! (r--0—R) 
e ikR =. = i ^ 
grad, n e dv x x 


r 


g* 
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In the present case, the differential operators act on the upper limit of integration 


Y, 
V2 (r+o—R) also. Here, of course, we have to take into account the de pendence of r, o 
on xj, yp or zy through the variables x, y, z also. In order to avoid repeating the same (in 
principle) caleulations as those already presented when deriving formula (3.10), we write 
at once that, e.g., 


R 
9 ODORE Nro NER 
e UE R|, rq pce Fak 
Oz, V 7t ( 9 V. Vr d- 0 -— R N 5 à i 


We now return to the formulae (3.5). We can carry out therein all the operations indi- 


cated and thus obtain the concrete form of the vector E (formula (3.5)). First of all, we state 


that, by formulae (3.2a), the two last terms in the formulae for E are of the order of AV A. 


o 
We now consider the term grad, 2 (...) appearing in the formula for E. We shall investi- 
ME Aw 
gate how the operator >— >— acts on the function 
Ox, Oz, 


V2kn— (g--r—R) 
eikr eg 
| R | ems a) . (3.15) 


It is readily seen that possible difficulties can arise when computing the following expressions: 


à Qk up Lip 2k 
Ja y- ee: 8x, Or, V ER En 


since the function V2kz-(r + o — R) vanishes at the limit of shadow (in this case, r + o = R); 
then, the functions are proportional to (r-- o — R)~” and (r+ o — R) ^ respectively, and 
infinities can appear. If the expressions considered had finite values at the boundary of 
shadow, then the largest term obtained on aic A of (3.15) ys be od 


to k2, since then the terms containing the derivatives = 2 +0 — 
AL 


92 


eee Q— — R would be of the order of k, k* and k*, respectively, and could be omitted 
OX, Oz, 
with respect to the term of order 4? appearing on differentiation of the function e FERIR 


This, however, would lead to the conclusion that the largest contribution to the value of 


the x-component of the vector £ in the neighbourhood of the boundary of shadow is due 
to the part of the "Kirchhoff term" which is proportional to 4?. The above conclusion 


would be correct for the remaining components 7 the vectors Æ and H also. 
=) 
Almen: KE 


OZL "Ox, OZ. 


AR have 


We now have to prove that the expressions 


finite values at the boundary. of shadow. 
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By (3.14a) we have 


^L Er e d. ep —^4L 
e / 60 IR 
— Vro =R= ——— | 
92] 2Vy$ —R » 
Le. 
d [mr l p s oe R | cos (R, x) — cos (o. x) 
an Re : 
Oz EC 4 ro d cos (3 ro) 


where cos (o. x) etc. denotes the cosine of the angle between the vectors 0 and z (%15 the 
vector directed along the z-axis). 

Now, applying the well known theorem for the trihedral angle formed by the vectors 
" R, x and expressing cos (0, R) and sin (o. R) through trigonometric functions of the 


angle subtended by the vectors o and r, we finally obtain 


9 , * 9 TO M ro 
3 = ARS 1 r+o+R Ar? cos (Rx) sın“ = cos D i 
Ae 4 ro R(R+Y R? —?? sin? (or) 
9. 3 
— = sin (z sin (Rz) cos (oR: Rz) (3.14b) 


In this formula cos (o R; Rx) denotes the cosine of the angle formed by the plane containing 
the vectors 0, R and the one containing the vectors R, x. On approaching the boundary 
of shadow this angle tends to the one formed by the plane tangent to the surface of the cone 


of reflection and the plane containing the vectors R and x. On this cone, the axis of which 


* ` . . . LI 
is directed along the vector ds and the apex of which is the point P, are situated all the 


points of observation having the following property: cos (r, ds) — — cos (9, ds) (see $2). 
The words "on approaching the boundary of shadow" are always to mean that the point 
r 0 2 x 
P glides over the cone, and — — — — ; the point P, remains constant. 
í Ln 0 


~~ € . c 
On the basis of formula (3.14b) we conclude that expression —— yr -- o — R possesses 
ez, 
a finite value at the boundary. of shadow. 


2 2 en 
5 2 : : e e S eg 
In order to prove that, in this region, the expression —— Vr+ o — R also has finite 
OXL SSL 


value, we diflerentiate the right hand side of formula (3.14b) with respect to xz. When 
diflerentiating, we naturally have to take into account the comments preceding 
formula (3.14a). 

It can be easily proved that all expressions obtained are finite with the exception of 
the derivative 4 


cos (oR; Rx); in this case nothing can be said a priori. We now have 
XL 


to prove that this derivative too, has a finite value at the boundary of shadow. This means 
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we must prove the following expression to be finite in this region: 


8 | öx? | 
9x, LIOXT| |Rxx| 


"Eck 2) p Vv = 
, TE. IE. cr 
(the vectors 7, o, R are coplanar). It is sufficient to prove the vector —— < = to be finite, 


M Oxr, |o X T| 

: 9 Rxx 
since the vector u— Eu 
o Xr, |R x x 


> 


4 
Qr 
Q X 


has this property certainly. Here we take note, moreover, 


that the vector n — 5 tends towards the vector perpendicular to the plane tangent 
r 


to the reflection cone as we reach the boundary of shadow. 


MEER ER, 
In order to prove that the vector —— Ede an = is finite, we define the following three 
xr |o X T| 
: : M uu ac f 
mutually perpendicular axes directed along the vectors —, 7, — x n. As can be easily seen, 
r r 
the following three formulae are correct 
4 On 
n = (0, Shae ge 
XL 
T 
> > e = > 
r on r > Tee 
since — n = 0 
r Ox Ox r 
(ie AA C E Ur T ET I ate + 090 
—xn npe nm ; 
r Oxr, [o xr Leer Ox Ixy 


since n | 


— + 


dw 
X 
St 
nl 
ll 
S 


5 n ó j zd Uf T 
The projections of the vector — on the directions n and — assume finite values at the 
OXL if 
eg ðn : 
boundary of shadow. Thus we have to prove that projection of the vector — on the di- 
XL 
r > a 6 : : 2 > 
rection — x n has also a finite value, i.e. that the right hand side of the third of the above 
r 


formulae is finite. 
or +ds dg > Qs > 


Making use of the fact that =t ; =t — i (where 7 is the unit vector 
IX Wb Ox T IX L ex L 
~> os 
in the direction of the x-axis), and adding and subtracting the term o(nt) 7 we can write 
235 
Po aN aa nt 9s ro 
TAn = ctg Ms de 
r Ixy Toe ur 2 
>> Ex Is at Is 
r(nt) — r(nt — o(nt 
(n 1) (nt) Dar o(n t) Oxy (3.16) 


| => > 


lg x r| 
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The first term on the right hand side of eqs. (3.16) vanishes at the boundary of shadow; 


thus we have to investigate the behaviour of the second term in this region. 


> > 
mM: : om ot 
Taking into account that a aS 
: 7 0 
oF ti (ot) (0-1) 
es ex 0 0 o? 
== 2 " = str 
Oxy ol s 
AN 
es 


(see eqs. (2.3)). the second term of the right hand side of eqs. (3.16) can be expressed thus 


^ ^ 
H „= Qs ~ = os 
r(n1) — r(nt) Z— — o(n- t) - 
i CXL oxy 
DAT 


r(n 1) [cos (o. K) + cos (r, K)] E 


DRE 


Bee gare o- t) (0 A, 
[(— nt) | (ti) — ey) (QUI , (ni) sin? (o. ds)] 
amo (r+ 0) Riz 
j er> n c ` x PS (3.17) 
GC lox TI 0 


The first term of the right hand side of eqs. (3.17) is finite at the boundary of shadow since, 


as can be easily verified by applying the theorem of cosines to triherdal angle, cos (o. A) + 


Y TD , $ : x : : 
- cos (r, K) cos =. Finally. the following expression remains to be investigated: 
= g ex] 
++ (0 t) oi) dem 
(ti) e tet) + (ni) sin? (o. ds) 
ee? 0 0” i 
(—nt)- — = 
| X r| 
ay € 
Taking into account that sin?(o, ds) = 1 — [ 7 | and using the vector relationship (a x b)(e x 
Q 


E 


x d) = (ac)(bd) — (ad)(be) we can rewrite this expression as follows: 


+ t)]} 


» A à E TE — > > ><, z P 
Now, applying twice the vector formula a x (b x c) = b(ac) — c(ab) we find the following 
value of the expression under consideration: 
> -= "ma "P" ~ 
: (i x t) o (ot) 


f ar 
20°r sin? = 


in transforming this equation use was made of the fact that "^ + 9^ ol. 


r 0 
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Thus, we have proved that at the limit of shadow and in its neighbourhood the expression 
E 2 


dx, 8 


of DERA of the terms in Kottler’s formulae have thus been proved, since we can likewise 


eh +o — — R has a finite value. All the foregoing conclusions concerning the orders 


verify that the situation. is the same in the case of the remaining components of the vectors 
E and H. 

These considerations lead to the following practical conclusion: in order to find the 
intensities of the electric and magnetic fields in the neighbourhood of the boundary of shadow 


on the basis of eqs. (3.5) it is sufficient to take into account only the expressions containing 


6) eikR 
© 
the operators grad, —— 2: and grad, and, in these expressions, to differentiate the function R 
"E 
only. Thus, the intensity of either field in the region under consideration expressed by the 
eikR 


corresponding derivatives of the function - multiplied by Fresnel’s integral. This fact 


R 
is in full agreement with experimental data on the basis of which the intensity in the diffrac- 
tion image agrees with the theoretical data as calculated on from the scalar theory by 
Kirchhoff. 

Thus, thanks to the application of the stationary phase method, we have been able 
to derive approximate formulae for the vectors E and H. These formulae are much simpler 
than those of eqs. (1.1). We note moreover that, although eqs. (1.1) and the approximations 
obtained thereform are derived on the assumption of a dipole vibrating in the direction 
of the z-axis, they can be generalized for the case of a dipole vibrating in any direction. 
Obviously, the corresponding formulae for a magnetic dipole with arbitrarily directed axis 
can be given immediately. In the case of a number of points P, situated on the edge Z’, the 


electromagnetie field at P is given by 


E=) E. 
H=>H,. 


The author is very greaily indebted to Professor A. Rubinowicz for suggesting the 
subject of the present investigation and for his numerous critical and friendly remarks, and 


for his extremely helpful discussions. 


REFERENCES 


Kottler, F., Ann. Phys. [Leipzig], 71, 457 (1923). 

Rubinowicz, A., Ann. Phys. [Leipzig], 53, 257 (1917). 

Rubinowicz, A., Ann. Phys. [Leipzig], 73. 339 (1924). 

Rubinowicz, A., Die Beugungswelle in der Kirchhoffschen Theorie der Beugung, PWN, Warszawa 1957. 
Sommerfeld, A., Vorlesungen über theoretische Physik, Bd. IV, Optik, Wiesbaden 1950. 

Stratton, J. A., Electromagnetic Theory, N. Y. 1941. 


Vol. XX (1961) ACTA PHYSICA POLONICA 


Fasc. 5—6 


ZUR FRAGE DER EXOELEKTRONENEMISSION VON GEPULVERTEN 
KCL-EINKRISTALLEN, DIE VORERST UNTERSCHIEDLICH 


VERFÄRBT WURDEN 


von B. SUJAK 
Institut f. Exp. Physik, Universitat Wrocław 
UND A. SZAYNOK 


Physikalisches Institut der Technischen Hochschule, Wroclaw. 
(Eingegangen am 30 November 1960) 


Die photo- sowie thermostimulierte Coelektronenemission von gepulverten KCL-Einkris- 
tallen, die vorerst unterschiedlich verfärbt wurden, wurde untersucht. Die Verfärbung der 
Kristalle wurde additiv mit Potassiumdampf oder längere Röntgen sowie ß-Strahlen Bestrahlung 
erzeugt. 

Die photo- und thermostimulierte Coelektronenemission von Pulvern aus additiv ver- 
färbten Kristallen zeigte ein unterschiedliches Verhalten gegenüber der Emission, die von den, 
aus mit X- oder /-Strahlen vorerst verfärbten Kristallen vorbereiteten Pulvern gemessen wurde. 

Bei den mit X- und besonders mit fStrahlen behandelten Kristallen ist die gemessene 
Anfangsintensität der photostimulierten Coelektronenemission sehr stark Zeitabhängig. Wenn 
der Zeitabstand zwischen Verfürbung und Pulverisierung relativ kurz und die Zeit der P-Bestrah- 
lung lang ist, so werden sehr kleine Anfangsintensitäten der photostimulierten Coelektronenemis- 
sion gemessen. Die messbare Emission ist in diesem Falle aber sehr kurzfristig. Die photostimu- 
lierte Coelektronenemission von Pulvern aus additiv verbärbten Kristallen ist dagegen relativ 
kaum Zeitabhängig. 

Die thermostimulierte Coelektronenemission von Pulvern aus mit X- und f-Strahlen zeigt 
gegenüber den aus additiv verfärbten Kristallen nur schwache Intensität bis zu den Temperaturen 


von etwa 350°C. 


Einführung 


Die Untersuchungen der Verteilung von elektrischen Aufladungen der Pulver-Kerne 
einer Pulver-Wolke, die aus vorerst verschieden verfärbten KCl-Einkristallen vorbereitet 
wurde, lieferten Ergebnisse (Szaynok 1954, 1958, 1960), die mit einer anderen Methode 
nachzuprüfen waren. 

In den Pulver-Wolken, die aus Pulvern im Schwebezustand entstehen, konnten nämlich 


kaum Pulver-Kerne gefunden werden, die einen Photoeffekt für das weisse Licht zeigen, 


(425) 


426 


wenn das Ausgangsmaterial (KCl-Einkristalle) vorerst mit f-Strahlen verfärbt wurde 
(Szavnok 1960). Wenn dagegen das Ausgangsmaterial zur Vorbereitung der Pulver vorerst 
mit X-Strahlen oder additiv im Potassiumdampf verfärbt wurde, so wurden stets Pulver- 
Kerne gefunden, die den Photoeflekt zeigten. 

Das Erscheinen des photoelektrischen Effektes aus verfärbten Alkalihalogeniden 
ist seit langem bekannt (Elster und Geitel 1896, Lukirsky, Gudris und Kulikowa 1926, 
Fleischmann 1933, Asmus 1936) und wurde erneuert studiert (Sujak 1953. 1955, 1956, 1957, 
Bohun 1960) in Zusammenhang mit der Erscheinung. die unter der Bezeichnung „Exoelek- 
tronenemission“ bekannt geworden ist (Exoelektronentagung 1957). Ausserdem wurden 
einige Studien des photoelektrischen Effektes von Alkalihalogeniden durchgeführt (Apker 
und Taft 1950. 1951, Arsenjewa-Geil 1957). 

Aus diesem Grunde schien es sehr interessant zu sein im Falle der Verfärbung der 
KCl-Kristalle mit -Strahlen das „Fehlen“ des photoelektrischen Effektes mit Hilfe der 
Exoelektronenmethode zu untersuchen. In dem vorgelegten Bericht beschäftigen wir uns 
also mit diesem Problem. Nach einer kurzen Beschreibung des Verhaltens der Pulver- 
Kerne im Schwebekondensator nach Kunkel und Hansen (1954) geben wir die Ergebnisse 
wieder, die mit Hilfe der Exoelektronenuntersuchungsmethoden für das gleiche Versuchsma- 


terial erhalten wurden. 


Verhalten der Pulver-Kerne im Schwebekondensator 


Die KCl oder NaCl-Einkristalle können nach der Pulverisierung eine Staubwolke 
im freien Fall geben. In dieser Staubwolke vorkommen nun Pulver-Kerne, die negativ, 
positiv oder elektrisch neutral sind. Die Grössen der Ladungen der einzelnen Pulver-Kerne 
können mit Hilfe einer Anordnung nach Kunkel und Hansen (1950) gemessen werden. 
In dieser Anordnung steht der Schwebekondensator senkrecht und die abfallenden Pulver- 
Kerne werden mit Wechsellicht beleuchtet und photographiert. Somit kann die Bahn der 
einzelnen Pulver-Kerne festgelegt und die Ladung errechnet werden, wenn die Spannung 
und die Geometrie des Kondensators, sowie die Grösse der Pulver-Kerne gegeben werden. 

Bei der Staub-Wolke, die aus nichtverfärbten Einkristallen von KCl und NaCl vorbereitet 
wurde, ist symmetrische Verteilung der Ladungen, die durch die Staub-Kerne getragen 
werden, festgestellt worden, Dieses bedeutet, dass in der untersuchten Staubwolke die gleiche 
Menge von negativer wie positiver Ladung vorkommt. 

Wenn dagegen die Ausgangskristalle vorerst verfärbt wurden (F-Zentren oder ähnliche 
Kristalldefekte eingeführt) so zeigt sich in der Ladung der Staubwolke eine klare Unsymme- 
trie, die von der Intensität und der Methode der Verfärbung abhängt (Szaynok 1958, 1960). 

Während der Messungen der Ladungsverteilungen an Staubkernen, tritt eine photo- 
elektrische Wirkung des Belichtungslichtes ein. Dieser Photoeffekt tritt besonders bei den 
Staubkernen klar hervor, die aus den additiv (in K-Dämpfen) verfärbten KCl Einkristallen 
gepulvert wurden. Ähnlich verhalten sich auch Staubkerne, die aus mit X-Strahlen ver- 
färbten KCI Einkristallen gepulvert wurden. In der benutzten Anordnung nach Hansen und 
Kunkel zeigt sich diese Photoemission in einer Krümmungsveränderung der Bahn des 
gegebenen Staubkernes innerhalb des homogenen Feldes des Kondensators hervor (Abb. 1.) 
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Mit der Zeit der photographischen Beobachtungen fällt die photoelektrische Einwirkung 

des Lichtes ab. In der Tabelle I wurde eine Messreihe für KCl-Staub in dieser Hinsicht 
zusammengestellt. 

In den Staubwolken, die aus vorerst mit ß-Strahlen verfärbten KCl-Einkristallen 

vorbereitet wurden, konnte Szaynok kaum Staubkerne finden, die einen Photoeffekt zeigten. 

Dieses hat die Vermutung hervorgerufen, dass das photoelektrische Verhalten der einzelnen 


Staubkernen (bei gleichen Versuchsbedingungen) von der Methode der Verfärbune abhänet. 


\bb. 1. Das photographische Bild der Abfall-Bahn des einzelnen Pulver-Kernes bei welchem die photoelektrische 
Emission stattfindet 
TABELLE 1 
KCI (verfärbt mit K-Dampf) F-Zentrenkonzentration etwa 1018 em ë 


Die Zahl der Staubkernen, Die Zahl der Staubkernen, | 


"nz a | | N 
Es | die einen Photoeffekt zeigen | die keinen Photoeffekt zeigen = m 
Nr ; | INES Esa Aj 
| Noh No | ph ? 
ESS | 15 5 0,75 
6—10 | 15 | | 0.94 
1115 24 10 0,71 
16—20 14 10 0,58 
gp | 15 | 12 | 0.56 
26—30 | 12 12 0,50 


Der vorgelegte Bericht enthält Ergebnisse der Untersuchungen der photoelektrischer 
nd thermischer Emission aus dem gleichen Versuchsmaterial, die mit der Exoelektronen- 


missionsmethode durchgeführt wurden. 


Die Exoelektronenemission 


Die photoelektrische Emission von kleiner Intensität, die erst nach einer vorherigen 
rregung statt findet und mit der Zeit der Messung abfällt, ist mit dem Namen der photosti- 
wulierten Co-Elektronenemission (oder Exoelektronenemission) bezeichnet (Exoelektro- 


entagung 1957). Wenn solche, erregte Emission, im Dunkeln während einer gleichmässigen 
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Erwärmung des Emitters gemessen wird, so wird sie thermostimulierte Exoelektronene- 
mission genannt. Im allgemeinen werden dazu Detektoren benutzt, die sehr kleine Intensità- 
ten der Emission von etwa 10-19? A zu messen erlauben. Deswegen sind die Messungen der 
“xoelektronenemission sehr oft mit Hilfe des offenen Luftspitzenzählers geführt. 

l. Die photostimulierte Exoelektronenemission aus den KCl-Pulver- 
kernen. 

Auch in unserem Falle wurden die Messungen der photostimulierten Exoelektro- 
nenemission (die weiterhin als Ph. Ex. Em. bezeichnet wird) der gepulverten Kristallen 
mit Hilfe eines offenen Luftspitzenzählers geführt. Die Proben kammen etwa 3 mm vor 
der Öffnung des Zählers zu liegen. Zwischen der Probenunterlage und dem Zylinder des 
Zählers wurde eine Spannung von etwa 70 V angelegt um die Exoelektronen zu beschleunigen. 
Die Impulsgeschwindigkeit N/t wurde wie üblich auf einem  Papierband automatisch 
registriert. Die Ph. Ex. Em. wurde unter der Belichtung mit weissem Licht einer 6 V, 
35 W Wolframlampe gemessen. 

Vorbereitung der Pulver. Die Pulver wurden aus synthetischen Einkristallen 
von KCl, die nach der Kvropoulos Methode gezogen wurden, vorbereitet. Bevor die Pul- 
verisierung mit einem Porzellan-Stock im Porzellan-Mörser erfolgte, wurden die Kristall- 
stücke entsprechend verfürbt. Die verfärbten Kristallstücke wurden immer Lichtdicht 
aufbewahrt. Auch die Pulverisierung erfolgte immer in Dunkelheit. 

Die Einkristallspaltstücke von KCl wurden auf die drei Arten verfärbt: durch Diffusion 
der A-Dümpfe im Vakuum, durch X-Strahlen Bestrahlung (40 kV. 8 mA, Cu-Anode) 
oder durch Bestrahlung mit einem Strontium-Präparat (3Sr, 0,54 MeV, 10 mC). 


Die Spektralanalyse der benutzten KCI-Einkristallen ergab die folgende Beimengungen: 


Na — etwa 0,01%, 


Fe Spuren, 
Ca weniger als 0.001%, 
Sr Spuren, 


M nicht feststellbar. 


Messergebnisse. Die Einkristallstücke von KCl, die im &-Dampf additiv verfürbt 
wurden, zeigten eine F-Zentrenkonzentration von etwa 1018 em-?, Dieses wurde mit Hilfe der 
Lichtabsorptionsmethode ermittelt. Die Verpulverung und gleich darauf folgende Messungen 
der Ph. Ex. Em. erfolgten mehrere Tage nach der Verfärbung. Der zeitliche Verlauf der 
Ph. Ex. Em. zeigt die Kurve 4 auf der Abb. 2. Die Intensität der Emission fällt mit der 
Zeit relativ sehr langsam ab, wie das zu sehen ist. 

Die Einkristallstücke von KCl, die mit N-Strahlen (40 kV, 8 mA, Cu-Anode) verfärbt 
wurden, lagen 1,5 Stunden der Strahlung (etwa 15 em von der Róntgenróhre entfernt) ausge- 
setzt. Die erziehlte F-Zentrenkonzentration ist leider nicht bekannt. Die Verpulverung 
der in Dunkelheit aufbewahrten Kristallstücke erfolgte etwa 3 Stunden nach der Beendigung 
der X-Strahlen Bestrahlung. Der Verlauf der Ph. Ex. Em. zeigt eine sehr schnelle Abnahme 
der Intensität mit der Zeit, wie dieses die Kurve B auf der Abb. 2 wiedergibt. 

Die Einkristallstücke von KCl, die mit P-Strahlen des Strontium-Präparates (er. 


0.54 MeV, 10 mC) 107 Stunden lang bestrahlt wurden. zeigten eine F-Zentrenkonzentration 
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von etwa 2,2 - 1016 em, Die Verpulverung und Messungen der Ph. Ex. Em. erfolgten aus 
technischen Gründen erst in etwa 20 Stunden nach der Beendigung der fi-Bestrahlung. 
Den zeitlichen Verlauf der gemessenen Ph. Ex. Em. zeigt die Kurve C auf der Abb. 2. ns 
Intensität ist sehr klein und ist schon nach etwa 30 sec nicht mehr messbar. Ein anderes 
Kristallstück von KCl wurde mit dem gleichen Strontium-Práüparat etwa 60 Stunden lang 
bestrahlt. Die dadurch erzwungene F-Zentrenkonzentration betrug etwa 1,8 - 1019 cm 73. 


A 
Nr | 


KCL* K 


intergrund 
jen T = —r T T T m 
3 4 5j 6 t [min] 


\bb. 2. Die Ph. Ex. Em. (in willkürlichen Einheiten) von Pulvern, die aus KCl-Kristallen vorbereitet wurden, 
die vorerst auf folgende Weise verfürbt wurden: 

curve A — additiv im Potassiumdampf, Kurve B — durch X-Strahlen Bestrahlung, Kurven C und D durch 

i-Strahlen Bestrahlung. Weitere Einzelheiten sehe Text. Die Intensitäten der einzelnen Kurven sind nicht 


untereinander zu vergleichen. 


Nachdem der Kristall etwa 70 Stunden in Dunkelheit aufbewahrt wurde, wurde er verpulvert 
nd die Ph. Ex. Em. gemessen. Die Kurve D (Abb. 2) gibt die gemessene Intensität der 
yh. Ex. Em. wieder. Wie es zu sehen ist, ist der Verlauf der Kurve D, sehr dem Verlauf 
er Kurve C ühnlich. Die Kurve D ist nur in der Richtung der höheren Werte der Intensität 
er Ph. Ex. Em. verschoben. 

Die Photoempfindlichkeit der Pulverkerne steigt mit der Zeit der Lagerung des Pulvers 
owie der Kristalle an. Diese Erscheinung ist wahrscheinlich mit dem Abfliessen der nega- 
ven Raumladung, die als Folge der /-Bestrahlung im Innern des Kristalles zum Erscheinen 
ommt. : 

2. Die thermostimulierte Exoelektronenenemission aus den KCI- 
ulverkernen. 

Es wurden auch Untersuchungen durchgeführt um festzustellen ob die Pulverkerne 
us mit verschiedenen Methoden verfärbten KCl-Einkristallen auch ein verschiedenes 
erhalten der thermostimulierten Co-Exoelektronenemission aufweisen. Diese thermosti- 
wulierte Co-Exoelektronenemission wird weiterhin als Th. Ex. Em. bezeichnet. 

Zu den Messungen der Th. Ex. Em. wurde ein begitterte Luftspitzenzähler benutzt. 
ie Temperatur des Zählermantels war konstant auf etwa 50°C gehalten. Zwischen dem 
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N/t 


4 
N/t ! F 
N/t °C 
KCL + K 
d + 440 
5 Einkristall 380 
+ 320 
7 260 
200 
1 
2 a) E 30 
0 T 20 
0 ep 3 2 1 0 
4 t [min] 
N A 
A T vi E 
N/t ec C 
KCL* XY 
4 440 
sus KCl + fs 440 2 A Eınkristall 
eS * Pulver 380 
N 
320 
260 
200 
c) 140 
80 
0 20 
3 EA 1 0 me 1 0 
t/min] t [min] 
Abb. 3. Abb. 4. 


Abb. 3. Die Th. Ex. Em. (in willkürlichen Einheiten) von Pulvern, die aus KCI-Kristallen vorbereitet wurden, . 
die vorerst auf folgende Weise verfürbt wurden: . 

a) — additiv im Potassiumdampf, b) — durch X-Strahlen Bestrahlung, c) — durch P-Strahlen Bestrahlung. 

Für weitere Einzelheiten sehe Text. Die Intensitäten der einzelnen Kurven sind nicht untereinander zu vergleichen. 


Abb. 4. Die Th. Ex. Em. von KCI-E 


Ankristallen, die vorerst auf folgende Weise verfürbt wurden: 1 
a) — additiv im Potassiumdampf, b) — durch X-Strahlen Bestrahlung. Die Intensitäten der einzelnen Kurven 
sind nicht untereinander zu vergleichen. Für weitere Einzelheiten sehe Text. 
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Gitter und der Heizunterlage der Proben (Entfernung etwa 3 mm) wurde eine Exoelektronen 
beschleunigende Spannung von etwa 70 V angelegt. Die Impulsgeschwindigkeit N/t wurde 
wiederum als Funktion der Zeit automatisch registriert. Die Temperatur der Heizunterlage 
wurde von etwa 20°C bis zur 450°C beinahe linear mit der Zeit. mit einer Geschwindigkeit 
von 2°C/sec erhöht. 

Das gegebene KCl-Kristallstück wurde in drei Teile gespaltet und diese mit je eine 
der 3 Methoden verfärbt. Danach wurden von den verschieden verfärbten Kristallstücken 
Teile abgespaltet und pulverisiert. Nach der Pulverisierung wurden die Th. Ex. Em. Kurven 
in voller Dunkelheit aufgenommen. Es stellte sich heraus, wie weiter gezeigt wird, dass die 
Th. Ex. Em. Kurven unterscheiden sich untereinander merklich. Die typische Kurven 
der Th. Ex. Em. sind auf den Abb. 3 wiedergegeben und zwar für: a) mit K-Dampf b) mit 
X -Strahlen und c) mit -Strahlen verfärbten KCl-Kristallen die verpulvert wurden. 

Das im A-Dampf verfärbte KCI-Kristallstück zeigte eine F-Zentrenkonzentration 
von etwa 1.01 - 1018 cm ?. Die Verpulverung und Messung der Th. Ex. Em. erfolgte mehrere 
Tage nach der Messung der F-Zentren Konzentration. 

Das mit f-Strahlen verfärbte KCI-Kristallstück zeigte eine F-Zentrenkonzentration 
von etwa 2.15 - 1016 em”. Diese Verfärbung wurde nach einer 104-Stundigen Bestrahlung 
mit /-Strahlen des Strontium-Präparates erreicht. Die Verpulverung erfolgte in etwa 20 
Stunden (aus technischen Gründen) nach der Beendigung der Bestrahlung. 

Die Verpulverung des KCI-Kristallstückes, das 2 Stunden lang mit X-Strahlen Belichtet 
wurde (40 kV, 8mA, Cu-Anode), erfolgte nach etwa 5 Stunden. 

Um festzustellen ob die Unterschiede in den Th. Ex. Em.- Kurven mehr den unter- 
schiedlichen Verfärbungsverfahren als dem Pulverisierungsprozess zuzuschreiben sind, 
wurden auch Messreihen an kleinen Spaltstücken der additiv und photochemisch verfärbten 
Kristalle durchgeführt. Das Volumen des jeweiligen Probestückes betrug etwa 2 mm? 
bei einer Dicke von etwa 0,5 mm. Es wurden jedoch die gleichen Verhältnisse gefunden 


(Abb. 4) wie im Falle der Pulver. 


Abschliessende Bemerkungen 


Aus den Ergebnissen ist die Folgerung zu ziehen, dass der verschiedenen Ladungsvertei- 
lung in der Pulverwolke aus auf verschiedene Weise verfärbten KCl-Kristallen, auch ge- 
wissermassend verschiedenes Verhalten der Exoelektronenemission beisteht. 

Die Unterschiede in Verläufen der Ph. Ex. Em.-und der Th. Ex. Em.-Kurven sind 
den Messergebnissen nach mehr der Vorgeschichte der Einkristalle als dem Pulverisierungs- 
process selbst zuzuschreiben. Dieses ist in guter Übereinstimmung mit den letzt publizierten 
Ergebnisse nvon Bohun und seinen Mitarbeiter (Bohun 1960, Bohun und anderen 1960) wo 
diese Autoren zu den Schlussfolgerungen gelangten, dass die thermische Vorgeschichte 
der Probekristallen die Ph. Ex. Em. und Th. Ex. Em. beeinflust sowohl wie die Dauer der 
Bestrahlung mit X-Strahlen. 

Es bestätigte sich auch Teilweise, dass in manchen Fällen die Photoemission aus den 
Pulverkernen, die aus mit ß-Strahlen verfärbten Kristallen hergestellten Pulvern stammen, 
sehr klein sein oder sogar ausbleiben kann (z.B. kurz nach der erfolgten -Bestrahlung). 
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Dieses, sehr Zeit abhängiges Phenomenon ist auf die erscheinende Raumladung zurück-- 
zuführen. Es ist wahrscheinlich, dass die negative Ladung, die die $-Strahlen in den Pulver-. 
kernen erzeugen, verhindert die Photoelektronen den Kristall zu verlassent. 

Die Messergebnisse zeigen auch, dass die für eine Verfärbungsmethode erhaltenen 
Exoelektronenemission Ergebnisse dürfen nicht einfach auf andere Verfärbungsmethoden 
extrapoliert werden, wie es in den Anfangsstadien der Untersuchungen gedacht wurde: 


(Exoelektronentagung 1957). j 
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The present paper deals with the theory of the second virial coefficient B(T) = EX x; Aj BG) 
of the equation of state for polar gas mixtures. General formulas for the virial coefficients B 
are derived, containing, in addition to the contribution arising from the presence of intermolecular 
central forces, others due to the various tensorial forces acting between polar molecules: dipole- 
dipole, dipole-quadrupole, dipole-octupole and quadrupole-quadrupole interaction; moreover, 
dipole-induced dipole and quadrupole-induced dipole inductive interaction. The formulas are 
given in two forms: the general form as derived by means of tensor formalism and relating to 
polar molecules of arbitrary symmetry and arbitrary central forces potential, and one obtained 
by specialization for molecules presenting the axial symmetry and a Lennard-Jones (6—12) 
potential. In particular, the latter formulas are applied to one-, two- and three-component gas 
mixtures. 


1. Introduction 


Kamerlingh Onnes proposed an equation of state for imperfect gases, of the form of 
the virial expansion 


ee Be) 


RT - y |y y3 oe 


wherein p, V, Rand T are the pressure, molar volume, gas constant and Kelvin tempera- 
ture, respectively. B(T), C(T), D(T), etc. are the second, third, fourth, ... virial coefficient, 
If the values of the virial coefficients B, C and D accounting for the experimental divergences 
of imperfect gases from the equation of state of a perfect gas are known, statistical mechanics 
will provide valuable information on the intermolecular forces. 

The classical theory and quantitative discussions of the second virial coefficient B(T) for 
intermolecular forces of various types are given in papers by Lennard-Jones (1924), Stock- 
mayer (1941), Rowlinson (1949), Pople (1954), Buckingham and Pople (1955), Stogryn 
ind Hirschfelder (1959), and others (see, Hirschfelder et al. 1954). 

The third and fourth virial coefficients, C(T) and D(T), are discussed by Rushbrooke 
und Scoins (1951), Kihara (1953), and the fourth — by Katsura (1959). 
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The present paper brings the theory and a discussion of the second virial coefficient 
for a mixture of polar gases. In computing B(T), account is taken of the intermolecular 
central forces and of various types of directional forces, the latter being considered to re- 
present a perturbation in the central forces. In order that the formulas be correct for mole- 
cules of arbitrary symmetry, all contributions to B(T) from the directional intermolecular 
forces are computed by means of the tensor formalism. The general formulas thus derived 
for B(T), which hold for arbitrary central forces, are discussed for the cases of dipole and 
quadrupole molecules having the axial symmetry, on the simplifying assumption of central 
forces given by a Lennard-Jones (6—12) potential. The latter formulas are applied to certain 
special cases, namely, a one-component gas, a two-component gas mixture and, finally, 


a three-component gas. 


2. The general expression for the second virial coefficient of a mixture of gases 


In the case of a gas mixture, the second virial coefficient B(T) in the equation of state (1.1) 
is given as follows: 


B(T) = >) ar”, (2.1) 
ill 


wherein » is the number of components of the mixture, and x; is the mole fraction of the 


3 
i-th component in the gas mixture, b x; = 1. 
i=1 
Classical statistical mechanics yields B(T) in the form of 


(12) 

N A- 
VS (12) i 
B(1 ) -—— 302 N e kT — a dı a do!dof? (2.29 
with N denoting Avogadro's number, k — Boltzmann's constant, T — the Kelvin tempera- 

1 . ms . 
ture, u£? — the total potential energy of interaction between the molecules "1" and “2” 
of species “i” and species “j”, r$» — the vector connecting their centres, and o, wo” — 
variables describing their orientation; Q = fdo; = Jdo; i 
The total potential energy can be expressed in the form 
(12) __ ,,(0) 7 (12) 

u; = u (rg) - Vi "(ra, O; 0) (2.3) 


0 . : 
uU (rj) being the potential energy due to the central forces (angle-independent potentials), 
: 12 ; ! š " 5 
and lf? — the energy due to the tensorial forces (angle-dependent potentials), Considering 
V2). represe mn ats x : 
Vy to represent a perturbation in the energy u(r), the second virial coefficient of eq. (2.2) 
can be put in the form of an expansion: 


oo 
BO — pi. M pan. (2.4) 
n=1 
wherein the zero term 
d» 7, _ 
By = — 2aN f le" IET — 1j "dr; : (2.5) 
0 
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is the second virial coefficient for a central force field u (r..), and the subsequent terms 
given by 


(UD N uU rij) 
By E 2025! (o (rij, @;, 0)" e - E S dri? do! dw? (2.6) 


are contributions thereto resulting from the tensorial (directional) forces. 
For the sake of simplicity, u® (rj) will be assumed throughout the present investiga- 
tion to be given by the following Lennard-Jones potential (1924) (see, Hirschfelder et al. 1954): 


le) - (82) an 


wherein the parameters o; and ¢;,, possessing the dimensions of a length and energy, respec- 
tively, are constants characteristic of the chemical species of the interacting molecules. The 
force constants between unlike molecules relate to those between like molecules by the 


empirical combining rules 
o 


g; = t (o; 9j); Ej = (e; £j) ^ (2.8) 


Applying the Lennard-Jones method (1924), and by eq. (2.7), B, is obtained in the form 
(see, Buckingham and Pople, 1955) 


B® = 2 3 No. 72267» (2.9) 
wherein 
1 1 | a ley 
Ayes AY nee AN Py ey Peu) a 
F(yy) m {Phat 2 Heri) s CT LT (2.10) 


s a function tabulated in the monograph by Hirschfelder et al. (1954), whereas 


oo 
oo 
: E [E n = 6m = n— i 
=g Er Ee ELSE 2 3 Yi 
H,(yi) = 12y$0ij J tyne lom nj! ) rgdrij = Vij 2x m! 


0 
ale 


ire functions introduced by Pople (1954) and tabulated by Buckingham and Pople (1955) 
or all n and for values of y from 0.6 to 3.2. 

In the subsequent sections, the terms BG? as given in general by eq. (2.6) will be comput- 
.d for various directional forces acting between two dipole or two quadrupole molecules. 


3. Dipole anisotropically polarizable molecules 


The potential energy of the electrostatic interaction between two unlike, dipolar, anisotro- 
ically polarizable molecules is given by the expression (cf., Barker, 1953) 


; Doa TDT sg) (0051) (253) 4,051) ,, aM PG) PO) 4. 
Vi» = e ie z TE (us o ut PF Ha Rhys Me Tas Tye 1 
Ee ee 
Ad a afb pP TP TD TD — 6a) 
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wherein u® is the -component of the permanent electric dipole moment of a molecule 
of species ‘i’, and od is its electric polarizability tensor. The first term in the foregoing 
expression accounts for dipole-dipole interaction (Keesom), whereas the subsequent ones 
are those of dipole-induced dipole interaction (Debye, Falkenhagen). Thetensor TOP account- 


ing for dipole-dipole interaction is defined as follows: 
(3) l 26 1 (3 2) 
Tap — — 5 {BTijaTijp — TüDap] » . 
ij 
where Ó,5 is the substitution tensor, unity if x =f and zero when «+f. 


By the method explained in the Appendix, substitution of (3.1) in eq. (2.6) yields the 


consecutive terms of the expansion (2.4)!) in the form: 
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where a; = 4a is the mean polarizability of a molecule of species "7", and the brackets < > 
denote the mean value dependent solely on the model of central forces, i.e. 


oo ul") (1; 5) 
—n\ _ =n Peas 
ty") = 4a [ r5 e M mdr. (3.4) 
0 


The expressions (3.3) are correct for a central forces potential u® (r; j) of arbitrary form. 
Assuming, for simplicity, Lennard-Jones (6—12) potentials, we have, by (2.11), 


3—n 


(i^ = Hua), (3.5) 


1) In computing the expressions (3.3), of all the terms due to the energy of inductive interaction between 
the molecules only the ones linear in the polarizability æ have been retained. 
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whence the second virial coefficient given by the expansion (2.4) can be written in the form 
O RGD : p 
pO IS BU a B, (3.6) 


wherein B® is given by eq. (2.9). 
The contribution to the second virial coefficient B®? from dipole-dipole interaction is 
of the form 


G TMG LGN 3 uou Ys 
i : Oe em 4 DE A A Zt. = 
dip—dip B 288 obe? He( yi) + 400 0$ e} ^1 (y) + 
29 oY; 
+ -940800 et cy 


The last contribution to B® in (3.3) results from inductive interaction between the 
molecules; if these have the axial symmetry, it assumes the form 


3 N | 
ans TPO esce nu — SIE in H = 
dip—inddipD Wore; y3 Leon + ioa) e (Yi) c 
B up ofl 
HG) Yay Le 2 HT a ô ) IE 5 
T 40 0,22, E (1 2 2 ia e Rave) Aj) 12 i) + 
29 ud ut; yS one 8 
| 134200 cde E ( +7 2 Mei + u (1 +7 à) | His (yi) + a (3.8) 
17 71, 
here, the quantities 
(i) G) 
1 i i «y — ay 
Mi) = 3 (aff + 20?) and dg = ee (3.9) 


lenote the mean polarizability and the anisotropy of electric polarizability of a molecule 
X species "7", respectively, o being its polarizability along the axis of symmetry, and 
po — that perpendicular to the axis. 

From (3.8), the anisotropy of polarizability of the molecule, dq, is seen to have no 
yearing on the second virial coefficient in the first approximation; it plays a part in higher 
nes only. ; 

If, in particular, the dipolar molecules are isotropically polarizable, ô=0, ey. (3.8) 
educes to 


T aN 5 F 
dip—ind dipB&) = — pr rm (xou + uia) m ya) + 
3.2. 4. 29 ud ub y 
Hay Gy X | MG) Gy H ? ae AT (3.10) 
+ 0 086 Pu09 + 1310959«$ Pis 09) 


4. Anisotropically polarizable quadrupole molecules 


In the case of a gas mixture consisting of anisotropic molecules which, though not 
ipolar, possess a permanent electric quadrupole moment, we have (see, Kielich 1960) 


1 ] jj 1, ; 2, j 2, j (2,5) a, 7) qd, D) G7) Gi) FE 4.1 
yg? =—5 O OT Pm n (ag 05507 Js Kap 9,8 Orn, ) Taya Then ior 
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i : Pre dee CET ; ! 
here, OË) is the tensor of the quadrupole moment of the molecule of species 7", given by 


02 3 3 2 e (3709 r9 — rá; ÖzB); (4.2) 
n 


with e? denoting the n-th electric charge of the molecule of species “7”, and 1(O — its radius 
vector. : 
The first term in eq. (4.1) accounts for the potential energy of quadrupole-quadrupole 
interaction, whereas the second is that of quadrupole-induced dipole interaction. The 
respective tensors T and ides describing these interactions are determined by eqs. 
(4.10) (see, Appendix). 
By eqs. (4.1) and (2.6), the consecutive terms in the second virial coefficient are obtained 


as follows (cf. eqs. (A.8) in SECUN 


BY hp E 5 RO + - 969 225) £u 
; IAN 
ee - ee 
Bi" —  4gpsTa C an 
8N {209 999 + OLOVOVPOPY rt». 
* 3538 35 k2T? : 
Dus z- 64N AET OT (4.3) 
945,373 - «8 dn POS . 


With respect to the foregoing expressions, the second virial coefficient of eq. (2.4) 
assumes the form 


pe» — BË Bo 4 ae Be (4.4). 


quadr-quadr quadr-ind dip 


In the case of molecules possessing the axial symmetry, and for a (6—12) Lennard-Jones 
potential, the contributions from quadrupole-quadrupole and quadrupole-induced dipole 
interaction are, herein, of the form 


MOHON | 180,99: 


us nando em = 480078} 7 E TET |^ ol) — 3430%¢ Sur Hisy ^y) T. H (4.5) i 
yer yeij 


R aN 
quadr-ind dipB) = (aaO + Oha) Hs (ya) — 


— 16oba: yt 


C), O PEETA 
— SPO (adat Obando) Hal) E (4.6) 
Or is the quadrupole moment of an axially symmetric molecule of species "7", which, by 
the general definition of eq. (4.2), is given by 
Oy = OF = 4 D O B22 — 7), (4.7) 
LI n 
where Z,,; is measured along the molecular axis, and 
o = — 269 = — 208, EL 


ANE 
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5. Multi-polar molecules 


Let us now consider the case of interacting molecules of a gas mixture presenting, 
quite generally, a dipole, quadrupole and octupole moment. The potential energy of multi- 
pole interaction between two unlike polar molecules is given by the following expression 


(see Kielich 1960): 


Bm DT +4 3t (u oP Lie On) TE + 


+5 BR — DOS 2965 + Mp un 1D) TO t I (5.1) 
wherein (20, is the tensor of the octupole moment of the molecule of species “i”, defined 


as (see, Buckingham 1959) 
ett pa en (rro — 72 (265, + 1009, 4 r3, 9). (5.2) 


apy — na! np! ny np” ya 
The terms in eq. (5.1) determine the potential energy of dipole-dipole (rz 2), dipole- 
quadrupole (r; T and, furthermore, dipole-octupole and quadrupole- pere (rg ?) inter- 
action, Pee poctively. 
Analogously, by (5.1), eq. (2.6) yields the following contributions to the second virial 
coefficient (see, formulas (A-9) in Appendix): 


BY = 0, 


3 uN Or dot 
BY = age (40 — a (50009 + 08 Ona uo) > — 


ET ER - 3.009,00, 4), 
BP m rg tO uPOGO RP A? Cr) + 
* ams uei 0901201208 + + 
Beer 03999999 99 9t? Gas, 
Eo E Ber vs + why ou ui uy) Cr > 


(5.3) 


On substitution of the foregoing expressions in eq. (2.4), the second. virial coefficient 


for polar gases is obtained in the form 
: "um " , 2 Gi 
BO = B 22 T dip—dip? E dip—quadr? Ho and Se dip—octup? z (5.4) 


where the first, second and fourth terms are given by eqs. (2.9), (3. 7) and (4.5), respectively. 
The two remaining ones, accounting for dipole-quadr upole and dipole-octupole interaction, 
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are, in the case of axially symmetrical molecules, given by the following expressions: 


dip— guadrD To — 1095 5c 1095.22. T uet: gu ud) D ij) = 
OijEij 


6050 yi | | 2456009» | 5:5 
DI Hl) + H, ie (5.5) 
350 "m ej Hysyij) «Sober (i j) + : 
Gj) aN 32708. B DRE ET eae 
dip-ortap US 14078 (ua) (7 OUG) 10(Y;) ur 
Opin ncn Y 
» 10008 a (o£) + Rouo) Has) + (5.6). 
wherein Qa is the octupole moment of the axially symmetrical molecule of species “i”: 
2a us = 0, = i D eg, | 5 gu "xs 3r2 i (5.7) 
with 
QU, =— 2003 =— 2085. (5.8) 


6. Applications to special cases 


We shall now proceed to apply the formulas for B(T) derived in sections 2—5 to some 
special cases. 

(1) One-component gas. Vor a gas consisting of a single species, the contributions to 
the second virial coefficient computed in sections 2—5 assume the form 


By = $a! NF(y), (6.1) 

aip—tipB =— E ME nU ges o M ah (6.2) 
dpa iB = — FOE H0) + 237, Maal) + e Ha) + 

1 2s d |i: P E m HOME | + a (6.3) 

quadr—quadr D. = — N | T(y)— ee Hy(y) m , (6.4) 

GARS RUE a {us m Ben kuh (6.5) 

dip-quadrD = — N sco — i Hal) — Se Hys( y) + s (6.6) 

toand ET Od — 00 (m 
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The expressions (6.1), (6.2), (6.4) and (6.6) are identical with the ones derived by 
3uckingham and Pople (1955). 

In an earlier paper (Kielich 1960), the second virial coefficient Bs of the molecular 
onstant of light scattering was computed for the molecular models considered in the fore- 
oing Sections; for molecules of constant polarizability, the isotropic contribution to B 
elated to B(T) by the formula 


sis 


y 4 
Bs =— nm, (6.8) 


vherein A is the light wavelength. 
(ti) Binary gas mixture. By eq. (2.1), the second virial coefficient of a binary gas mixture 
s of the form 
B= a4 BOD + 2 x x, BO? + pe. (6.9) 


n particular, the molecules of one component of the mixture being spherical and those 
f the other quadrupolar, the expressions derived in Sections 2—5 assume the form 


B= $n, NF(y1), (6.10) 
2 3a 0% 
B% = 3 mo} {Fld e Er ME (6.11) 


3 (2) Of) | | ee 6 0200522 
16 032£22V32 iv: 35 032699 


2 
B22) — 3 str Yoo) Du Hala) | Pe 


70%, 18 Ot y32 | 
onen — Ht ‘ 6.12 
320 019.3, Hiofyog) 343 O82. 15 23) (6.12) 
Eq. (6.11) was derived previously and applied to an argon-carbon dioxide mixture 
y Buckingham (1959). 
(iii) Ternary gas mixture. In the case of a gas mixture consisting of three unlike com- 
onents, the second virial coefficient has the form 


B(T) — ier F 22 D T pn ais 20D d" ADU T meh) (6.13) 


E.g. assume the molecules of component 1 to be spherical (%q) 4 0, jt — Og — 0), those 
f component 2 to be quadrupolar (u= 0, 6,5, #9), and those of component 3— dipolar 
4g) #0, Og = 0). The expressions derived in sections 2—5 for BD, B? and B®) once 
jore yield eqs. (6.10), (6.11) and (6.12), whereas those for BO), BED and BG? lead to the 


llowing formulas: 


2 ADU) | 
(18) — NiF a2 2.202 E ; (6.14) 
B 3 nots | (X13) 16 03561, Ys e(13) 
2 ADU) 
(28) — STNE Fi =. = Sie SA Man) — 
B 3 1023 vy (a3) 16 05565, V33 6( V23) 
30%) | E l 
re S H , (6.15) 
3208; £33 1463) sil 12 £53 sly. 23) ( 
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4 „4 4 
2 Gal) E 1 5 ) Ho Nee 
BR) == 3 noža N I (¥33) — BollsessY 2, D ar 400 ZW | d 5 (3, 13€ 33) 


29583 8 | IE NM E H 
: -[1+— ôg] H +... (¥33) + 
- 134400022685 7 3 a) 192083835 | °M 7 


3ula)y33 29u) y33 i | | 6.16 
gie DD LM er Vos) s obe (6.16). 
! 40068,:3, P 12 aU a5) * 940800 940800043233 Hast 33) 


The expressions derived in sections 2—6 for the second virial coefficient can be easily 
applied to more complex gas mixtures. We refrain here, however, from considering these 
special cases. 


APPENDIX 
The method of isotropic averageing of Vi; 


In order to compute the terms in B® as given by eq. (2.6), it is convenient to carry 
out integration over all possible configurations of the molecules first, i.e. to compute the 
mean isotropic values of the respective powers of their tensorial interaction energy: 


vg = E [ i (V(r... w "dodo. (A.1) 


On inspection of the various forms of the energy V; in sections 3 — 5, it will be seen 


that computation of BR implies averageing of the Sede products of the tensors My» 
Aap»? Ip and Q,,,. This can be effected as follows: Let T, ay... be a tensor of arbitrary order, 
given in a fixed coordinate system in three-dimensional space (X,), &—1, 2, 3, and Tow 
one given in a coordinate system in motion with respect to the first one, (X), o=1, 2, 3. The 
lensors are related by the transformation formula 

Ten 5 De OV Da c E (A.2) 
wherein, for rectangular coordinate systems (X,) and (X,), the transformation coefficients 
Osa» +++ have the meaning of cosines of the angles subtended by the respective coordinate 


axes of the systems. The directional cosines c, and @,, fulfill the conditions of normalization 
and orthogonality: 


Doaap = Ogg. — Cg Du = Oar: (A.3) 


Thus, the problem of averageing the tensor T,,,...reduces to that of averageing the 
corresponding functions of the directional cosines over all possible orientations of the system 
(X,) with respect to (X,). Tensorial summation over the indices c, t v,... leads to the fol- 
lowing mean values of the functions of the directional cosines: 


Oog xB = 3 O gp ar > 


ry 06 ae 36 (4 OapOve ae Gl Og py) 65,0 ME 
+ (40,05, — 04405, — 0,50,,,) 05,0, + (4 05405, — an y — Say Òga) omom , 


Veg xp 
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= WIE. ] 
D oq rp vy® 950,0 un — 570 1530,,0,,0; , + X30 0,,0 


B X50, 00 5,0 xa xa X05, OL, + dd T. Xg orð 22% on zs X0 5,0 euro i 000 


ap 


or ou jg 330,0, 0, = X40 PTS aa 3s 


2g 
ii X110 yO oO oA Fi: 120 990 rv = X130 5,0 uÔ Fs %149 0% oA vy xs 3150 ,,0 au vas ? (A.4) 
with the notation 
x; = 16a, — 5 (ag + a3 + Og + a5 + ag + a7) + 2 (ag + ag + aio + ayy + a + yg + 444 + 235), 


X3 = 16a, — 5 (ay + a3 + ag + ayo + Gia + ais) + 2 (a4 + a5 + ag + a; + ag + an + aq + 0134), 


Xs = 16a4 — 5 (ay + a4 + ay + du + ai + aa) + 2 (a4 + ag + ag + Gy + Gg + ai F ais + a) 
x, = 16a, — 5 (a4 + a5 + a + a49 taz + aia) + 2 (Ag + ag + ag + Gy + ag + an + Gy + di5)» 
x; = 16a; — 5 (a4 + a4 + ag + an + a + ais) + 2 (as + ag + dg + a; + dg + d49 + 033 + aig), 


Xy = l6ag — 5 (a4 + a4 + ag + ag + ay + is) + 2 (a3 + ag + a4 + ag + aio + au + d43 + 3), 
Hq = l6a, — 5 (a4 + ag + aio + y + ara + ara) + 2 (ag + ag + a4 + a5 + ag + a + ayy + as), 


Xg = 16ag — 5 (ay + ag + ag + ag + aio + t12) + 2 (a4 + ag + a4 + à + an + tis + a44 + a): 

4 = 16a, — 5 (ag + a4 F ae + Gg + an + ar) + 2 (ay Fat as + a, + 039 + ta H- tia + a) 

x10 = 16a, — 5 (dg + a4 + Gy + ag + a44 + ayy) + 2 (a4 + az + ag + ag + ay + ai + 43 + ai) 

X = 1624, I ta F Gz + ag + tio + 035) + 2 (a, F ast a, + ae + Gg Hat ait 034); 

= lOd — 5 (as +, + Gz + Gg + Qs + 044) + 2 (a4 + Gg + G4 + Gg + dy + 049 + Oy, + 235). 

4497 16043 9; Gy + ay + ag Ft. + d45) +2 T ag + as tu + Gg F 049 + 033 T 034); 

41,7 1605, — 5 (ag + dy t Gs = GoT 013 + t) + 2 (3 + a, + a, + a, + Gg + Gy + 033 ++ as) 

X15 = 160, — 5 (a5 + a5 + ag + au + trs + 234) + 2 (a4 + a3 +, + à, + Ag+ ag + ayo + ao), (A.5) 


nd 
(AR 59 p%y6%en lg = Os pO yey ag = Opn? 5e 
Ga = Oy aspe as = 0, ann per Ag = ben Say Bor 
a, = Od Ag = ÖgyOßedon: Ag == Say pn de> 
Ao = OgepyOdn> 4j = Sqn py de Ah, = CUNT PW yn 


dj, = 0989p ye: 14 = ÖgedmOßo» a> Sqn Oye ps (A.6) 


With the formulas (A.4), the mean isotropic values of the respective powers of the 
nergy Vy given in sections 3—5 are readily computed. 

Hence, if the interaction potential energy of a pair of dipolar molecules is given by 
q. (3.1), the following mean values result (in computing the terms accounting for the part 
f the energy due to inductive interaction, only the ones linear in œ have been retained) : 


ER 2 2 Gi) par) 
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ri (21 (aout + iio) Boc + 


9 AG) , GG), G GG) ON cee DES OUEST ee 
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Dey One APA TAT TN TE) 
Ty aPuP ul?) uct IRB: 12 T Ti, l 
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TS utut) sj) rp GI) rp G3) rp G3) nr qr) 
Ke = Pte {7 OT fee (dee do Tin Ti + 
POLO TO TOT OT TG RETT teed T (A.7) 
Likewise, for quadrupole molecules, eq. (4.1) yields 


Vg = — & (009094 00092.) TOTS, 
Vi ais OROORO T9. TS + 
+ as (089000009 + ALLAN TET T: 
V3 = — d$. 009999ot et eO TG TG) TG X. (A.8) 


For the case of polar molecules we have, by eqs. (5.1) and (A.4), 


V, — 0, 
V} = Vit TT + i (ORO + OEL TETE + 
"m Ta BLR 29), "up 705200100200) 3 
+30 t 2 s sy) TO T ? 
V§ = — hp uPOQORup{TOT OTS, + 2 TOT OTD — 


RT (5,490 (GO) ad G s 3 ».Gg 3 es i) oe 
1575 {uz / D [o 0 "e x 9% t 02) gout nae TES i T aie 


—_ 16 (9) go (9G 9 qqQoOY9GO m E Q7) a 
99225 II às bn en I Ty TE 


TURION in mn man mas ia RT EE. 
RR EEE ese 


1J 
Aes ee CAG mq NN EM 
gos Vail Ib s MG) + UDR aaa My) Loe Lon Lee Tome + s (A.9) 
. The tensors T®, TO and T 9 appearing in the foregoing expressions and characteriz- 
ing dipole-dipole, dipole-quadrupole and quadrupole-quadrupole interaction, respectively, 
are defined as follows (see, Pople 1952): 
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+ rg igo.) + 15 (05,0, 4- 0,0, + 0,,0,)]. (A.10) 


With the foregoing definitions, we have, for the various products of the tensors TO 
DO) and T™ appearing in eqs. (A.7) — (A.9), 


or otvo 
GG) _ 6,76 Gp Gi mp anm o = 
2. Re = 67 pe ESAS ON I Ue = lar; et 


TOT qn Ih p po 566r. 18 IE, DTH) — 90r 29 


TPTÜTS- cU TPTPTS,-— 265% 
Dian OA T ese 10807, %, 

TOTOT OTS = DI lores, 

TË TO TO, =— 194407; 55. (A.1) 


In applying the expressions of eqs. (A.7) — (A.9) to molecules having the axial symme- 
ry, use is made of the following identities: 


3Aapllallf— Xa gg = OHO", 
apap = 3x00, — 0,,0,, = 40°, 
Hals = UO, Ogg gO py = OP, 
Halt ply Qapy= WO, Ban apy = 382°. (A.12) 
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TWO-PHASE FERROELECTRIC SYSTEMS 
(I. BARIUM TITANATE-METAL) 


By Mns. KAROLINE LEIBLER AND Mr. Wojciech BRANSKI 
Polytechnic High School, Chair B of General Physics, Warsaw 
(Received January 2, 1961) 


In this research the dielectric properties of the two-phase ferroelectric-metal systems 
have been investigated. It has been found that a considerable increase of the dielectric constant 
can be obtained corresponding to the increase of the concentration of the metal admixture. The 
shape of the increase curve depends upon the kind of admixture and the size and shape of its 
grains. The Curie temperature does not change within the limits of experimental errors. However, 
changes in the Curie-Weiss temperature can be observed, accompanied by a change of the 
Curie-Weiss constant, ie. by a change in the slope of the paraelectric part of the curve 1/e,,. 


x 


Investigations of the dielectric properties of the two-phase dispersed systems has revealed 
he dependence of the dielectric constant of the system upon: 

l. The dielectric constant of the medium 

2. The voluminal concentration of the admixture 
nd 3. The size and the shape of the dispersed particles. 

In the dielectric-metal systems a considerable increase of the dielectric constant is 
btained corresponding to the increase in the concentration of the metal, whose dielectric 
onstant can be taken as e= oo. At a certain critical value of concentration the specimen 
ecomes a conductor. In previous research work such dielectric plastics as paraffin and poly- 
yrene have been used for the medium (Leibler 1959). In this research the dielectric proper- ` 
es of the two-phase ferroelectric systems have been investigated. Barium titanate was 
sed as the medium with the admixture of these powdered metals: silver, nickel and iron. 


T. The experimental method of measurement 
1. The technology of the preparation of specimens 


he metal powder with its grains irregular in shape, was passed through a 0.1 mm sieve. 

Since it was necessary to bake the barium titanate ceramics at the high temperature 
° 1350°, the greatest difficulty experienced in the preparation of the barium-titanate and 
etal system was the prevention of the oxidation of the metal. Thus the specimens, after 
ey had first been compressed under a pressure of. 3.5 thousand atm. had to be baked 
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in the vacuum oven under a pressure of 107? to 107* mmHg. The specimens were baked | 
twice, since after the first baking they were still rather brittle. The plates were sprayed with | 
silver. 

A defect of the method consisted in the necessity of a rapid cooling of the oven. 

The dielectric constants of the specimens obtained by this method were not optimum | 


ones. They amounted to 1000 for pure titanate. 


2. Measurement of the dielectric constant 


The measurements of the dielectrie constant and those of the dielectric phase angle were 
taken by the bridge method at frequencies of 1 MHz and 200 KHz. The measuring condenser 
was of the parallel plate type with micrometric control. The accuracy of measuring was of 
the order of 0.5 pF, which for high dielectrie constants gives an accuracy of over 1% for the 
real component of the dielectric constant. The dielectric phase angle was established. with 
an error not exceeding 5%. 

'The measurements of temperature were taken in the condenser placed in the ultrather- 
mostat. The accuracy of temperature measurement was 0.2^ C. 


3. Determination of the voluminal concentration. 


'The voluminal concentration of metal was determined by measuring the masses of the two 
components and the specific density of the specimens. The measuring error, however, must 
be estimated at 6—8%. 


4. Examination of the structural condition of the systems 


For the purpose of establishing the structural condition of the system and of making sure 
that a mixture of barium titanate and metal powder has really been obtained the specimens 
were examined by the X-ray method. This X-ray analysis confirmed the existence, in the 
specimen, of barium titanate and metal. 


II. Measuring results 


|. Investigation of the dielectric constant. at the fixed temperature of 20°C has shown 
an increase in the dieleetrie constant corresponding to the increase in the concentration 
of the metal as is the case in systems with ordinary dielectrics. At certain values of the 
concentration. characteristic of the kind of admixture the specimen becomes a conductor. 
In fig. I are shown the curves for the systems: 


l. barium titanate-nickel 


2. barium titanate-silver 


2. Investigation of the dielectric constant as a function of temperature has shown 
a normal increase of the dielectrie constant at temperatures near the Curie point. The Curie 
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point does not shift visibly. In fig. II, 1 are shown the curves as affected by temperature: 
&=/(t), for the system barium titanate-silver. Fig. II, 2 illustrates the same for the 
system barium titanate-nickel. In fig. III 1, 2 changes can be seen in the Cue 
temperature accompanied, as the experiment has shown, by a change of the Curie-Weiss 
constant, i.e. by a change in the slope of the paraelectrie part of the curve. 
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Fig. I. — Relative dielectric constant &,, versus Ag and Ni concentration in BaTiO, dashed lines — according 
the theoretical formula 


III. Experimental data and theory 


Theoretical evaluation of the dielectric constant of the system as a function of the volu- 
ninal concentration of the metal powder is based on the use of a series of simplifying 
nodels. 

In an earlier research it has been shown that Piekara’s (Piekara 1932) formula gives 
he best approximation for the dielectric metal systems. Half-empirical formulae have also 
jeen arrived at containing a modifying coefficient corresponding to the critical concentration 
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at which the specimen becomes a conductor. As can be seen from fig. I, the obtained results 
show that the best approximation is achieved with Odolewski’s (Odolewski 1957) modified | 


formula 


1 
ëz . where k = E 
cr 


£,, — dielectric constant of the mixture, £, — d.c. of the medium and is the concentration 


of the metal at which the dielectric becomes a conductor. 
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Fig. II, 1. — Relative dielectric constant Ew versus temperature for BaTiO,-Ag 
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IV. Conclusions 
It now appears that two-phase dielectric can be obtained with high dielectrie constants 


if mixtures ferroelectric-metal are used. The dielectric properties of such mixtures depend 
upon the voluminal concentration of the admixture and the shape of its grains. 
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Fig. II, 2. — Relative dielectric constant &,, versus temperature for BaTiO,-Ni two-phase system 


The admixture does not affect the Curie temperature, but the Curie-Weiss constant 
changes within fairly wide limits. At a certain critical value of the concentration of the 


admixture the material becomes a conductor. 
Further research is being carried out with a view to improving the technology of the 
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Fig. III, 1. — 1/e,, versus temperature for BaTiO,-Ag two phase-system 
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Fig. II, 2. — 1/e,, versus temperature for BaTiO,-Ni two phase-system 
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two-phase ferroelectrics, to finding how this technology affects the constants of the system, 
and to learning the ageing process of the two-phase ferroelectrics. — . 

It is hoped that this research work will result in the obtaining of a new group of ferro- 
electric materials with properties to suit requirements. 
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ABSORPTION IN THALLIUM BROMIDE VAPOUR 


By E. OsraAszEWICZ 
Evening School of Engineering at Bialystok 
(Received December 23, 1959; revised manuscript received November 23, 1960) 


Absorption in thallium bromide vapour was investigated versus 1. the vapour density, 
2. the temperature, and 
3. admixtures of inert 

gases. 

Measurements were carried out with a grid spectrograph of grid curvature radius 3 m. 
The light source consisted of a hydrogen lamp. The time of irradiation was 8 hours. Absorption 
was found to depend on the vapour density and to be independent of the temperature of the 
absorption tube at constant pressure and of admixtures of foreign inert gases, and to consist 
of continuous absorption in the ultraviolet range and of bands situated in the long-wave region. 


. . . . . 4 
Isotope splitting due to the presence of two bromine isotopes yielded the magnitude of o — 52 ; 
M» 
the latter was compared with the figures obtained theoretically. From the system of bands, 
the heat of dissociation in the normal and excited states was computed, and the thallium bromide 
molecule was shown to present an atomic compound when in the state of vapour as it then disso- 


ciates optically into a thallium atom in the normal state and an excited atom of bromine. 


1. Ionic and atomic compounds 


With respect to the chemical bondings they present, molecules fall in two classes: 
onic and atomic (Franck et al. 1927). Ionic compounds, when in the normal state, dissociate 
‚ptically into ions, whereas, when in the excited state, they yield two normal atoms. On the 
ther hand, atomic compounds; when in the normal state, dissociate optically into normal 
toms and, when in the excited state, into one normal and one excited atom. Ionic compounds 
re exemplified by the alkali metal halides (Sommermayer 1929). Characteristic examples 
f atomic compounds are provided by the halides of borium, aluminium, gallium, indium 
nd thallium. Spectral analysis of its absorption provides the clue for classifying a molecule 
s belonging to the one or the other group. The absorption spectra of thallium bromide 
apour were investigated by Butkow (1929), whose results, however, are incomplete. 

It was the aim of the present investigation to fill in the void left by that paper and to 
rove that the thallium bromide molecule in the gaseous state presents an atomic com- 
ound. 
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2. Experimental Part 


Spectrographically pure thallium bromide salt is placed in the quartz vessel 4 connected 
to quartz tubing B of 30 cm length and 2 em diameter sealed off at high vacuum. The tubing 
is disposed within an electric oven. with plane quartz windows, thermally stabilized. The 
vessel could be heated with a separate small oven (Fig. 1). The apparatus was arranged. 
according to the diagram of Fig. 2. A hydrogen lamp provided the continuous light source, 


Fig. 1. Absorption tube in oven. 4 — vessel, B — absorption tube, C — oven, D — auxiliary heater (oven). 
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Fig. 2. Diagram showing arrangement of apparatus, Q — light source, L,, L, — quartz condenser lenses, 


B — absorption tube, S — spectrograph slit. 


A medium quartz Hilger spectrograph and a glass Fuess spectrograph were used for the 
investigation, together with a grid Hilger spectrograph for separating the isotopes. Irradiation 
ranged from 30 minutes to 8 hours. The wave numbers of the band heads measured were 
reduced to vacum (Kayser 1925). 


3. Dependence of absorption on vapour density 


The temperature of the absorption tube was kept constant at 1000^C, whereas that 
of the vessel was varied. The dependence of the vapour pressure on the temperature was 
computed from eq. (3.1) (Wartenberg and Bosse 1922): 


25164 atm. gr: 
log p CREE LED i 4- 5.059 atm (3.1) 


wherein pa, is the pressure in atmospheres, and T — the Kelvin temperature. 
Appreciable absorption requires at least 300°C and consists of the following regions: 
à. lhe system of bands G, 
b. three regions of continuous absorption situated in the ultraviolet range, 
€. a system of diffuse bands. 


p 
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a. Band system C 


This contains bands that are shaded towards the red end of the spectrum and present 
heads that are sharp and intense. Their sharpness decreases rapidly towards the red end. 
Table 1 Bows the coordination of the bands relating to the T1,55Brz9 molecule in Deslandres 
table. From Table 1, the graphs of the potential energy U versus the distance r separating 
the nuclei in the normal and first excited state of the molecule are found to have the shape 


shown in Fig. 3. 
LU 


Fig. 3. Potential energy graphs of thallium bromide molecule, for band system C. N — graph in normal state 


P — in excited state. E,; — energy of electron transition, E,, — excitation energy of bromine atom 


D” — dissociation work in the fundamental state, D’ — dissociation work in the excited state. 


b. Regions of continuous absorption 


Towards the short-wave end, the band system C is contiguous to three regions of contin- 
ious absorption and one narrow band having very diffuse edges and an absorption maximum 
t 3384 À. The first of these regions ranges from 3340 to 3310, the second — from 2595 
o 2450 A, whereas the third one lies within the range of shortest wavelengths. This is the 
icture presented by the continuous absorption at the vapour density corresponding to 
00°C. As the vapour density increases, the absorption, both continuous and band-form, 


istends in either direction. 
c. System of diffuse bands 


Another system of bands, lying in the region of 3910—4627 A, approaches the system 
‘from the long-wave end. The measurements are assembled in Table 2. This system differs 
'om C radically as to the intensity and character of the bands. The latter are of very low inten- 
ty and diffuse, without shading. 


4. Effect of foreign inert gases 


Measurements were carried out with two tubes filled with a mixture of thallium bromide 
ipour and argon, at the pressures of 2 mm and 190 mm Hg. The absorption of the vapour 
as investigated versus the pressure and temperature. Photographs failed to detect an effect 
f argon on the intensity distribution and general character of the bands. 
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TABLE 2 
System of diffuse bands 
No. | Zin A |vin cm” No. Ain A |vincm! No. Ain A |vincm^! 
— E um p i 
1 3909.05 | 25574.43 20 4132.50 | 24191.62 39 4325.00 | 23114.90 
2 3912.00 | 25555.16 21 4140.03 | 24147.64 40 4336.70 | 23052.54 
3 3923.00 | 25483.50 22 4151.20 | 24082.65 41 4346.80 | 22998.98 
4 3936.88 | 25393.66 23 4161.50 | 24023.04 42 4358.00 | 22939.88 
5 3940.40 | 25370.98 24 4170.50 | 23971.20 43 4379.90 | 22625.18 
6 3951.38 | 25300.48 25 4181.70 | 23907.00 44 4400.80 | 22716.78 
7 3963.60 | 25222.47 26 4191.00 | 23853.95 45 4413.30 | 22652.44 
8 3967.48 | 25197.81 27 4200.62 | 23799.76 46 4433.75 | 22547.96 
9 3979.00 | 25124.86 28 4213.00 | 23729.40 47 4447.27 | 22479.40 
10 3990.90 | 25049.94 29 4221.70 | 23680.50 48 4466.61 | 22382.80 
11 3995.08 | 25023.73 30 4230.00 | 23634.03 49 4480.37 | 22313.34 
iby 4014.09 | 24905.22 31 4242.40 | 23565.52 50. 4501.57 | 22208.26 
13 4023.21 | 24848.78 32 4251.30 | 23515.60 51 4523.00 | 22103.04 - 
14 4043.02 | 24727.02 33 4261.00 | 23462.10 2 4536.92 | 22035.22 
15 4051.50 | 24675.27 34 4272.20 | 23400.58 53 4557.10 | 21937.65 
16 4073.80 | 24540.20 35 4282.93 | 23341.96 54 4572.12 | 21865.58 
17 4080.90 | 24497.50 36 4293.00 | 23287.20 55 4592.89 | 21766.70 
18 4102.60 | 24367.93 37 4308.48 | 23225.05 56 4608.90 | 21691.75 
19 4110.45 | 24321.40 38 4314.90 | 23169.00 57 4627.50 | 21603.90 


5. Effect of the temperature 


For investigating the temperature dependence of the absorption, the vessel was heated 
to 400°C, while the tube was raised to a still higher temperature, when the two parts were 
separated by melting the quartz and sealing off. Thus, vapour at the pressure corresponding 


34534 


3384 À 


a) 


3453 Å 
3384 4 


b) 


Fig. 4a and 4b. Dependence of absorption on temperature of tube 
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>) a temperature of 400°C was obtained in the tube. Absorption was measured at 650° 
nd 1250°C. From the microphotograms of Figs. 4a and 4b, the temperature cannot be 
üd to affect the absorption. 


6. Isotope effect 


Double band heads of equal intensity appeared on the plate. These were due to the 
otopes of bromine. It seems reasonable that T1,95;, which is the more numerous, is effec- 
ve in producing the spectra. In the case of the (0,0) band, isotope splitting vanishes 
l'able 3). Towards greater distances from the origin, splitting increases. From the isotope 
plitting »,—v, as appearing on the plate, the value of o can be computed. using the formula 
Herzberg 1939) 

yy — *, — (p — 1) (Wy! — Wav") (6.1) 


TABLE 3 


i Splitting of band heads of molecules Thy; Br;, and Tlj; Br gı- 


; » Wave'number for Tlos Br,, and | Isotope splitting, in 
x d Tl; Bra, molecules, in cm ! cm 04—1 
0 3 28572.0—28566.4 6.2 0.0108 
0 4 28383.4—28376.2 7.2 0.0094 
0 5 28192.4—28184.0 8.4 0.0089 
0 6 28005.6—2 1994. 1 11.5 0.0101 
0 7 27818.6—27805.8 12.8 0.0097 
0 8 27634.4—27619.5 14.9 0.0100 
1 3 28668.9—28664.0 4.9 0.0101 
1 4 28479.9—28437.0 6.9 0.0101 
1 5 28292.2—28283.3 8.9 0.0103 
1 6 28105.0—28095.5 9,5 0.0101 
1 7 27918.8—27907.2 16.6 0.0094: 
2 if 28005.6—27994. 1 11.5 0.0099 
2 8 27818.6—27805.8 12.8 0.0096 
2 9 27634.4—27619.5 14.9 0.0098 
2 10 27449.6— 27433.4 16.2 0.0094. 
2 11 27264.0—27246.1 17.9 0.0095 
2 12 27080.1—27060.5 19.6 0.0095 
3 Tu 27337.0—27320.5 16.5 0.0099 
3 12 27154.1—27135.7 18.4 0.0092 
3 13 26971.3—26951.4 19.9 0.0093 
3 14 26791.3—26769.4 21.9 0.0094 
3 15 26612.2—26587.5 24.7 0.0099 
3 16 26432.2—206406.5 25.7 0.0097 
3 17 26254.5—26228.3 26.2 0.0092 
3 18 26072.3—26045.8 26.5 0.0089 
3 19 25896.7—25867.0 29.7 0.0092 
average Q4—1 — 0.0097 
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with IF, and Wo standing for the oscillation quanta in the excited and normal states, and 


v, v" — for the corresponding oscillation quantum numbers, whereas 
o=1/ 2 (6.2) 
Ha 


and fy, My are the reduced masses of the Tl 9;Brg; and Tly9;Brzg molecules, respectively. 
Denoting o as computed from eq. (6.1) by o, and as obtained from (6.2) by o, we have 
0, — 1.0090 and o, = 1.0097, i.e. og — 0, = 0.0007. The divergence is seen to be tolerably small, 
x > i ® E 7 * LI L * 

amounting to about 0.0795. Hence, the experimental results are in satisfactorily good agre- 


ement with the theory. 


7. Evaluation of energy of dissociation 


From Fig. 3, we have 


D” =E,,+D'—E,, (0.1 
wherein: E, is the energy of the electron transition (0, 0), 
E „ — the excitation energy of a bromine atom, 
D” — the energy of dissociation in the normal state, and 
D' — that in the excited state. 


Eq. (7.1) yields the value of D”, all other quantities being known, namely: E, from plate 
measurements amounts to 29143 em!, D' obtained by linear extrapolation to the point 
of convergence is 566.5 cm-!, and E, as computed by Turner (1926) amounts to 3700 em-!. 
Hence, D” =26009.6 cm^lor 73.9 kcal per mole. The value obtained by chemical methods 
(Butkow and Terenin 1927) was 73.0 +0.5 kcal per mole. This proves the thallium bromide 
molecule in the excited state to represent an atomie compound, as in this state it dissociates 


into a normal atom of thallium and an excited atom of bromine. 


8. Conclusions 


The present investigation completes that of Butkow. By raising the pressure in the tube, 
considerably vaster experimental results were obtained. From the dissociation energies in. 
the normal and excited states as evaluated by the present author, the thallium bromide 
molecule is found to represent an atomie compound. No thallium isotopes were detected, 
neither was any rotational structure obtained even when using a grid spectrograph. The 
temperature effect failed to yield the results expected. The presence of argon was found 
to have no effect on the absorption spectrum. 

The present investigation was carried out during the years 1936 and 1937 at the Institute 
of Experimental Physies of the King Stefan Batory University at Vilna, under the direction 
of Professor Józef Patkowski. 


9. Review of subsequent litterature 


Subsequently, investigation of absorption in thallium bromide vapour was carried out 
by Howell and Coulson (1941), and by Tiruvenganna Rao (two papers, 1949). Howell and 
Coulson applied conditions similar to those of the present paper, obtaining similar or iden- 
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ical results. Tiruvenganna Rao, using a grid spectrograph of dispersion 5.7 Ä/mm, separated 
he isotopes of thallium in the spectral region about 4200 Ä; in his second paper, he comple- 
ed Deslandres table for the Tl;;Brg; molecule to 4050 A by detecting some new bands. 
None of the foregoing authors succeeded in separating rotational structure. No paper dealing 
vith the subject appeared after 1949. 
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PHOTO-STIMULATED EMISSION OF EXO-ELECTRONS FROM THE 
SURFACE LAYER OF AMPHOTERIC METALS REACTING WITH 
SODIUM HYDROXIDE OR POTASSIUM HYDROXIDE 


By B. SUJAK AND J. WAWRZYNIAK 


Institute of Physics, Polish Academy of Sciences, Wroclaw, Chair of Experimental Physics of the University, 


Wroclaw 
(Received December 7, 1960) 


When reacting with a thin layer of hydroxide sprayed onto their surface, the amphoteric 
metals Al, Zn, Sn and Pb have been found to emit negative charge carriers into an atmosphere 
of air, if the reacting surface is simultaneously irradiated with white light (Amin ~3600 À). 

The intensity versus time graphs of the emission observed show an increase during the 
first phase of the experiment and, subsequently, an approximately exponential decrease. 

'The emission observed is related to the momentary shift of the effective long-wave boundary 
of the photo-effect during the reaction, until the wavelengths of the light waves of the incident 
beam are attained. The long-wave boundary of this photo-effeet is displaced towards shorter 
wavelenghts as we proceed from Al through Zn and Sn to Pb, i.e. as the affinity of the amphoteric 
metal for the hydroxide diminishes. 

The photo-stimulated chemo-emission of exo-electrons was detected and measured with 


the aid of a point counter with steering grid, in atmospheric air. 


Introduction 


For over 70 years it has been known that e.g. oxidation of phosphorus in air is accompa- 
nied by ionisation of the surrounding atmospheric air. On the other hand, various metals 
in contact with active gases or vapour of chemically active liquids emit electrons and/or 
negative or positive ions. The reader is referred to the detailed review of the problem as 
published recently by Roich and Yarpovietsky (1959). 

It will be remembered that the phenomena of emission of electrons of very low density 
;ecurring exclusively subsequent to excitation and vanishing gradually with the time are 
known jointly as exo-electron emission. 

A review of the opinions expressed during a discussion dealing with the nomenclature 
und the state of information on the subject as of 1956 will be found in the Report on the 
Innsbruck Conference (Tagungsbericht ""Exoelektronen" 1957). 

All known exo-electron emission effects are at present generally brought under one 
f the two following headings with respect to their basic mechanism (see, e.g. Sujak 1959 a): 

(463) 
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1) Electron emission from hitherto electron-invested electron traps (crystal lattice 
defects) of the emitter; 

2) Electron (or ion) emission accompanying chemical reactions of the emitter. 

It should be stressed (for earlier mention, see Sujak 1959 b) that cases can occur wherein 
the mechanism of emission is one resulting from superposition of the two foregoing processes, 
thus making a decision difficult. Moreover, effects related to the work of the exo-electron 
detector can also oceur which, although resembling exo-electron emission in their time- 
dependence, are due e.g. to diffusion of gases or water vapour volatilizing from the specimen 
(Sujak 1958). 

The present paper deals with photo-stimulated chemo-emission of exo-electrons accom- 
panying chemical reactions between the deformed surfaces of amphoteric metals and con- 
centrated sodium or potassium hydroxide. No data on the subject are available hitherto. 
Indeed, it is a new finding that metals such as Al, Zn, Sn, technical cadmium, brass or 
soldering alloys containing Sn, whose surface had been pre-processed, show an increased 
emission of exo-electrons on illumination with white light when their surfaces are reacting 
with sodium or potassium hydroxide (Wawrzyniak and Sujak 1960). 


Device 


The measuring device constisted of a point counter (diameter 1.4 em) closed with a metal 
grid, in air. By applying an accelerating or retarding potential between the specimen and 
the grid, the sign of the particles emitted could be determined. The light source consisted 


Fig. l. Diagram of measuring device: 1 — counter cathode, 2 — counter anode, 3 — plexiglass lense, 4 — 

grid, 5 — specimen, 6 — interference filter, 7 — bulb, 8 — high voltage source, 9 — amplifier, 10 — oscil- 

loscope, 11 — amplitude discriminator, 12 — pulse rate counter, 13 — recording milliamperemetre, 14 — 
specimen polarizer. 


of a tungsten 6 V/50 W bulb with glass condenser. Interference filters were interposed be- 
tween source and specimen for investigating the spectral distribution of the exo-electron 
emission. 

The electronic recording device consisted of a wide-band amplifier, an oscilloscope 
for visual observation of the counter pulses, an amplitude discriminator and rate meter 
coupled to the recording device. 


The device is shown schematically in Fig. 1. 


465 


Experimental Results 


The newly abraded surface of certain metals is found to be the source of photo -stimula- 
ted. exo-electrons emitted during illumination with visible light. Especially intense sources 
are aluminium and magnesium. The intensity of photo-stimulated exo-electron emission 
in aluminium, zinc, tin, technical cadmium, brass, and soldering alloys (containing tin), 
can be appreciably enhanced by spraying the surface with a thin layer of NaOH or KOH (see, 
Wawrzyniak and Sujak 1960). 

Some other metals, such as: magnesium, copper, iron, bismuth, antimony, nickel, 
mercury, molybdenium and tungsten failed to exhibit this property under identical condi- 
tions. Originally (Wawrzyniak and Sujak 1960), lead had been found to show no increase 
in emission on being sprayed with a thin layer of hydroxide, notwithstanding the fact that 
it is a weakly amphoteric metal. More recent investigations, however, proved that lead, too, 
exhibited the effect, provided the glass condenser was removed from the path of the light 
beam incident on the specimen. Obviously the condenser had absorbed the light quanta of 
the energy required. 

Photo-stimulated emission of exo-electrons can also be obtained by spraying a thin 
layer of hydroxide on aluminium already covered with a natural layer of oxide. On the 
other hand, no such emission could be obtained from aluminium covered with a thick, 
thermally formed oxide layer in a specimen that had been heated to about 500°C in atmo- 
spherie air. According to Hunter and Fowle (1956), aluminium heated to this temperature 
should be covered with a layer of oxide about 140 À thick. 

Neither was it possible to excite aluminium oxide (in the form of y —Al,O,) to exo- 
electron emission by spraying it with NaOH or KOH. The same is true of ZnO. It should 
be stressed that, even in the case of NaOH and KOH reacting with Al, Zn. 5n and Pb, no 
emission of exo-electrons is detected unless the surface is illuminated with a light beam. 

Further investigation dealt chiefly with emission from aluminium, zinc and tin, as the 
intensity here is a relatively high one, so that the emission is easily accessible to detection 
when white light is used. The surface of the specimens was always abraded mechanically, 
steel being used in the process. The hydroxide was brought into contact with surface either 
by using cotton-wool, a small brush, or by spraying from a gun. NaOH and KOH solutions 
ranging from In to saturated at room temperature were used. 

The intensity of emission and its time dependence varied strongly according to the 
nanner of spraying the hydroxide onto the surface of the specimen. In most cases, following 
ın original increase in intensity, the latter gradually decreased to zero. Hence, the intensity 
versus time graphs exhibit a maximum. In shape, the graphs strongly resemble those obtained 
ny Lohff (1956) for zinc oxidated in an atmosphere of rarefied oxygen. Generally, the time 
luring which the intensity of emission increases, and that during which it vanishes, differ 
rom one specimen to another. These differences reside in the difficulty of mechanical sur- 
ace processing (Mader and Sujak 1960) and of reproducing the thickness and degree of 
iomogeneity of the hydroxide layer. 

Fig. 2 shows two typical intensity versus time graphs for two Zn specimens sprayed 
vith a layer of NaOH. Graph 1 represents a case of relatively slow increase and decay of 
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the intensity, whereas graph 2 exemplifies one of rapid increase and decay. In addition to the 
graphs of Fig. 2, others showing a high rate of increase foliowed by very slow decay were 
obtained, whereas the inverse case of slow increase and rapid decay of the emission is a rare 
one. 

The intensity of emission corresponding to the maximum on the graph varied from 
specimen to specimen of the same metal. However, the various metals can be said to differ 
in the intensity of their emission. Thus, the ones investigated were found to present the 
following sequence of decreasing intensity: Al, Sn, Zn, Pb. The inverse sequence is obtained 


————————— 
time 


Fig. 2. Emission intensity N/t from two Zn specimens versus the time. NaOH layer obtained by spraying 


on applying the criterion of the time required for attaining the maximum on the intensity 
versus time graph. In general, aluminium exhibits a very short period of increasing intensity, 
whereas e.g. zinc requires longer time, although cases of the inverse cannot be ruled out. 
For this reason, a great number of measurements had to be carried out, in order to obtain 
statistical insight. The time during which the intensity of emission increased ranged to 
some tens of seconds, whereas that required for total decay amounted to as little as several 
seconds or as much as some tens of minutes. 

Efforts to establish a quantitative relationship between the intensity of emission and 
the thickness of the hydroxide layer proved unfruitful. The same is true with regard to the 
choice an concentriaton of the hydroxide used. In the extreme cases, Le. for a very 
"thin? hydroxide layer, the emission exhibits low intensity and. decays rapidly. As the 
thickness of the hydroxide layer is increased, the intensity corresponding 


emission on the graph rises, as also the time if its decay. For Schick 


to maximum 


Ä ayers produced e.g. 
by immersion of the specimen in the hydroxide, no emission can be detected. This is 


probably due to the fact that, as the source of emission is situated in the layer of the 
metal contacting the hydroxide, a homogeneous, thick layer of the latter transmits no 
exo-electrons to the detecting device. 

Dy applying an electric field decelerating particles of negative charge, it can be shown 
that those participating in the present effect are indeed of this kind. Inasmuch as the primary 
particles are photo-electrons, a counter surrounded by atmospheric air at normal pressure 
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records negative oxygen ions, since the electrons possess a very short free path and are very 
soon captured by the electro-negative oxygen (cf. also Hrbek and Vlasakova 1957). 

By using appropriate Zeiss interference filters, the relative spectral distribution of 
the intensity was obtained. To a certain degree, the time variations of the intensity could 
be rendered negligible by changing the filters rapidly several times. The intensity of emission 
reduced to constant incident energy. N/t/E, is plotted in Fig. 3 versus the light wave-length. 
The intensity is seen to fall as the wavelength increases. The long-wave limit of the photo- 
effect shifts towards shorter wavelengths as one proceeds from Al to Zn and to Sn. 

After the reaction between the hydroxide and the metals Al, Zn, Sn and Pb has taken 
place (even in complete darkness), illumination with white light fails to produce emission 
of exo-electrons. 

In order to make sure that the property of emitting exo-electrons is not related to an 
arbitrary exothermal reaction, the metals Al, Zn, Sn, and Pb, and, moreover, Sb and Bi 
were subjected to the action of the acids HCl, H;PO,, HNO,, HClO, and H,50,. In no 
case, however, did any photo-stimulated emission of exo-electrons above the average counter 
background occur. The foregoing statement holds in the case of white light as used when 
investigating the effect of the hydroxides NaOH and KOH. In some cases, when the reaction 
between metal and acid proceeded very turbulently, disturbances occurred in the counter 
background due to tiny drops of the acid penetratinginto the counter or to vapour of the 
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Fig. 3. Relative spectral distribution of exo-electron emission 


icid producing work pulses therein (Sujak, in print). Such disturbances, however, could 
be easily identified by extinguishing the light, by visual observation of the pulse forms 
yn the oscilloscope screen, and by varying the voltage between counter grid and specimen. 


Discusston 


As already stated, the foregoing emission of exo-electrons and/or its increase due to 
he action of hydroxides could be detected only in Al, Zn, Sn, Pb, technical cadmium, 
ring alloys and brass. No such effect could be detected in Mg, Cu, Fe, Ni, Sb, Bi, Hg, Mo,W. 


solde 
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Chemical analysis proved the cadmium to contain admixtures of zinc, as expected. 
This would restrict photo-stimulated emission of exo-electrons as due to chemical reactions 
with the hydroxides NaOH or KOH to the metals Al, Zn, Sn and Pb, and to alloys containing 
these metals. 

Magnesium, like aluminium. yields a high intensity of photo-stimulated exo-electron 
emission immediately after abrasion of its surface. However, if the magnesium surface 
emitting exo-electrons is acted on by hydroxide e.g. in the form of a layer consisting of 
a great number of tiny drops, the intensity of the exo-electrons emitted. is found to fall 
immediately. Aluminium, in contrast, exhibits a steep increase in exo-electron emission 
under the effect of hydroxide. The decrease in intensity as observed in magnesium can 
be explained by the absorptive effect of the hydroxide which, to an extent, covers up the 
emitting surface. 

The absence of excitation to emission in the oxides Al,O, and ZnO would seem to 
prove that such emission requires the presence of both hydroxide and amphoterie metal, 
or rather of their direct contact. 

Indeed, Al, Zn, Sn and Pb are amphoteric metals which, on reacting with a strong 


hydroxide, vield specific chemical compounds, as e.g.: 


Al + KOH + 3H,0 — K[AL (OH),] + 2 H, 


Zn + 2KOH + 2H,0 — K,[Zn (OH),] + H} 
Sn + 2KOH + 4H,0 > K,|Sn (OH),] + 2H,. 


NaOH reacts with the amphoteric metals similarly. Obviously, the foregoing compounds 
are by no means the sole ones that can be formed, as e.g. Al and Zn also yield compounds 
of the type of KAIO, or K,ZnO,. 

The following facts point to the chemical reaction between the amphoteric metals 
and the hydroxides NaOH and KOH as accounting for the photo-stimulated exo-electron 
emission observed: 

l. Emission occurs only when an amphoteric metal is present and in contact with 
NaOH or KOH. 

2. Such contact produces no emission in the oxides of these metals. 

3. As the reaction. between the hydroxide and metal proceeds to termination, the 
emission. vanishes. 

4. No such emission could be obtained for a number of reactions with acids (for the 
geometry of the present experiments). 

It is improbable that the pulses in the counter should be due to the hydrogen set free 
in the reaction, as is the case for various vapours (Sujak, in print), since hydrogen is liberated 
also without the intervention of white light, whereas the exo-electron emission observed 
requires that the specimen be thus illuminated. 

Hence, the authors consider that the emission of exo-electrons as observed in the pres- 
ent investigation is due to a photo-electric effect accompanying the chemical reaction 
between the amphoteric metals subjected to deformation and the hydroxides. The chemical 
affinity for hydroxides is known to increase from Pb to Al. The experimental shift of the 
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long-wave limit of the photo-effect towards shorter wavelengths occurs more or less in accord- 
ance with the decrease in their chemical affinity. Thus, the less active metals require light 
quanta of higher energy to produce exo-electron emission than the more active ones. 

Among the metals investigated, Sb also is characterized by low affinity for hydroxides. 
The authors consider that the reason why no emission is observed here resides in the low 
rate of the reaction or in the absence of quanta of sufficient energy from the light used for 
illumination, as was the case of Pb in the first stage of the present investigation. 

A more detailed investigation is proceeding. 
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PION INTERACTION IN FERMI STATISTICAL THEORY 
By V. S. BARASHENKOV 

Joint Institute for Nuclear Research, Dubna, Laboratory of Theoretical Physics 
(Received December 8, 1960) 


By assuming a 7 7- resonant interaction the energy spectra of pions and nucleons generated 
in inelastic NN collisions at E=9 BeV have been calculated. The results of the calculations are 
compared with experiment. 


The effect of a possible, rather frequent, resonant zr-interaction on calculations 
made according to the theory of multiple production of particles [1—8] has been recently 
discussed in the litterature. In Refs. [3, 6, 7] it has been shown that, for xN- and NN- 
interactions, taking such interaction into account within the framework of the Fermi theory 
leads to worse agreement between the calculated and experimental distributions of stars 
in multiplicity. On the other hand, in Refs. [2, 5] it became possible to obtain a better agree- 
ment between experiment and the theoretical momentum distributions of particles produced 
in æ p-collisions at H=1.0 BeV and E=1.4BeV. However, at energies E~1 BeV, 
when a very small number of secondary particles is produced, the statistical calculations 
should be made very carefully. Indeed, a more accurate analysis performed in [5] has shown 
that taking into account of the resonant gz-interaction fails to explain the experiment 
so well as was pointed out earlier in [2]; the agreement with experiment can be obtained 
without taking az-interaction into account if all the statistical calculations are made more 
accurately. 

Since the influence of the resonant interaction of pions becomes greater as the energy 
increases, we calculated the momentum spectra of secondary pions and protons for pp- 
collisions at E=9 BeV. The results of the calculations are given in Figs. 1 and 2. In these 
caleulations, the spin and isotopic spin of a pion isobar were assumed equal to unity, and 
its mass u* — the fourfold mass of pion: u*=4u. All other assumptions and the method 
of calculation are the same as in [3]. The experimental data given in Figs. 1 and 2 are taken 
from papers [9, 10]. The complete experimental data are presented for protons. 


1 [n Refs. [9,10] 28 pion tracks and 106 proton tracks have been analysed. In Figs. 1 and 2, the mean sta- 
tistical errors are indicated. At the same time, the total number of protons was considered to be 164, since in 


[10] 40 protons are added due to various corrections. 
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Fig. 1. Momentum distributions of pions in the c.m.s. The continuous curve shows the results calculated 
the results obtained with z-interaction. The dashed histo- 


without zrz-interaction; the dotted curve 
gram is that of the experimental data. The values of the momentum p are given in units of luc, where 4 is the 


mass of a pion. 


Fig. 2. The momentum distributions of protons in the c.m.s. The continuous curve shows the results calculated 
without zrr-interaction; the dashed line — these calculated with TUrdnteraction. The dashed histogram 
is that of the experimental data. The values of the momentum p are given in units of uc where u is the mass 


of a pion. 
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As we see, in both cases, the momentum distributions in the theory without mm - 
and with zz -interaction are found to be close to one another. Eo according to 
the results of papers [3, 6], taking into account of the resonant interaction of pions dita 
‘not on the whole improve the agreement between theory and experiment (Cf. [11]). This 
is also seen from the data of Table, wherein the theoretical and experimental values of the 
mean momenta of the particles generated are assembled. According to the theory which 
takes into account az - interaction, the average number of particles generated per act 
exceeds the experimental value by approximately 20%. 


TABLE 


Mean momenta of particles in the c.m.s. 


Kind of particles 
Variant of calculation | ya protons 

J 

Theory without (zz)-interac- 

tion p (BeV/c) 0.57 0.79 
| 

Theory with (a)-interaction | 

p (BeV/c) | 0.42 | 0.71 

Experiment? p (BeV/e) | 0.42- 0.1 | 1.1+0.3 


? The mean dispersion A p = Va EIN (Pos p N,, is given. 
nm 


On the other hand, in [1] it has been shown that without considering az -interactions 
it is difficult to account for the experimental data on annihilation processes. Subsequently, 
this conclusion was confirmed by the more accurate calculations of [4. 7. 8]. 

Thus, within the framework of the Fermi statistical theory, a consistent consideration 
of the resonant zz -interaction presents difficulties. 

The problem of resonant interaction between pions in the theory taking into account 
peripheral collisions (see, e.g. [11, 12]) is solved in a new way. The difference in the multipli- 
cities and the momentum spectra of particles generated in the central and peripheral colli- 
sions and values of the cross sections yet unknown allows to take into account sz - interac- 
tion without contradiction with experiment. However, final conclusions will require more 


detailed experimental data and numerical calculations. 
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ON THE INFLUENCE OF A MAGNETIC DIPOLE UPON THE GRAVITA- 
TIONAL FIELD 


By JarosLav PACHNER 
Physical Institute, Technical Universitety, Pracha CSSR 
(Received December 28, 1960) 


The influence of a magnetic dipole upon the gravitational field is investigated under the 
assumptions of a body with spherically symmetrical distribution of mass and with a weak magnetic 
moment. The resultant relations prove that the magnetic dipole exhibits also a weak gravitatio- 
nal quadrupole. 

The discussion of the results shows that the general relativity theory can play a certain role 
in the future quantum theory of elementary particles, and that there exist no measurable differences 
between the general relativity theory of Einstein-Maxwell and the unified field theory considered 
in this paper in the domains outside the elementary particles. 

A concluding remark deals with the influence of the magnetic fields in interstellar matter 


on the curvature of the universe. 


Introduction 


There is no doubt that the general relativity theory of pure gravitational field is the 
nly genuine field theory hitherto known. From the logical point of view there exists, how- 
"ver, a grave objection against its extension to the electromagnetic field, for this field enters 
nto it through the energy tensor which has been always considered by Einstein as a provi- 
ional means of representing matter. The four-dimensional asymmetric unified field theory 
ff gravitation and electricity proposed by the author is one of several attempts which intend 
o remove also the mentioned objection. 

The aim of the author in the present paper is not to discuss the physical and logical 
oundations of his theory, but to investigate the influence of a magnetic dipole upon the 
‚ravitational field and to find by this way the domain where the differences between his unified 
ield theory and the general relativity theory of Einstein-Maxwell take measurable values. 

In our investigation we shall assume a body with spherically symmetrical distribution 
f mass and with a magnetic dipole. Further we shall suppose that this body does not rotate 
n order not to complicate the calculation by unimportant details for our problem. 

In the first part of the paper the problem is treated exactly; in the second one the 
ield equations are simplified and then integrated under the assumption of a strong gravita- 
ional and a weak magnetic field. A third and last part contains a discussion of the results 


btained. 
(415) 
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Exact theory 
1. General form of the field equations 


The field equations of the four-dimensional asymmetric unified field theory of gravitation 


and electricity proposed in paper [1] are as follows: 


S v SA oz DY Be (1.1) 
uv t ue 
g ,—4n DS (1.2) 
k 
Fuse] His] 0, (1.3) 
Fun ren E (or Ls -— 0. (1.4) 


In empty space the field equations (1.1) and (1.2) become: 


Ry + Hw = 0 (1.1a) 
ur 
G3 =. (1.2a) 
and there exist then the following identities among the field equations - 
gr 

(Fu TOS I’ 0, (1.5) 
w 

O „= (1.6) 

[L(R + Ho) ger] =0 1.7 

12\ [so] [aro] 2 ].= y (4.1) 


These equations are to be supplemented by the definitions 


SU TE | a m Ya TB N 1B ta 
Iu ar Mi w,a | ua,” DI aß 1 p Br? (1.8) 
ap ap 
he age eee oe 
H,, = — 2d (g Supßay — 48 Gap) (1.9) 
ap ap 
—2 ^ = " p . 
Hu d [2 (Ear T SupSar) » 4 8 EabE uv A Eu): (1.10) 
~~ > „oc ^ FT 
Tu = (t, — 3 Re aftu) + (Hyp — 1 Hg, )]l/ — det g,, (1.11) 
2 
of? — — » ç 
s Lye log det gap, (1.12) 
Ges 24 ge V c det Sap: (1.13) 
veo à : 7 N, -— 5 k č . : = 
By £77" we denote here the well-known fully skew-symmetrical tensor density with the 


components +1, 0, and by d the new universal constant of the unified field theory under. 


consideration. All the physical quantities are measured in geometrical units. 
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The partial derivative is indicated by a comma, the pta derivative with respect 


to I by a semicolon, the covariant derivative with respect to { us] \ (the Riemannian derivativ e) 


by a colon, and the covariant derivative with respect to Z " by a triple-point (see Eq. (1.5)). 
The cyclical divergence of an arbitrary covariant skew-symmetrical tensor a, is written as 


v 


follows: 


— af | Hi ; 
Cur: o] = Uso] C uv, C vost FO any (1.14) 


The symmetrical and. skew-symmetrical parts of an arbitrary tensor a, 
E^ are denoted by aps Lii and a,,, I"! 


vu yv? uv ? ut 9 Ihe ? 


or affinity 


z 
respectively. The same notation holds also for 
the contravariant components of a tensor or of a tensor density. 

With the help of the definition relations (1.8)—(1.13) we can solve the system of the 
18 differential equations (1.1), (1.2). (1.3) and of the 64 algebraic equations (1.4) and de- 
termine the 16 unknown components of the fundamental tensor g,, and the 64 unknown 
components of the affinity Jp. The 6 identities (1.5), (1.6), (1.7) reduce the number of 
the independent differential fold equations from 18 to 12 as requested by the postulate 
of general covariance of the field equations. 


2. The rotationally symmetrical fundamental tensor 


In our investigation we shall use troughout the system of spherical coordinates (xj. xs, 
E Xa) =r, 2, p, t). 

Since the symmetrical and the skew-symmetrical parts of the fundamental tensor 
transform independently from each other, we determine the rotationally symmetrical form 
of the fundamental tensor for each of its parts separately. 

By a suitable choice of coordinates the symmetrical part of the fundamental tensor 


becomes diagonal and does not depend on the angle g [2, 3]. We write then 


du D 0 0 
ae 0 er? 0 0 (2) 
= 0 0. eran? 0 
0 0 0 aoe 


Since we have assumed a body with spherically symmetrical distribution of mass, the 
rotational symmetry of the field is caused by the magnetic dipole only. Its influence on 
the symmetrical part of the fundamental tensor is expressed by the functions N, 0, 0, c. 


[f the body has no magnetic dipole, then 


n=0=o=0=0. (22) 
The explicit form of the functions y, v is known from Schwarzschild’s exact solution 
Xf the spherically symmetrical gravitational field: 


ef =e = l-2mir, (2.3) 
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m being the mass of the body. From Eq. (2.3) we get 
utr—=l. (2.4) 


. . , , , , J res. 
In the "cartesian"" coordinates (x3. xa, x4. x4.) = (x, y, Z, 7) defined by the condition 


that at infinity we have the pseudoeuclidean metric (—1, —1, —1, +1), all the skew-symmet- 
rical components of the rotationally symmetrical fundamental tensor differ in general 
from zero. After transforming to spherical coordinates, two of the components vanish and 
the remaining four do not depend on the angle g if the axis of the rotation symmetry coin- 
cides with the z-axis. 

The transformation formulae are 


x = r sin Ò cosp, z=rcos®, 


y-—rsinósng, t=t; 


With their help, we find 


, , . 
812 = (— 4&5 COS P — Bgg sin p)r, 


, . , . , \ . 
813 = [812 sin Ò + (e35 sin p — 83, cos g) cos Yr sin 9, 


, ; , . , . © . 
823 = [812 cos Ü — (gi, sin p — £54 cos g) sin Or? sin J, 


, , . . 
gia = [(g4 COS p + ga sin g) sin d+ £g4 cos OI, 


, , . , . 
Soa = (Sia COS P+ gg, sin p) cos Ò — gay sin Ölr, 


v 


, . / . 
831 = (— gy Sn P+ 234 cos gr sin 9. 


v 


It follows hence that the components 812 and g,, vanish if the body is rotationally symmetrical 


about the z-axis. We may then write 


0 0 ur sin d w 
0 0 vr? sin 9 ST 

imn ur sin Ù — vr? sin & 0 quet n 
— 10 — sr 0 0 


Ihe functions u, v, w, s depend merely on the coordinates (r, 0). In our case, when we 
consider only the magnetic components of the field, we have 


t0 -- sm). (2.6) 
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Ihe covariant components of our rotationally symmetrie fundamental tensor are thus 


t Ly ^ 5 n =. 
— ct 0 ur sin d 0 
0 — eer? vr? sin J 0 
o = 
eu * Dt gio org 
— ur sn 9 — vr? sin? — e?r? sin? 0 0 
0 0 0 dena 


In order to determine its contravariant components, we compute first 


— det g,5 = D exp (uty - n-- o - o —- e) r* sin? ð, 


where 


2,—0—96G 


D=1+u2e "TE ve 


With the help of formula (1.12), we obtain 


pile eiD, 
g2? = — ect (L4 ue") | pg, 
g? = e ?[Dr? sin? 9, 

git era 


m gene ED cams M i J . 
E mue IDrsan o 


23 
a = oe € Dr sin e. 


)ther components vanish. 


The contravariant components of the tensor density given by (1.13) are 


quo exp [uty ++ oo 0) 7? sin 9 VD. 


3. Solution of the system of the algebraic equations (1.4) 


(2.10) 


The solution of this system of the 64 algebraic linear equations for the unknowns 777, 
vas carried out by Hlavaty [4] and by Mme Tonnelat [5]. Since the formulae found by Hlavaty 
vith the help of tensor calculus are rather complicated, we use here those given by Mme 


"onnelat, which are more suitable for our purpose. 
The unknowns J", are determined by the relation 


uv 


€ € € € 


(3.1) 


The Christoffel symbols are to be computed either by the well-known. expression 


il 
€ € ae 
i $ 2 h : (uns LP ha)» 


(3.2) 
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or, more advantageously, by means of Dingle’s formulae [3]. In Eq. (3.2) 


^ 


e 
hu = Sw m Er log det hes. (3.3) 


uv 


The symmetrical tensor VE, is given by the equation 


Ww 
AT é 0 to o © 
NS = he gi no or - ro au)" (3.4) 
Since in our case 
det Ey = 0, (3.5) 
5 


the skew-symmetrical tensor /77, is given, according to Mme Tonnelat, by the relation 


uv 


ema otin ¢ 
Dy =h [A ms Soo Er jh h Aum (2 i D)]. (3.6) 
where 
/ — T | P MM aw Br ‚orin 
A o Dur Styw,0] T Suro T 5 Saah h Bast ^ SAnS[or,o] 
D v v v v vi 


px (Sv). ==) By T we NE )/2D tC uro (3.7) 


The tensor C o \ anishes if (1.2a) is satisfied. Beeause we shall always suppose the fulfilment 
of this condition, we do not quote here the explicit expression for Cures €45,54 and geP? qe. 
vv 


note here the fully skew-symmetrical indicators of Levi-Civita with the components +1, 0. 

In our case of the rotationally symmetric field, where the fundamental tensor has the 
form given by Eq. (2.7), the determinant of ha; and all the contravariant components A" 
which are different from zero are 


det 4,5 = exp(u+ + y+ 0+6+4 o)rsin? 0; (3.8) 
1 Sen 29 _ _ ot]! 
hil = —¢ S h22 — — e-elp, 
A93 = — e-"|r* sin? Q, h= e-*-". (3.9) 


Ihe Christoffel symbols that do not vanish are explicitely given by the formulae 


fi) 1 [u ên (1) e 
In A \ore ears? 22] = — |r+ 5 "3 e-"-nte, 
J | | l co jl 1 Qv o 

la] "ADEL v m Te: (2 = or oa 
[2| - ERAS T le 

lf 389 8i l| 7 3 85 

J2 | 12 y Mss 

Bie sin d cos Ô + 339 sinio ee P =; 2 chen 
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V ta e e LL 

ux 29° Bl ner 

Den {ay = cote 94+ 2 22 

tit - l1 [ovr _ [4 l ew 

14 2 \dr or }’ put o5 (5.10) 


The computation of the components of the tensor ue by means of formulae (3.6) 


and (3.7) is rather tedious. The non-vanishing components are 


r= =h" Ust &13:1 — 813D3/D) [1 + 8358,55 1 A99/(2 — D)] + 


+ [3 (823,1 — 813,2) + 813: 2 (S13) 2 + 823D 1)/2D] eygers 597. 


2 96 : 
In = h?2 113 (823,1 813,2) + 813:2 — (213) 2 + £3) 1)] |l 


+ D)] + (€13.4 813 D,/D) et — Di 


TES = | 1 (22 1 813,2) 823.1 (S23) 1 = £13) 9) ][1 — Suse uos a ei D)| ae: 


xis (823:2 — S93 DD) esos 99, 
zo 32 : 
J os — 422 [23.2 = S23) |D) [1 se Besos)” 2h33] (2 Er D)| 4 


(S23) 1 al S13) 2)/2D| 823813 h (2 = D)}, 


Ta T M i (e 823,1 813,2) um (S23) = S13) »)} [ier 
+ (8358151 1h? + Sosogh??h??)/(2 — D)]. (3.11) 


Since we know the components 777, that do not vanish, we easily determine the compo- 


nents of the tensor N, which are different from zero: 


Ni =— 2 inea, Ne = — 2h (I 8154 I Neo). 
Na, = -2[ igo, N = = DA 23813 T Nas): 
NA = Tha NS — (Th gya + Ree) 


Nig = — Tigh”, NE ER. 23813 + Tg): (3.12) 
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Let us quote here the following two formulae (where the usual summation rule holds) 

which will be useful for our further calculations: | 


Ds = (log V— det z,,),, : (3.13) 
P — 0, if Eq. (1.2a) is satisfied. (3.14) 


Substituting Eq. (2.8) into (3.13) we have 


pe 2. 108 j 1 aD 
SSS 2c Vis Qut 
M TA a ihe 2D Ər 
1 28 1 OD 
T0. NL. HS d — Un wt 1 E 
Ta = cotg J 2 95 (2j 2- 0 —- 8 4- 0) 3D 35 (3.15) 


The P with ~=3, 4 vanish. 


4. The differential field equations in the rotationally symmetric case 


In this case the fundamental tensor has the form given by Eq. (2.7). From the 10 diffe- 
rential field equations (1.1a) only the following 5 do not vanish identically. They have been 
computed with the help of Eqs. (1.8) and (1.9). We have found: 
ior up = 11: 

b= Daan Ca e DRUSI DD TE r 


— y 1 
le,1 ILA 


2 (€ 1 7 D 13 2 —1 22 | 
In (215 —1 2a) — 2D] = 0, (81345511 499 — Besos ?* 33) = 0, (4.1) 


~ 


for uv = 22: 


oggi S 1 9 12 — gl 3 "3 r13 2 ta E 2 2 12 "€ *1 7'1 
[ D$ 1 I3 (23 I) l al I PE Tsa, + Toa US — Us) Tal T 


4 ‘ ' SA 29132: ^ 
p Pal +2] Tal 23] — d? h,4 (— ggg h1A33 + EosGogh?*h33) = 0, (4.2) 
for uy = 33: 
= Lr 11 (9713 Se ey rm ‚2 jc 2 ’ 1711 972 
[ uv kJ 33 es T'ia) My, Hl s3 CDs a IS) Yd Lys! 1370 2D 3 io 


rey 2 EH y 9€ 
— Lest os] — d^? has (13813h h33 + 8238234 h33) = 0, (4.3) 


for uv = 44: 


[- Tant TAQTA - T3) - T TAGPA — TE] 


-9 E r 
— dhg (= B13g13h!h33 — Egg og 22/33) = 0, (4.4) 


and for uv —12: 


ils Tal That DADA Dh-D8)— DÀ 4 PAT 8 + qu pae 


a 


HTAA Peet LRU TE — 24-389... —0. zu 
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From the 4 identities (1.5), we see that among the preceding equations only 2 (i.e. for e=1, 2) 
are non-trivial. They reduce the number of independent field equations (4.1)— (4.5) from 
to 3. 
From the 4 field equations (1.2a) only the following one (for u=3) does not vanish 
identically 
31 233 
$3 FJa — 0. (4.6) 
The identity (1.6) is here trivial. 
From the 4 field equations (1.3) only the following one remains as non-trivial in our case 


5 
\ ae. e Ji 11 77:2 Ub 3 79 7 732 
(Rosa + Fio) ar (Ha Hg 2) = = D 23 EI 12 t Ds T3 DA. ! TUI 
2 
SENE pic TIS prier: temas rd = Va Tl ge eT us 1 ae 
F T33 Do) Lial s3 i T'i p L'il 22] zg l isa Hl is 11 ur I 1 13,27 


12 t | ı3 AE E 1 12 | Dap TP Re ¢ ea, 
i3 Dig | Dos I3) HJ id 33 | 1 25] Ta 1 sal a +2d %(Sozı — 813 2) == (4.7) 


v 


In the considered case the identity (1.7) is also trivial. 


Maxwell’s approximation 


5. Definition of the Maxwell approximation and its relation to the general relativity 
theory of Einstein-Maxwell 


The field equations in their exact form are too complicated for an exact solution. The 

most natural approximation is that in which the symmetric (gravitational) part of the unified. 
field is assumed to be strong, and the skew-symmetrical (electromagnetic) part to be weak. 
In this approximation, which was called the Maxwell approximation [1, 6], all the squares 
and higher powers of g,, are neglected with the only exception of the terms where the 


square of g,,is divided by the square of the very small universal constant d. The field equations 


keep then their form, but the definition equations change. We have now 


en let on 
Run = — i hP Bruna]: at LO gto: 8» (5.2) 

H,, = 2d (— g, 48, hi + T Eagle h hur): (5.3) 

Hp = 24 78, (5.4) 


Fn = Ps — 3 Pggh? Pn) + 2d *(— gs, 9-3 arto h 7} V— det h,,, (5.5) 


um Erat oe $ 
g = AETIA (5.6) 
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The solution of the 64 algebraic field equations (1.4) is very simple [6]: 


Í l 
) € ' 
jl pe = ja J " = h” (- 3 ftl m seu]. (3.7) 
Multiplying Eq. (1.3) by 127, we obtain in the Maxwell approximation 
1 72 TAF A 
2 d Rune) O[uw,0] -— 0. (5.8) 


This relation differs from the second set of the Maxwell equations only by the first term 
which may be neglected in many cases because of the smallness of the universal constant d. 
This is the only term in which the field equations in the Maxwell approximation are different 
from those of the general relativity theory of Einstein-Maxwell. 

In our case of a rotationally symmetric field, where the fundamental tensor is given by 
Eq. (2.7), we consider the functions u and v as the strong ones and the functions u, v, 1j. 0,0, @ 
as so weak that their squares and higher powers may be neglected if they do not appear 
in Hy ; H, , multiplied by the very large constant d~*. It is quite natural to apply this rule 
through all the further calculations, where the magnetic field is considered as a small pertur- 
bation. The integration of the differential field equations becomes then possible, but the 


results will be only approximative even within the scope of the general relativity theory. 


6. The solution of Eq. (4.6) 


We begin our computation with Eq. (4.6). In the Maxwell aproximation this equation 
takes the simple form 


m 
- (ure—*) ag 0, 6.1 
eo d 


JE 


thanks lo Eqs. 2 10), 2 11) and (2: 4). 
With regard to the equations of the Maxwell the ory describing in spherical coordinates 


the distribution of the magnetie field produced by a magnetic dipole, we assume 
u(r, 0) = U(r) sin 9, — v(r, 0) = F (r) cos 0. (6.2) 


After having substituted them into Eq. (6.1) we obtain the relation between the unknown 
functions U(r), V(r): 


2 1U 
Vm = (Tre?) ee (os Onl ar: 
(Ure) = (r — 2m) 3i U (6.3) 


The latter expression follows by applying Eq. (2,3): 


The ? T2 »Ynpre ] 6 
The explicit expressions for DS 


As it follows from Eqs. (5.7) and (3.10) the explicit expressions for /'5, are given by 
the formulae ; 


l [ou , on 
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I3 cer [enm 2 e re^" sin? 9, rh = E (2 ento + = entr. 
i=- ga og TÀ- 35 

I$ = —e-et9 sin 3 cos à — 95 sin? 9, A = E 25 e”, 

TA d DÀ. cotg 0+ 5. 

nana (5+ 3) Ph sh om 


The non-vanishing components of /7, computed with the help of Eqs. (6.2), (6.3) are 


uv 
v 


T =— sn?ó L 2m) qu Lx | , 
Y dr r 


— — sin 9 cos 9 : p gu 
—— s 0: olm d s 


277 , 
1% =— sin À cos à G — Amr? + 4m?r) > + (r? — 4m?) | ; 


T: 


er SD) 


T € 
T = sin? L 2m) a 2m | 
y dr "y 


l d?U 2m\ dU 

a e » 

Ti = cotg J l 97 m) IR c (2 > ) 2d : (7.2) 
We give here also the Maxwell approximation of Eqs. (3.15): 
roe PNE ES. 
Tja = Sa 39, Ht 15-0 + o-+ m), 

a Id " 

Dog = coig O t 5-55 (gc et et a). (7.3) 


8, The solution of Eqs. (4.7) 
In the Maxwell approximation, Eq. (4.7) takes due to (6.2), (6.3) the following form 


2m due. m\ dU AB 2m |" d*U 
eim ELM. (1 En d e. 1— = ga 


r 


pafa 10e, ts mero 


1 Ta | 
r r2 | dr? r re j dr? 
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3: 2m  2m?*| dU Ly #7? i (8.1) 
UNDE 4 208° r? dr ro mar | 


It is an ordinary linear differential equation for the function Ufr). We shall integrate it by 
successive approximations. 
In the first approximation d —0. Eq. (8.1) is then a linear differential equation of the 


first order for dU/dr. By repeated. integration we find 


fita = dr fe Ky | —2mr+ 2m? ls ( " =] +K, 
em ! I r(r= 2m)? Am®| r(r — 2m) ; r 


ES | jb 2 ES (8.2) 


2 
E r^ 


K,, K, being integration constants. We now put the first term of the approximate solution 


for U} into the right-hand side of Eq. (8.1). and get 


2m d2U, ; m at 11K,d? f , 
3 E 2 |. A72 a —(. 8.3 
j | : r ) dr? a | : ridr 2r4 ^ 0 


Integrating this linear differential equation of the first order for d, dr. we obtain 


dU, SAP ‚4m | E lld? 


— E i = * y 
dr r* ^: 75 Ar® 


Integrating again, we get the second approximation 


Ki 1% 3m a 144m? | 4 K.. (8.4) 
Pa 20r? E 


We now compute by Eq. (6.3) the explicit form of the function V(r) using the function 


U(r) in the second approximation 


(8.5) 


For m=d=0 Eqs. (8.4) and (8.5) describe the distribution of the magnetic feld in 
Maxwell’s theory. The mass of the body is not the source of the magnetic field in the general 
relativity theory, but it does influence it by the terms with m0. The terms with d40 
are due to the existence of the first term in Eq. (5.8), by which the field equations of the 
unified field theory under consideration differ in the Maxwell approximation from the field 
equations of the general relativity theory. 

Denoting the magnetic moment of the dipole expressed in the geometrical units by Sand 
being aware of the fact that the magnetic field vanishes at infinity, the integration constants 
Ki. A, have the values! 

K,=—3Sd, K,=0. (8.6) 


) See Eq. (69) in paper [1] expressing the relation of the skew-symmetric components of the fundamental 
ensor to the intensities of the electromagnetic field. 
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In the further computation we shall substitute for U(r), V(r) the functions 
U(r) = Sd|r*, | V(r) = — 26d[r*, (8.7) 


neglecting all the correction terms, as their influence is negligibly small outside the regions 
occupied by the elementary particles. 


9. The solution of Eqs. (4.1)—(4.5) 


On transcribing the system of the 5 non-linear partial differential equations (4.1)—(4.5) 
E. the 4 unknown functions n, 0, o, œ into the Maxwell aproximation, we take into account 

s. (2.3). Substituting then the linear terms from the power series for the expone ntial 
E and using Eqs. (6.2), we obtain 


2m \ 2? E M ee : EM ges de 
2e + (Rut w= (1 5 IF 5 (0 FG o)d 3 m 2n + 20 + 20) — 


r r er 
mo, l[9?n Mm 
ee Ges 0) d n otg 9 | + 
2 5, | 37 + 30 + 30 — 3w) E: [en 20 cotg ù | 
+ 2d-? |o (1 — a sin? V — V2 cos? o| - 0, (9.1) 
r 
2m\2*o 10 as MB m 2 | 
(2/r?) (Rog + Ha) = (1 7 | Ir? : 3, | H+ 90 3- 0 3- @) 2 Or Zn 4 
€ € l et | f | e LI G E , 
+ 4o + 20 + 2) 4 | op (9 4- 6 4- eJ 59 | 0 + 20) cotg 9 — 2n+ 20 | + 
+ 2d-? |- U? (1 — a sin? 9 -- V? cos? 2 = 0, (9.2) 
r 
om on A 
(2/r? sin? 9) (R3, + Ha) = (1 zr = + = (-—n+0o+30+0) 
: ar Ton 
e 25 + 20 + 40 + 20) + Sha E o +20 + w) cotg à — 29 + 20 | + 
— — 7] ) I r2 Pond ^ $9 ] Uo ah 
+ 2d? |^ n = 25) sin? 9 + V? cos? J =(), (9.3) 
F 
DaN ecm a m 2 
22607 (Ra t Hs (1 > ) 572 + : 3; Qo) E: z, 0—20 4- 0) 4 


1/020 | 9o | | 4 2 72 ns? : zs I 
E 20 os a —U2 (1 — =) sin? 9? cos? 9] =0, (94) 
s Te VAT d so) Y ; 


2m 


2m\ 9? 9 en 
2e”(Ri + H) = (1 : Js w) + (1 a Jic - o) cotg 0 + 


12, 
JOE UM 


He + a + Ad-? UV (r — 2m) sin 0 cos 0 — 0. (9.5) 


e) + 


nd Aen 
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Putting here for U(r), V(r) the functions (8.7) and at the same time neglecting, of course, 


all the terms with m, we may suppose the unknown functions to have the form 

n(r, 0) = (S?/r) (A, sin? 9 4- B, cos? 9), 

olr, 9) = (S?[r*) (A, sin? O + B, cos? 9), 

a(r, à) = (5?/r*) (A, sin? + B, cos? 0), 

o(r, 9) = (S?/r*) (A, sin? 0+ B, cos? 0). (9.6) 
By this way the system of the 5 linear partial differential equations (9.1) —(9.5) reduces 


to the following system of 10 linear algebraic equations for the 8 unknowns ER 


B, (v—1, 2, 3, 4) 


SADR FEA PEAL SNO (9.7) 
A, E B, 4-8D,4- 3B, + 5B, 42, (9.8) 
Boda P Boa Ae ea (9.9) 
A, — A, + 6B, + 34, — 5B, -- A, — 3B, — — 4, (9.10) 
Ay A Achaea (9.11) 
A, Ag+ 34; + BEA = ME (9.12) 
SAD (9.13) 
AIR, (9.14) 
4, — B, — 24,4- 64, — 4B,-- 54, — 5B, = — 4, (9.15) 
p By ne (9.16) 


Since there exist 2 non-trivial identities among Eqs. (4.1)—(4.5) as well as among Eqs. 
(9.1)—(9.5), there are merely 6 independent equations among (9.7)—(9.16). Denoting by 
a, P two arbitrary real numbers, the unknown coefficients in Eqs. (9.6) are given by the 

à I S ) 
formulae obtained by solving the system of Eqs. (9.7)—(9.16): 


=g B,=a+? $ 

l =—}a— 4p B,=— 4—3% +}, 

a=—-4+—1a B-—i—ia—i, 

"RENS dr (9.17) 


If we introduce into the functions (9.6) the Legendre polynomial of the first kind 


P. (9) = M cos? 9 — i d 
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these functions take the form 
mr, 9) = (S*/r*) | (« 4 i b)+ PP;(9)], 
et D) = (CA —1a—15) 1,09]. 
o(r, 8) = (St [=a 4p) + (—2 te], 


or, vy = (S*/r*) ie + $ P,(0)]. (9.18) 
Choosing with advantage 
a=0, p=, (9.19) 
we obtain 
n(r, 8) = or, 8) = + (S*/r*) [1+ 2P,(9)], 
o(r, 9) = a(r, 0) = — 3 (S*/r*) [1+ P(9)]. (9.20) 


Substituting at last these functions and relations (2.3) into the symmetric components 
of the fundamental tensor, i.e., into the gravitational potentials given by Eq. (2.1), we find 


with the same degree of approximation as used hitherto 


2m \~* ; 
(1-22) ped enara 


r 


M — N + (2m/r) + E (637) ay. 


mr Fy + A]. 


1 
£33 = sin?) h == (S?/r4) [L+ Do : 


| fas — 3 (1 er is ro! 
~ + h- (2m/r) + z G3 [1+ 2b. (9.21) 


The resultant relations (9.21) show that the magnetic dipole exhibits also a gravita- 


tional quadrupole, even in the case when the mass of the body should tend to zero. 


Discussion 
10. On the quantitative relation of the unified field theory to the general relativity theory 


The qualitative relations of the unified field theory under consideration to the general 
relativity theory have been already mentioned in Sec. 5. In order to be able to discuss the 
quantitative aspect of the problem we determine first the critical radius rom at which the 
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influence of the magnetic dipole on the gravitational potentials takes approximately the 


sune value as that of the mass. From Eqs. (9.21). we get 


rom = VS2/6m. (10. 1a) 


The numerical estimation will be made easier by replacing formula (10.1a) written in the 
geometrical units by the following one, where S* is measured in (gauss x cm) and m* in 


grams [1]. 
3 Ye 9/7 k 92 
Tom = VS*2/6m*e?, : (10.1b) 


c being the velocity of light 3 x 10! cm/sec. : 
Putting into Eq. (10.1b) the mass of the Sun m*~=2 x 108% grams and its magnetic 


moment S* x 8.4 x 1033 (gauss x cm), we find 


. 


Tom = 9x 10? cm. (10.2) 
The critical radius lies here below the radius of Schwarzschild’s singularity 


r, 2m zx 3x 10° cm. (10.3) 


In celestial mechanics the influence of the magnetic field on the gravitational effects is thus 
wholly negligible. 
Putting into Eq. (10.1b) the mass of the neutron m*~ 1.7 x107% grams and its magnetic 


moment 5*2 0.95 x 10 (gauss x em), we obtain 
Tom = ax l0ri5em. (10.4) 
In this case the critical radius highly surpasses the radius of the Schwarzschild's singularity 


n=2m23Xx ]079* om. (10.5) 


5 


Therefore, in the future theory of elementary. particles the gravitational effects caused. by 
the magnetic moments of the particles may be neglected only outside the regions occupied 
by these particles, 

Till now we have restricted ourselves to general relativity theory. The differences 
between the unified field theory in the Maxwell approximation and the general relativity 
theory manifest, themselves in the distribution of the magnetic field by the term containing 
the new universal constant d. Neglecting the mass of the body (which is allowed with a very 
good approximation in the case of elementary particles), we have from Eqs. (8.4). (8.5) 
and (8.6): 


U(r) = (Sd/r?) [1— (33/20) (d/r)?}, 

V(r) = — Gsdir) [1— (83/10) (d]n]. (10.6) 
According to the hypothesis formulated in paper [1] the universal constant d lies in the 
interval 


10-4 em < d —10-13 em. | (10.7) 
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If the value of d approaches its upper limit, the distribution of the magnetic field. inside 
the elementary particles determined by the unified field theory differs from the one given 
by the general relativity theory. In these domains we must not use, however, the Maxwell 
approximation, but we must carry out the computation by means of the exact theory. (Let 
us remark that the equation (8.1) has been integrated under the assumption d<r). 

In the case of an electrically charged mass-point, the field equations of the unified 
field theory were integrated exactly [7, 1]. The electric charge influences here the following 
two of the four non-vanishing gravitational potentials (the cosmological constant A being 
excluded for the reasons given by Einstein [8] and recently by the author [9]) . 


£u = — [1 — (2m/r) + (g]r)*] 3, 
gs = [1 + (a[r)*(d]r)*][1 — (2m]r) + (q/r)?| (10.8) 


The gravitational effects of the electric charge are comparable (within the scope of the 


general relativity theory) with those of the mass in the regions below the critical radius 
Tog = q?|2m. (10.9) 


For the electron with mz:9 x 1073 grams ~6.7 x 109 cm and q24.8x 10710 gauss 


~1.4 x 10-34 cm we have 


lle 


For l5 10713 em. (10.10) 


For the proton with mz 1.7 x 10~*4 grams ~ 1.3 x 1079? em and with the same electric 


charge, we get 


lop poc Wem. (10.11) 


Also these results show the importance of the general relativity theory for the future theory 
of elementary particles. À 

In the distribution of the electric field due to an electrically charged mass-point, there 
is no difference between our unified field theory and general relativity theory. This difference 
appears only in the gravitational potential 24, (second term in the first bracket of Eq. (10.8)) 


and manifests itself in the regions 


r, < V qd. (10.12) 
Putting into this formula the elementary electric charge and the upper limit of d, we have 
rp 0229 em. (10.13) 


From this numerical discussion we may conclude: 

1. The general relativity theory may play a certain role in the future quantum theory 
Xf elementary particles. 

2. There are no measurable differences between the general relativity theory of Einstein- 
Maxwell and the unified field theory under consideration in the domains outside the ele- 
nentary particles. The latter is justifiable in empty space merely from the logical point 


of view [10]. 
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11. On the cosmological problem 


In the universe there are extensive regions of magnetic fields caused by interstellar 
matter. The question arises whether these fields influence the curvature of the universe 
and how this problem should be treated mathematically, for the astronomical observations 
testify the universe to be isotropic, while the magnetic fields show a characteristic axial 
symmetry. 

The answer is quite simple. We have to apply to the universe containing such polar 
fields the methods used in the theory of paramagnetic matter. The axes of these “elementary 
magnetic fields" are statistically so distributed that the universe as the whole is isotropic. 
This condition of isotropy can be satisfied if the radial components of the magnetic intensity 
vanish at the boundary of each "elementary magnetic cell’. 

If we describe in the first approximation these magnetic fields as caused by magnetic 
dipoles, the integration constant A, in Eqs. (8.4) and (8.5) does not vanish and is determined 
by the condition 


Ba. — 0, ie Plr)=0, (11.1) 
ry being the radius of the magnetic cell. The gravitational potentials are given also in this 
case by the functions in Eqs. (9.6) with the only difference that the coefficients 4,, B, (y= 
=], 2, 3, 4) are no constants more, but functions of the form 

/ $ " ; ^ 
A, =a, +a, (rir)? + a, (r/r,)®, 
, \ g 7 N > 
B, =b, + b, (rjr) + b, (r/ro)®, (11.2 


, 7 7 : = N : S i 
Q,,Q,,a,, and b b b, being real constants of the order of magnitude of the units. 


a diy doe 

In his recent paper [9] the author assumed the universe to be divided into N equal 
cells which are so large that the gravitational field, and due to it also the curvature of the 
space, may be considered with a sufficient accuracy as spherically symmetric at their bounda- 
ries. If we wish to take into account the magnetic fields of interstellar matter, we must 
suppose that each of the N cells contains a large number of “elementary magnetic cells’ 
which contribute to the curvature of the space. It seems, however, that only the rest mass I 


: ART : : 
and the *Hubble's mass* have a decisive influence on the curvature of the universe, 
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Fase. 


EINWIRKUNG VON FLUSSIGKEITSDAMPFEN AUF DEN LUFT- 
SPITZENZAHLER I 


* 


Von B. SuJAK 


Institut für Experimentalphysik der Universität, Wroclaw, Institut der Physik, Polnische Akademie der Wissen- 


schaften, Wroclaw 
( Eingegangen am 31 Januar 1961) 


Die Einwirkung von Flüssigkeitsdimpfen auf den Luftspitzenzühler wurde untersucht. 

Die verschiedene Flüssigkeiten, die unter das begitterte Ende eines Luftspitzenzählers 
unterstellt wurden, bewirkten das Auslösen von Arbeitsimpulsen im Zähler deren Geschwindigkeit 
und Gestallt zeigten sich von der Zeit der Einwirkung sowie von der Natur der Flüssigkeiten 
abhängig zu sein. 

Die sogenannten Zeit-Kurven, die die Zeitabhüngigkeit der Impulsengeschwindigkeit für 
verschiedene Flüssigkeitsdämpfe darstellen wurden untersucht. Auch die Spanungsabhüngigkeit 
dieser Zeit-Kurven sowie die Impulsgeschwindigkeit-Abhängigkeit von den Arbeitsspannung 


des Zählers wurde einführend gemessen. 


Einführung 


In einigen bisherigen Berichten habe mich mit dem Einfluss von Wasserdampf auf 
den Luftspitzenzühler befasst (Sujak 1958a, b, c; 1959; Sujak und Bohun 1959; Lewowski, 
Sujak und Wojciechowski—im Druck). 

Der vorliegende Bericht befasst sich mit der Einwirkung von anderen Flüssigkeitsdämp- 
fen auf den Luftspitzenzähler, der mit einer scharfen Spitze bestückt ist. 

Dieser Bericht soll auch weitere Erkenntnisse zur Frage der Anwendungsmöglichkeiten 
ines Luftspitzenzählers bei physikalischen Untersuchungen liefern, wobei keine radio- 
ıktive Indikatoren zu benutzen wären (Sujak 1957a, b, 1959c; Lewowski und Sujak 1959; 
Mader und Sujak 1960). 

Untersuchungen von dieser Art (Hrbek 1959) sind auch von grosser Bedeutung für 
lie Exoelektronen-Detektionstechnique Fragen, besonders in dem Falle wo die Einflüsse 
ler Flüssigkeiten oder deren Dämpfe auf die Untersuchten Exoelektronenemissions- 
örscheinungen, der Gegenstand der Studien sind (Lewowski und Sujak 1961; Lewo- 
vski 1960, 1961; Wawrzyniak und Sujak 1960; Seidel und Roubinek 1953, 1954) oder für 
lie Fälle wo Exoelektronenemissions-Erscheinungen mit einer Entgasung oder auch ka- 
alitischen Wirkung der Proben zusammenhängen (Bathow und Gobrecht 1956). 
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Impulsgeschwindigkeit als Funktion der Einwirkungszeit der Flüssigkeitsdämpfe 


Man beobachtet im allgemeinen eine zeitweilige Änderung der Impulsgeschwindigkeit 
bei der sonst konstant gehaltener Arbeitsspannung des Luftspitzenzählers, wenn die Luftat- 
mosphäre in welcher der Zähler betrieben ist, durch einen Flüssigkeitsdampf zeitweise 
beeinflüsst wird!. 

Um das Verhalten des Luftspitzenzählers gegenüber verschiedener Dämpfe überhaupt 
vergleichen zu können, muss ein Vorgehen festgelegt werden, das praktische Aussichten 
für weitere Untersuchungen verspricht. Die „Zeit-Kurve“ scheint ein solches Vorgehen 
festzulegen. Durch die Zeit-Kurve wird eine Abhängigkeit der Zählgeschwindigkeit Nft 
von der Zeit t, während welcher die Flüssigkeitsdämpfe auf den Versuchs-Luftspitzenzähler 
einwirken, gedacht. Die Zeit-Kurve kann unter Umständen z.B. zu der Zeitabhängigkeits- 
Kurve eines radioaktiven Zerfalles, der mit Hilfe eines Geigerzählers gemessen wäre, vergli- 
chen werden. 

Die Zeit-Kurve kann wie folgt bezeichnet werden: 

| 
X Sf D 

Wenn ein Behälter mit der zur Untersuchung gelangter Flüssigkeit unter das begitterte 
Ende des Luftspitzenzählers unterschoben wird, so werden für verschiedene Flüssigkeiten 
im allgemeinen verschiedene Verläufe der Zeit-Kurve (1) beobachtet. 

Der zeitliche Verlauf der Kurve (1) muss offensichtlich wenigstens von den folgenden 
Parameter abhängen: 

a) Geometrie der Versuchsanordnung (G), 

b) Temperatur (T), 

c) Arbeitsspannung des Luftspitzenzühlers (V), 

d) Eigenschaften der Flüssigkeiten (E), 

e) Flüssigkeitsmenge die zur Verfügung steht (Molekülenzahl) (M). 
und folgende symbolische Abhängigkeit kann dieses abspiegeln 

N 4 
E = f(t, GC, T, V, E, M) (2) 


Die ganz allgemeine Beziehung kann experimentell untersucht werden indem alle 
Parameter bis auf einen festgelegt werden. Wenn z.B. die Geometrie der Messanordnung 
festgelegt wird, so kann (G) als eine Konstante der Apparatur angesehen werden. Dasselbe 
gilt für (M) wenn entsprechende Menge der Flüssigkeiten benutzt wird, die ausreicht um 
z.B. in dem Messvolumen einen gesätigten Dampf für längere Zeit aufrecht zu erhalten. 
Wenn weiterhin die Temperatur (T) in welcher die Experimente durchgeführt werden, 
sowie die Arbeitsspannung (V) des Zählers für gegebene Messreihen nicht geändert bleibt, 
so kann eine einfache Form der Beziehung (2) aufgestellt werden: 


N j 
t]cnvEgM; ^0: (3) 


1 Dieses würde z.B. von dem Kernphysiker als „Untergrundänderung“ des Zählers bezeichnet. 
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wobei die Funktion f; verschieden verlaufen kann für verschiedene Flüssigkeiten mit den 
verschiedenen Eigenschaften E; 

Diese Form der Beziehung (3) scheint eine der wichtigsten Abhängigkeiten zu sein, 
besonders bei den einführenden Untersuchungen, die zeigen sollen, ob eine versprechende 
Aussicht besteht den Luftspitzenzähler als elektrische Nase zu verwenden. Aus diesem 
Grunde wird dieser Abhängigkeit die grösste Aufmerksamkeit in dem vorgelegten Bericht 
zeschenkt. 


Bereiche der Arbeitsspannung des Luftspitzenzählers 


Die Arbeitsspannungsbereiche des Luftspitzenzählers können hinsichtlich dessen 
Arbeitsweise in zwei Hauptbereiche eingeteilt werden: 

a) proportional- und Geiger-Bereich, 

b) Bereich der sogenannten Dauerentladung (koronierender Zähler). 

In dem proportional und Geiger-Bereich arbeitender Luftspitzenzähler ist im Stande, 
wie bekannt, die einzelnen «. f sowie die langsamen Ladungsträgern und Elektronen zu 
'egistrieren. 

Der koronierende Luftspitzenzähler kann aber nicht mehr zur Detektion der einzelnen 
-Teilchen, langsamer Ladungsträgern sowie Elektronen benutzt werden wegen der hohen 
etwa 10* Imp/Sek) Impulsgeschwindigkeit, die durch die in dem Zählvolumen herschende 
Dauerentladung verursacht ist. Es wird auch der Ausdruck benutzt, dass der koronierende 
Luftspitzenzühler auf die einzelnen Elektronen nicht mehr anspricht. Die koronierenden 
[onen-Zähler werden jedoch zur Detektion von z.B. x Teilchen auf einem hohen Elektronen 
Hintergrund benutzt. Sie werden in diesem Falle als Funken-Zähler (Spark counters) 
bezeichnet und waren bereits der Gegenstand mehrerer Untersuchungen. Wie es weiterhin 
pekannt ist, spricht so ein koronierender Zähler auch mit Leuchtfadenimpulse auf Wasser- 
lampf und wie letzte Untersuchungen zeigen auch auf mehrere andere Flüssigkeitsdämpfe an. 

Wenn man also über durch einen Flüssigkeitsdampf verursachte Impulse sprechen 
will, müssen die beiden, oben erwähnten Arbeitsspannungsbereiche des Luftspitzenzühlers 
‚charf voneinander unterschieden werden. Um den Missverständnissen aus dem Wege 
u gehen, wird die durch einen Flüssigkeitsdampf ausgelóste Impulszahl, weiterhin mit N 
yezeichnet, wenn der Zähler im proportional- oder Geiger-Bereich betrieben wird. Mit Ny, 
lagegen, wenn die Impulse in Form von Leuchtfadenimpulsen in dem koronierenden Zähler 
Funken-Zähler) ausgelöst werden, wie dieses schon im Falle der durch den Wasserdampf 
iusgelósten Impulse vorgeschlagen wurde (Sujak 1958b). 


Versuchsanordnung 


Die benutzte Versuchsanordnung bestand aus einem mit Gitter abgeschlossenem 
‚uftspitzenzähler, Glasbehälter mit der zur Untersuchung gelangter F lüssigkeit (Abb. 1.) 
md den elektronischen Kreisen, die zur Registration der Impulse dienten, zusammen. 
Der Zähler wurde mit einem dichten Cu-Gitter von etwa 900 Maschen pro.cm? von 


ler Seite der Flüssigkeit abgeschlossen. Er bestand aus einem Messingzylinder von etwa 
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1.8 cm Innendurchmesser und etwa 3,5 cm Länge. Der Zähler wurde mit einer scharfen 


Stahlspitze bestückt. Die Entfernung Spitzenende-Gitter betrug etwa lem. Ein Kühlungs- 


mantel mit zirkulierendem Wasser ermöglichte die Temperatur der Zählers konstant auf 


etwa 22°C zu halten (Hoepler-Thermostat). Ein Glasbehälter mit der zur Untersuchung 


gelangten Flüssigkeit konnte unter das begitterte Ende des Luftspitzenzählers unterschoben 


werden. Die Entfernung Gitter-Flüssigkeitsoberfläche betrug etwa 5 cm. Der Durchmesser 


des Glasbehälters war etwa 4 cm. 


Der elektronische Teil der Apparatur bestand wie üblich aus einem Hochspannungsgerät, 


Verstärker und Impulsmittelmesser mit Selbstschreibegerät. Zur Kontrolle und Beobachtun- 


gen der Impulse wurde dauernd ein Kathodenstrahloszillograph benutzt. 


Abb. 1. Versuchsanordnung: 


Flüssigkeit 
Wasser 
Benzol 
Toluol 
Methylalkohol 
Ä thylalkohol 
Essigsäure 80%, 
Azeton 
Amylazetat 


Chlorkohlenstoff 


Chloroform 
Trichloräthylen 


* Bei 16°C 
** Bei 31,4°C 


Chemische 


Formel 


H,O 
CH 
C,H, 
CH,O 
C,H,O 
C,H,0, 
C,H,0 
C;H,40, 
C 
CHCI, 
C,HCI, 


] — Luftspitzenzähler, 2 — Kühlwassermantel, 


4 — Flüssigkeit. 


TABELLE | 


3 — Flüssigkeitsglasbehälter, 


[Siedepunkt bei) Dielektrische | Dipolmoment 
760 mm Hg | Konstante bei |der Flüssigkeit 
°G | 1. GER GO OSS 
100 80,8 1,84 
80,1 2,28 0,08 
110,7 23 0,39 
51 912 1,68 
18,4 25 1,67 
118,1 9,7 0,83 
56,2 21,4 2,71 
141 4,75 1,91 
76,7 2,24 0 
61.1 4.81 1,18 
87,2 3,42* 0,9 


Druck der 


gesätigten 
Dämpfe bei 
20°C, mm Hg 


160,0 
100,0** 


CYA 
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Untersuchte Flüssigkeiten 


Die Flüssigkeiten, die zur Untersuchung auf die Einwirkung deren Dämpfe auf den 
Luftspitzenzähler gelangten, wurden in der Tabelle I zusammengestellt. 


Versuchsergebnisse 


I. Proportional- und Geiger-Bereich 


AZETON. Nachdem Azeton unter den Luftspitzenzähler unterschoben wird, werden 
in einigen Sekunden Arbeitsimpulse im Zähler ausgelöst. Die Impulsgeschwindigkeit steigt 
mit der Zeit um über ein Maximum zu gehen und wiederum bis zu den Nullwerten abzufallen. 
Die ausgelöste Impulsgeschwindigkeit laufend auf dem Papierband des Selbstschreibege- 
rätes registrierte, ist auf der Abb. 2 als eine typische Zeit-Kurve für Azeton wiedergegeben. 
In diesem Falle wurde der Zähler im Proportional-Bereich bei einer Arbeitsspannung von 
etwa 3120 V betrieben. 

Zur gleichen Zeit mit der Aufnahme der Zeit-Kurve wurden die Änderungen der 
einzelnen Impulse mit Hilfe eines Kathodenoszillographen beobachtet. Die Gestalten der 
einzelnen Azeion-Impulse, die bei gegebener Zeit der Beobachtung auf dem Schirm des 
Kathodenstrahloszillographen erschienen, wurden schematisch auf Abb. 2. auch wiedergege- 
ben. Die Zeit-Momente, denen die entsprechende Gestalten der Impulse zuzuschreiben 
sind, wurden mit den Buchstaben A, B, C oberhalb der Zeit-Kurve markiert. 


A B 5 - 
A AB OC GO C/8A. A 
| all | 
[ on I ^ 
Luzeton -l 
x : Í r 
kH f 
0 
t 180 120 60 0 


sec 

Abb. 2. Zeit-Kurve fiir Azeton. Arbeitsspannung des Zahlers V 4 aus dem Proportional-Bereich etwa 3120 V. 

Oberhalb der Zeit-Kurve sind Zeitpunkte durch A, B, C, angedeutet, denen in Kreisen schematisch angegebene 

auf dem Oszillographen-Schirm beobachtete Impulsgestalten entsprechen. Durch y ist die Zeit markiert während 

welcher der Zähler y-Strahlen registrierte. Durch „Azeton“ ist die Zeit bezeichnet während welcher Azeton 
unter dem Zähler unterstellt war. Weitere Einzelheiten im Texte. 
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Nachdem der Behälter mit Azeton unter den Zähler unterschoben wurde (t, auf Abb. 2.). 


änderte sich mit der Zeit der Beobachtung die Amplitude und Gestalt der Impulse so als 


ob das elektrische Feld in der Spitzenumgebung wachsen würde. Die dabei am Zähler 
anliegende Arbeitsspannung bleibte aber unverändert. Nach etwa 20 Sekunden (bei der 
benutzten Versuchsgeometrie) wurden die Impulse „Kammartig“ gefunden (Zeit-Mo- 
ment C) und die Impulsgeschwindigkeit sank mit der Zeit bis zum Nullwert obwohl auch 


der Behälter mit Azeton in der unveränderten Lage bleibte. Dieses führte zur Ausbildung. 


des ersten Maximum auf der Azeton Zeit-Kurve. Nachdem die Impulsgeschwindigkeit - 


bis zum Nullwert sank, wurde der Behälter mit Azeton zur Seite geschoben. Dieser Zeit- 
Moment wurde auf der Zeit-Kurve mit t, bezeichnet (Abb. 2.). Wie zu sehen ist, steigte 
danach die Impulsgeschwindigkeit wieder, laüfte über ein zweites Maximum um wieder 
zum Nullwert abzufallen. Nachdem das zweite Maximum überschritten wurde, rückte auch 
die Amplitude und Gestalt der Einzelnen Impulse zu den Anfangswerten zurück (Zeit- 
Momente A). Dieses konnte beliebig oft wiederholt werden. 

Schon auf dieser Stelle muss darauf aufmerksam gemacht werden, dass sowohl die 


Zeit-Momente bei welchen die erwähnten Maxima auf der Azeton Zeit-Kurve erreicht 


werden, wie auch der Verlauf der Kurven selbst, sehr von der am Zühler angelegter Arbeits- 
spannung abhängten. Diese Arbeitsspannungsabhängigkeit wird am Ende des Berichtes 
noch einmall allgemein besprochen. 

AMYLAZETAT. Sehr ähnliches Verhalten zeigten die, durch Amylazetatdampf aus- 
gelösten Impulse, nur waren die Verhältnisse nicht so ausgeprägt wie das bei Azeton der 
Fall gewesen ist. 


Trienleräthylen 


320V — t it 
| 
————— E 0 
t. 180 120 60 0 


Abb. 3. Zeit-Kurve für Trichloräthylen. Arbeitsspannung des Zählers V 4 aus dem Geiger-Bereich etwa 3220 V. 


Bei allen Abbildungen wird durch ¢, der Zeitmoment der Unterstellung des Behälters mit der Flüssigkeit be- 
zeichnet, dagegen mit ty der Zeitmoment der Abstellung des Behälters mit der Flüssigkeit markiert wird. 


TRICHLORATHYLEN. Die Amplitude sowie die Impulsgeschwindigkeit wächste 
mit der Zeit der Einwirkung des Trichloräthylendampfes bis zu einer Sättigung auf. Eine 
typische Zeit-Kurve für diese Flüssigkeit ist auf der Abb. 3. wiedergegeben. Dabei ist es 
zu bemerken, dass nachdem der Behälter mit dem Trichlorüthylen zur Seite geschoben 
wurde, zeichnete sich auf der Zeit-Kurve ein kleines Maximum an und die Impulsgeschwin- 
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ligkeit sank sehr schnell zum Nullwert zurück. Kurz vor dem vollständigen Ausbleiben 
ler Impulse ‚konnte das Kleinerwerden der Amplitude klar bemerkt werden. 
ESSIGSAURE-80%. Völlig anderes Verhalten der ausgelösten Impulse wurde für 
ssigsäure gefunden. Die Zeit-Kurve steigte mit der Zeit bis zu einer Sättigung an nach 
lem der Behälter mit Essigsäure unter den Zähler unterstellt wurde. Dies zeigt die Abb. 4. 


Essigsäure i E i 


⁄ = 3150 V 


cz 


360 300 240 180 120 60 0 


‚bb. 4. Zeit-Kurve für Essigsaüre. Arbeitsspannung des Zählers V 4 aus dem Proportional-Bereich etwa 3150 V. 


\uf der angegebenen Zeit-Kurve sind scharfe Knicke zu sehen. Die Gestalt der Impulse 
nd deren Amplitude änderten sich nicht so ausgeprägt mit der Zeit wie das bei Azeton, 
unylazetat oder sogar Trichloräthylen der Fall war. Die ausgelöste Impulsgeschwindigkeiten 
üngten aber sehr stark von der am Zähler angelegten Arbeitsspannungen ab. 

TOLUOL. Toluol zeigte eine der Essigsäure ähnliche Wirkung auf den im proportional- 
der Geiger-Bereich betriebenen Luftspitzenzähler auf. Die hohe Impulsgeschwindigkeit, 
ie durch Toluoldämpfe ausgelöst wird muss wohl einem übergehen des Zählers in die 
Jauerentladung zugeschrieben werden. Dieses hatte dann zur Folge, dass nachdem, der 
ehálter mit Toluol abgestellt wurde, die Entladung reiste auf einmal scharf ab. Eine typische 
eit-Kurve für Toluoldampf ist auf der Abb. 5. wiedergegeben worden. 

CHLOROFORM. Die ausgelöste Impulsgeschwindigkeiten waren so gross, dass bei den 
enutzten Versuchsbedingungen konnten kaum übersichtliche Messergebnisse erhalten 
erden. 

BENZOL. Benzol bewirkte keine Impulse im Luftspitzenzähler wenn dieser im Pro- 
ortional- oder Geiger-Bereich betrieben wurde. 

CHLORKOHLENSTOFF. Ebenso bewirkte Chlorkohlenstoff keine Impulse in dem 


uftspitzenzähler wenn dieser im Proportional- oder Geiger-Bereich betrieben wurde. 


| Bereich der Dauerentladung (Koronierender Zähler) 


Wenn an den Luftspitzenzähler eine Arbeitsspannung von etwa über 3250 V angelegt 
ird, so laüft er in die sogennante Dauerentladung über. In diesem Zustand spricht der 
ähler auf mehrere Flüssigkeiten mit Leuchtfadenimpulse als Funken-Zähler an. Wenn ein 
ıtsprechend gewähltes Diskriminationspotential am Amplituden-Diskriminator gelegt 
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wird, das dafür sorgt, dass die viel kleineren Townsend’schen Korona-Lawinenimpulse: 
nicht mehr registriert werden, so können auch Zeit-Kurven für Leuchtfadenimpulse auf-. 
genommen werden. Diese Zeit-Kurven werden aber folgend bezeichnet: 


(Fr), =f a) 


t 


Dabei ist wie vorher gemeint, dass alle Parameter bis auf die Flüssigkeit (£,) konstant 
gehalten sind. 

Einige der typischen Zeit-Kurven der Leuchtfadenimpulse wurden auf den Abb. 6. 
bis Abb. 10. wiedergegeben und zwar für: Azeton — Abb. 6., Methylalkohol — Abb. 7 


Toluol 
V4 = 3150 V 


T EFE 


Abb. 5. Zeit-Kurve für Toluol. Arbeitsspannung des Zählers V4 aus dem Proportional-Bereich etwa 3150 V. 


Methylalkohol 
Vy = 3710 V 


tQ 180 120 60 0 BERN 120 50 0 
Abb. 6. Abb. 7. 
Abb. 6. Zeit-Kurve für Azeton. Arbeitsspannung des Zühlers V 4 aus dem Bereich der Dauerentladung etwa 3710 V. 


Abb. 7. Zeit-Kurve für Methylalkohol. Arbeitsspannung ey Zählers V 4 aus dem Bereich der Dauerentladung 


etwa 3710 V, 


Athylalkohol 
VIS ad E Amylazetat oy 
5 E y = 37104 
t? t 
Pa SET UE TT on —— 0 
Er 180 120 60 0 
Abb. 8. NSH, ©: 


Abb. 8. Zeit-Kurve für Äthylalkohol. Arbeitsspannung des Zählers V4 aus dem Bereich der Dauerentladung 
etwa 3710 V. 
Abb. 9. Zeit-Kurve für Amylazetat. Arbeitsspannung des Zählers V4 aus dem Bereich der Dauerentladung 


etwa 3710 V. 
t; t 
| 
Nw 
Wasse SW 
m nau i 
i 
= n 0 


300 240 180 120 60 a 
sec 


Abb. 10. Zeit-Kurve für Wasser Arbeitsspannung des Zählers V4 aus dem Bereich der Dauerentladung 
etwa 3710 V. 


Äthylalkohol — Abb. 8., Amylazetat — Abb. 9., Wasser — Abb. 10. Alle angeführte 
Zeit-Kurven wurden bei einer gleichen Arbeitsspannung des Luftspitzenzählers von etwa 
3710 V aufgenommen. 

Zu den Zeit-Kurven der Leuchtfadenimpulse ist folgendes zu bemerken. Nachdem 
las erste Maximum auf der Zeit-Kurve für Azeton überschritten ist, zeigt diese Zeit-Kurve 
inem unstätigen Charakter auf. Diese Unstütigkeiten beruhen auf kurzzeitigem Auslóschen 
ler Dauerentladung im Zähler, was mit Hilfe des Kathodenstrahloszillographen zu folgen 
st. Nachdem der Behälter mit Azeton zur Seite geschoben wird, zeigt sich auf der Zeit- 
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A B 
e = 
Sn Chiorkohlenstoff t 
lrichlorathyLen 
Va = 3570 V 


180 120 60 0 t E 180 120 60 0 
Abb. 11. Abb. 12. 
Abb. 11. Zeit-Kurve für Trichloräthylen. Arbeitsspannung des Zählers V 4 aus dem Bereich der Dauerentladung 
etwa 3570 V. Während der Zeit, die durch B oberhalb der Zeit-Kurve angedeutet wurde, konnte auf dem 
Leuchtschirm des Oszilloscopes das Bild B (im Kreis), und während der Zeit A das Bild 4 (im Kreis) beobachtet 
werden. 

Abb. 12. Zeit-Kurve für Chlorkohlenstoff. Arbeitsspannung des Zühlers aus dem Bereich der Dauerentladung 
etwa 3570 V. 


sec 


t 420 360 300 240 180 120 60 0 


Abb. 13. Zeit-Kurve fiir Benzol. Arbeitsspannung des Zählers aus dem Bereich der Dauerentladung etwa 3710 V, 


Kurve fiir Azeton und Methylalkohol ein scharfes Maximum an, dem ein scharfer Abfall 
bis zu dem Nullwert folgt. Dieses Maximum kommt noch ausdrücklicher zum Vorschein 
bei den Trichloräthylen Zeit-Kurven (Abb. 11.). Auch für Chlorkohlenstoff zeigt die Zeit- 
Kurve so ein Maximum (Abb. 12.). Bei der Trichloräthylen Zeit-Kurve (Abb. 11.) ist 
das Entstehen des ersten und des zweiten Maximum auf das Kleinerwerden der Amplituden 
der einzelnen Impulse (etwa 3 mal kleiner), die einfach nicht mehr registriert werden, 
zurückzuführen. Die wahre Impulsgeschwindigkeit, wie das mit Hilfe des Kathodenstrahl- 
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t 120 60 0 


u 180 120 60 0 ER 240 180 120 60 0 


V = 4430 V 


7er 
BR 2) AQ Ast) ye) “pO oto) 300 


e 


hlers V4 aus dem Bereich der Dauerentladung auf den Verlauf 


Abb. 14. Einfluss der Arbeitsspannung des Zà 
) V 47-4150 V, e) V 4—4430 V. 


der Zeit-Kurve für Azeton: a) V 43290 V, b) V 4=3570 V, c) V4=3860 V, d 
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oszillographen zu sehen ist, bleibt dabei sehr gross. Bei der Zeit-Kurve für Chlorkohlensto 
(Abb. 12.) ist dagegen der rasche Abfall der Impulsgeschwindigkeit bis zum Nullwert in 
den ersten Sekunden nach dem der Flüssigkeitsbehälter unterschoben wurde, auf den 
realen Abfall der Zahl der ausgelösten Impulse im Zähler zurückzuführen. Die Amplitude 
der Impulse ändert sich dabei nur wenig. 

Zu der Zeit-Kurve für Benzol (Abb. 13) ist noch hinzuzufügen, dass Benzol nur sehr 
kleine Geschwindigkeiten der Leuchtfadenimpulse bei Arbeitsspannungen des Zählers 


E 


P =F) 


— 
L2 i 
fe) a 
x e Benzol 
4 © Athylalkohol 
X Azeton 
A Methylalkoho! 


3 NES 36 40 44 48 57 


Abb. 15. Die durch einige Dämpfe ausgelöste Impulsgeschwindigkeit als Funktion der Arbeitsspannung V 4 des 
sich im koronierenden Zustand befindenden Zühlers: @ — Benzol, O — Athylalkohol, x — Azeton, A — 


Met hylalkohol. 


bei welchen andere Flüssigkeiten sehr grosse Leuchtfadenimpulsgeschwindigkeiten zeigten 
(siehe z.B. Abb. 15.), bewirkte. Es muss noch weiterhin bemerkt werden, dass die Benzol- 
dämpfe das Grösserwerden der Amplituden der Dauerentladungsimpulse mit der Zeit 
der Einwirkung bewirken. Die Amplitude wird mit der Zeit der Einwirkung etwa 2 mal 
grösser als im Falle, wo keine Benzoldämpfe in das Zählvolumen des Zählers eindringen 
können. Dieses Verhalten ist jedoch von der Leuchtfadenbildung verschieden. Es ist mehr 
als eine Auswirkung einer Erniedrigung des elektrischen Feldes in der Spitzenumgebung 
zu betrachten. 

Toluen bewirkte nur vereinzelte Leuchtfadenimpulse. 

Essigsäure bewirkte dagegen überhaupt keine Leuchtfadenimpulse für die benutzte 
Messanordnung und Arbeitsspannungen bis sogar etwa 5000 Volt. Die Amplitude der 
einzelnen Korona-Impulse ist immer während der Einwirkung von Essigsäure kleiner 
gefunden worden, im Gegenteil zum Benzol-Fall, wo die Amplitude der einzelnen Korona- 
Impulse durch den Benzoldampf vergrössert wurde. 

Am Beispiel des Azetons wird der Einfluss der Arbeitsspannung des Zählers auf die 
Zeit-Kurven der Leuchtfadenimpulse illustriert. Auf der Abb. 14, sind Zeit-Kurven gegeben, 
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die für verschiedene Arbeitsspannungen des koronierenden Spitzenzählers aufgenommen 
wurden und zwar: Kurve 1 — 3290 V, Kurve 2 — 3570 V, Kurve 3 — 3860 V, Kurve 4 — 
4150 V, Kurve 5 — 4430 V. Sehr ähnliche Einfluss der Arbeitsspannung des Zählers auf 
die Zeit-Kurven wurde für Methylalkohol gefunden. In dem letzten Falle fehlten nur die 
für das Azeton charakteristische Unregelmässigkeiten. Da die Untersuchungen des Einflusses 
der Arbeitsspannung des Zählers auf das Ansprechvermögen auf verschiedene Dümpfe 
nur einen Einführenden Charakter trugen sind nur einige Ergebnisse auf der Abb. 15. 
wiedergegeben. 

Die Messreihen lieferten oft nicht ganz reproduzierbare Ergebnisse, besonders wenn 
es sich um die Spannungsabhängigkeiten handelte. Die Messkurven zeigten oft Histeresis- 
Effekte, die sehr von der Zeit, die zur Aufnahme der Kurven benötigt wurde, abhängten. 
Es muss noch darauf hingewiesen werden, dass für den Benzoldampf erst bei höheren 
Arbeitsspannungen von etwa 4,5 kV grössere Impulsgeschwindigkeiten der Leuchtfadenim- 
pulse gemessen wurden (Abb. 15.). 


Abschliessende Bemerkungen 


Wie aus den vorgelegten Messreihen hervorgeht, können die Dämpfe der polaren wie 
auch der nichtpolaren Flüssigkeiten im offenen Luftspitzenzähler Impulse auslösen, je 
nachdem, ob der Zähler im Proportional- und Geiger-Bereich oder im Dauerentladung- 
Bereich betrieben wird. 

Wenn der Zähler im Proportional- oder Geiger-Bereich betrieben wird, so werden 
Impulse im Zähler nur im Falle der polaren Flüssigkeiten ausgelöst. Die nichtpolaren 
Flüssigkeiten dagegen wie z.B. die untersuchten Benzol und Chlorkohlenstoff lösen in diesem 
Falle keine Impulse aus. 

Wenn der Zähler im Bereich der Dauerentladung betrieben wird so wirkt auch Benzol 
und Chlorkohlenstoff auf den Zähler ein. In diesem Falle können auch Zeit-Kurven für 
die ausgelösten Impulse aufgenommen werden. 

Die Amplitude der ausgelösten Impulse ändert sich mit der Zeit der Einwirkung der 
Flüssigkeiten, was scheint mit den Änderungen der Dampfkonzentration in der Zählspit- 


zenumgebung zusammenzuhängen. 
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